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ACOUSTIC  SCATTERING  CROSS-SECTIONS  FOR 
TRULY  COMPOSITE  WIND-WAVE  SURFACES: 
SCATTERING  WITHOUT  BUBBLES* 

by 

David  Middleton** 


1.  Introduction 

The  study  of  scattering  from  random  rough  surfaces  and  in  particu¬ 
lar  here,  acoustic  scattering  from  moving  random,  wind-generated  wave 
surfaces,  is  a  venerable  subject  [l]-[4],  [6],  [7]-[l3],  which  is  of 
continuing  interest  and  importance.  This  is  true  scientifically,  be¬ 
cause  of  both  analytical  and  physical  problems.  The  former  center 
around  methods  of  approximation  needed  for  explicit  theoretical  results; 
the  latter,  in  choosing  physical  models  appropriate  to  the  wave-surface 
in  question.  Scattering  from  such  surfaces  is  also  of  critical  concern 
in  applications,  where  reverberation  is  often  a  controlling  factor, 
including  under-ice,  ocean  surface,  bottom  and  similar  active  underwater 
acoustic  environments.  Analogous  problems  and  applications  arise  in  the 
scattering  of  electromagnetic  waves  from  the  atmosphere-ocean  interface. 

Here  we  are  concerned  primarily  with  determining  the  intensity  of 
underwater  acoustic  radiation  scattered  from  random  moving  wave  surfaces. 
In  particular,  our  first  aim  is  to  derive  explicit  relations  for  the 
coherent  and  incoherent  (acoustic)  scattering  cross-sections  for  such 
surfaces,  by  somewhat  different  procedures  and  assumptions  from  those 
used  in  previous  work.  Both  mono-  and  bi-static  cross-sections  are 
included. 

Our  second  goal  is  to  construct  a  surface  wave  model  which,  when 
applied  to  our  cross-sectional  results,  has  the  potential  for  explaining 


*Work  supported  under  contracts  with  Naval  Underwater  Systems  Center, 
New  London,  Conn.,  1983,  Contract  N00140-83-M-NA11  (1983),  and  in  part 
under  Contracts  N00140-83-C-KA23,  and  N00140-84-M-LZ39  (1983,  1984). 

♦♦Contractor,  Physics  and  Applied  Mathematics,  127  E.  91  St.,  New 
York,  NY  10128. 
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the  discrepancies  between  theory  and  experiment  in  the  important  regimes  in¬ 
volving  backscattering  at  high  frequencies  and  small  grazing  angles  [l]-[4], 
without  recourse  to  a  bubble  mechanism  [l]-[4]. 

Here  the  dominant  scatter  contribution  is  no  longer  provided  by  the 
specular-point  (or  "facet")  scatter  [5],  but  by  the  small-scale  surface  com¬ 
ponents  (e.g.,  Bragg  scatter),  which  now  strongly  dominates  the  former,  pro¬ 
vided,  of  course,  that  there  is  an  adequate  small-scale  component  in  the  first 
place.  This  second  aim  is  stimulated  by  the  fact  that  backscatter  data  ex¬ 
hibiting  the  above  discrepancies  can  occur  without  bubbles  [6],  [7].  It  is 
also  encouraged  by  the  additional  observation  that  in  most  of  the  earlier  re¬ 
ported  data  [1],  [3],  [4]  there  appears  to  be  no  quantitative  verification 
of  the  presence  (or  absence)  of  near-surface  bubble  layers,  sufficient  to 
account  for  the  noted  discrepancy.  In  addition,  the  empirical  scatter  data 
vary  with  (input  signal)  frequency  in  a  way  not  adequately  predicted  by  ear¬ 
lier  theoretical  models.  Furthermore,  comparisons  of  (above-surface)  radar 
and  (below-surface)  acoustic  (backscatter)  data  [4]  are  not  really  convincing, 
because  they  were  obtained  for  different  oceans  (at  different  times),  and  with 


unknown  or  inadequate  "ground-truth":  for  example,  without  simultaneous  deter¬ 


mination  of  bubble  densities  below  the  surface  and  above-surface  water  droplet 
populations,  which  latter  can  noticeably  affect  refraction  and  the  scatter 
intensity. 

In  short,  we  consider  it  reasonable  to  state  that,  at  the  least,  bubble 
mechanisms  are  not  necessarily  the  explanation  of  the  aforementioned  dis¬ 
crepancies  between  theory  and  observation  (particularly  in  the  absence  of 
bubbles!).  In  fact,  there  is  evidence  for  a  wave  surface  mechanism  which  can 
possibly  account  for  the  observed  scatter  levels  (at  high  frequencies  and 
small  grazing  angles),  within  the  generally  available  theoretical  methods 
employed  earlier,  and  with  the  generalizations  employed  here.  The  postulated 
scattering  mechanism  (suggested  by  R.  H.  Mellen  [14])  is  the  ensemble  of  soli- 
ton-ripples  [cf.  Sec.  2.1  ff . ] .  These  are  generated  by  the  near-surface  wind 
impinging  on  the  (irregular)  gravity-capillary  wave  surface,  upon  which  ride 
the  soli  ton-ripples.  (These  soli  tons  are  found  to  be  the  limiting  solution 
of  the  Korteweg-de  Vries  equation  and  are  discussed  in  some  detail  by  Light- 
hill  [16],  see  also  [47].)  "Solitons,"  here  "hydraulic  jumps,"  are,  among 
other  things,  characterized  by  the  fact  that  they  preserve  their  shape  and 
speed  [16];  [47];  Sec.  III.  Preliminary  investigations  [17]-[21]  suggest 
that  these  hydraulic  jumps  not  only  embody  the  non-linear  mechanisms 
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whereby  the  wave  surface  itself  can  be  generated  by  the  wind,  but 
may  also  be  sufficiently  large  to  produce  the  observed  scattering 
levels.  What  we  have  here  is  a  truly  two-component  surface:  the  soli  ton 
ripples  ride  upon  the  gravity-capillary  wave  surface  which  is  their 
(quasi )-equilibrium  result.  Unlike  the  latter,  these  ripples*  are  non- 
dispersive  and  essentially  move  down-wind  locally,  and  travel  with  the  same 
velocity  o(0.4  m/sec,  [15],  [47]).  For  acoustic  and  electromagnetic  analyses 
of  two-  (or  more-)  component  surfaces,  see  [2],  [4],  [21]- [26],  [29]. 

The  usual  difficulties  with  the  theory  arise  in  the  inevitable 
approximations:  the  Tangent-Plane  or  Kirchoff  approach  used  for  the 
(moving)  large-scale  wave  surface  components  [34],  and  perturbations 

[1] -[3],  [8]-[ll],  [29],  or  other  modifications  (as  used  by  the  author 
here,  cf.  Sec.  7  ff.)  for  the  small-scale  surface  elements.  However, 
the  resulting  discrepancies  can  be  shown  to  be  quite  small  o(  1-2 . 5  db) 
vis-a-vis  the  uncorrected  Kirchoff  method  ([4],  for  example).  Diffraction 
terms  are  usually  neglected  in  the  analysis,  but  can  also  be  shown  to  be 
negligible  in  most  cases  with  respect  to  the  principal,  direct-scatter 
contributions  [cf.  Sec.  3.3  ff.].  As  we  shall  see  in  the  text  following, 
the  various  current  theoretical  procedures  yield  analytical  results  which 
are  formally  very  similar.  However,  they  can  diverge  numerically  because 
of  different  choices  of  the  physical  scatter  mechanism:  capillary  waves 

do  not  provide  sufficient  cross-section,  for  instance,  to  account  for 
the  high-frequency,  low-grazing  angle  discrepancies  in  (backscatter) 
cross  sections,  which  have  been  observed  empirically.  Their  (Bragg- 
scatter)  contributions  are  o( 10-20  db)  too  small.  Thus,  we  shall  con¬ 
clude  that  the  observed  difficulties  stem  not  from  fundamental  inade¬ 
quacies  of  theory  and  approximation  thereto,  but  rather  from  the  selection 
of  a  suitable  physical  model.  As  noted  above,  our  proposed  solution 
is  the  soli  ton-ripple  mechanism  for  the  principal  Bragg  scatter  contri¬ 
butions,  without  the  intervention  of  a  bubble  mechanism. 

Our  analytical  innovations  here  are:  (1)  the  use  of  a  physically 
independent  scatter  surface  (the  solution-ripples)  riding  on  the  gravity- 
capillary  wave  surface  which  in  turn  is  generated  by  the  energy  trans¬ 
ferred  nonlinearly  via  these  ripples  [14]-[12],  with  a  non-vanishing 
"til t-factor"--produced  by  the  large-scale  surface--now  determined  by: 

(2)  the  Kirchoff  approximation  applied  to  the  development  of  the  Bragg 


*0ften  called  "cat's  paw"  in  the  initial  stages  of  wind-wave  surface 
excitation,  [15]. 


scatter  terms  (instead  of  a  perturbational  approach  [4]).  This  removes 
the  restrictions  of  small  rms  heights  and  slopes  required  in  the  per¬ 
turbation  procedures  when  applied  to  the  small-scale  components,  and 
permits  calculation  of  higher  order  Bragg-scatter  terms  if  needed  [24]. 
(3)  We  treat  the  gravity-capillary  wave  surface  component  as  the  single 
surface  it  is,  and  use  the  approximate  but  justifiable  technique  of 
splitting  this  surface  into  low-  and  high-frequency  components.  (The 
associated  "tilt- factor"  here,  however,  vanishes  with  vanishing  grazing 
angle  (eoT  -*•  ir/2.)  An  insightful  justification  of  this  procedure, 
from  a  more  general  theoretical  viewpoint,  has  been  given  recently  by 
Bahar,  Barrick,  and  Fitzwalter  [24],  [25],  in  the  analogous  but  more 
complex  situation  of  electromagnetic  scattering  from  (simple,  i.e., 
single  component)  wave  surfaces.  (4)  We  develop,  and  employ,  a  corre¬ 
lation  function  and  corresponding  intensity  wave  number  spectrum  for  the 
soliton  waves,  obtaining  estimates  which  indicate  the  potentiality  of  the 
new  mechanisms  here. 

This  Report  is  organized  as  follows:  Section  2  provides  a  de¬ 
tailed  analytical  summary  of  the  mono-  and  bi-static  scattering  cross 
sections,  for  high  and  low  frequencies  (i.e.,  large  and  small  Rayleigh 
numbers),  and  includes  comparisons  with  the  results  of  Tolstoy  and 
Clay  [9],  McDaniel  and  Gorman  [4],  Kur'yanov  [21],  [1],  and  Bahar  et  al. 
[25].  In  Section  3  some  preliminary  numerical  results  are  obtained, 
based  on  the  soliton  wave  model  described  in  Sec.  2.1.  Included  also 
are  comparisons  with  backscatter  from  capillary  waves  and  a  short 
account  of  the  potential  contributions  from  diffraction  effects,  which 
are  seen  to  be  negligible.  Section  4  completes  Part  I  with  a  concise 
review  of  the  principal  new  results  and  some  next  steps  in  the  experi¬ 
mental  and  theoretical  analysis  of  this  class  of  scattering  problems. 

Part  II  consists  of  Sections  5-8,  where  the  mathematical  details  of 
this  study  are  developed; (see  the  Table  of  Contents). 

Finally,  we  emphasize  that  bubble  mechanisms  can  (and  in  some 
instances,  do)  account  for  the  observed  larger  (back-)  scattering 
cross  sections  observed  at  high  frequencies  and  small  grazing  angles. 
However,  here  we  offer  a  competing  mechanism,  in  the  form  of  soliton 
ripples  [11]-[19] .  Preliminary  observations  and  the  present  analysis 
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indicate  that  these  ripples  can  potentially  account  for  the  observed 
data  when  bubbles  are  absent  [6],  [7].  Further  experimental  results 
are  needed  to  establish  this  possibility  conclusively,  but  already  the 
work  of  Roderick  ([6],  and  Table  3.2  ff.)  indicates  encouraging 
potential  support  for  this  mechanism. 

Part  I.  Acoustic  Scattering  Cross  Sections 

Scattering  Cross  .Sections  for  Composite  Wave  Surfaces 
As  noted  above  (and  subsequently  in  Sections  7,  8)  the  principal 
aim  of  this  initial  study  is  to  obtain  analytic  expressions  for  the 
acoustic  scattering  cross  sections  of  a  truly  two-component  or  composite 
moving  wave  surface.  This  includes  an  initial  numerical  demonstration 
that  values  comparable  to  those  observed  experimentally  can  be  obtained 
with  the  introduction  of  a  suitable  physical  model,  without  having  to 
postulate  a  bubble-scattering  mechanism  [3],  [4].  This  is  particularly 
important  in  situations  where  no  (or  negligible)  bubbles  are  found  to 
be  present  [6]. 


2.1  The  New  Wave-Surface  Model:  Surface  Elevation 

We  postulate  the  following  potential  mechanism  for  the  small-scale 
scattering  component  of  a  typical  wind-wave  surface:  ensembles  of  soli  tons, 
or  hydraulic  jumps,  produced,  as  shock  waves,  by  the  impact  of  the  local  wind 
on  the  gravity-capillary  wave  structure.  These  solutions,  which  mathe¬ 
matically  are  the  limiting  solution  of  the  Korteweg-de  Vries  equation 
for  cnoidal  waves  ([16],  pp.  465-467),  travel  in  all  directions  on  the 
rough,  large-scale  surface,  with  constant  speed.  They  travel  on  a  thin 
viscous  layer  and  probably  represent  the  principal  (and  nonlinear) 
mechanism  whereby  wind  energy  is  transferred  to  the  overall  wave  sur¬ 
face.  On  initially  still  surfaces  they  appear  as  the  familiar  "cats-paw" 
effect  and  have  rather  many  ripples  in  a  typical  "wave"  packet.  As  the 
surface  builds  a  wave  structure,  with  the  transfer  of  energy  from  these 
soliton  trains  progressively  into  the  small  (capillary  and  then  gravity) 
wave  numbers,  the  wave  surface  becomes  quite  irregular  and  the  soliton 
packets  lose  their  directionality  and  structure,  becoming  essentially 
a  collection  of  individual  soli  tons  traveling  over  the  new  rough  gravity- 


capillary  surface.  (See,  for  example  the  illustration  in  [26]  opposite 
the  inside  book  title  page.)  Here  we  are  concerned  with  the  amplitudes 
and  durations  of  these  soli  tons  (as  discussed  in  Sec.  3.1  ff.)  from  the 
present  viewpoint  of  acoustical  scattering,  rather  than  with  their 
oceanographic  implications  for  wave  generation.  (See  Mellen  [14],  [15 3; 
also  [17 ],  [20],  and  in  particular,  Lighthill  [163;  also,  generally,  [47].) 

Unlike  the  soli  ton-ripple  component,  which  is  nondispersi ve,  rides 
upon,  and  is  essentially  independent  of  the  underlying  gravity-capillary 
wave  surface,  this  latter  (or  "G"-surface)  is  a  single,  directional 
surface,  with  a  wave  structure  obeying  the  dispersion  law  u>s  = 

[qK.  +  (3/Pw)K^f2,  with  K  =  2*/ A  where  f  (=  <*>  /2n)  is  the  frequency 
of  a  typical  fourier  component  of  the  (moving)  G-wave  surface  [26]- [28] 
and  3  is  the  surface  tension  (force  x  distance)  coefficient  of  water, 
with  pw  the  water  density.  It  is  convenient  to  divide  this  wave  surface 
spectrally  into  a  large-scale  or  “gravity-wave"  component  (g)  and  an 
essentially  independent  small-scale  "capillary"  component  (c),  at  some 
representative  wave  number,  kg.  This,  of  course,  is  a  mathematical 
device,  since  there  is  no  abrupt  transition  between  the  gravity  wave 
number  domain  and  the  capillary  domain.  This  surface  is  really  a 
single  rough  surface,  where  the  "capillary"  component  consequently  does 
not  ride  upon  the  gravity  component. 

With  the  above  in  mind  we  can  write  for  the  wave  surface  elevation 


iK’t)  =  lz  CcgCr.t)  +  cc(r,t)]  +  nG(r,t)t$(r,t),  r  =  jxx  + 


where  G  (=  g+c)  denotes  the  single  capillary-gravity  wave  surface  and 
S  indicates  the  elevation  of  the  soli  ton-ripples,  which  ride  on  i^Cg 
+?c)»  cf.  Fig.  2.1.  Here  nQ  is  the  normal  [at  ( r., t ) 3  to  the 
G-surface,  viz.. 


"6  =  «<ix'x+iyV  •]z)/<1  +  tx  *  CJ)I>,G.S+C  t  V  U» 


UL  e 
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with  generally 


n  =  (1  +  4?  +  .*-n  =  (Ltv  +  -  i,)n, , 


w 


(<C>  -  <5sXng>  =  <Cg>  =  0) 


Cross  section  of  a  moving  wind-wave  surface,  showing  soliton  ripples 
(greatly  exaggerated)  riding  on  the  gravity-capillary  wave 
surface  zr. 


since  the  rms  capillary  elevation  and  slopes  of  the  gravity  components 
are  small,  cf.  Eqs.  (7. 18)-  ( 7. 18a).  Here  c  and  cc  are,  respectively, 
the  elevation  of  the  "gravity"  and  "capillary"  components,  in  our 
artificial  division  of  the  single  wave  surface  Cg.  The  equilibrium 
surface*  is  <c>  3  0. 

The  spectral  division  of  the  G-surface  is  formally  indicated  by  the 
(surface)  wave-number-time  (intensity)  spectrum  W(kjx),  here  with  T  =  0. 

WG(k|°)  =  W  (k|0),  0  <  |k|  <  Ikpl, 

(2.3) 

=  Wc(k|0),  Ikjjl  <  |k|. 

The  two  components  of  Wq  have  a  common  directional  character,  determined 
by 


WQ(k  |0)  =  i  /  WG(fs)<5[k/2ir  -  ,Ks(a)/27T]X  df$,  f$  >  0, 


(2.4a) 


♦Since  the  solitons,  or  hydraulic  jumps,  are  one-sided  [16 ],  we 
see  that  <X>  =  <hQX?5>,  <s$>  >  0:  the  equilibrium  surface  is  <cq>  =  0. 
For  practical  purposes,  <ts>  <<  <C^>S,  so  that  we  may  regard  <£>  -  0 
here,  as  we  shall  do  henceforth. 


SsSwwV*.1 
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L  Vfs)<4  (i-4(“)/2-)>adfs 


+  /  wc(fs)<  *(*  -  JSs(«)/2^)>a  df$. 

fsD 


(2.4b 


where  =  Jc/2tt  and  JC$ (ot)  =  ^s(^sl2w*  ™  which  jw  is  the  unit  vector  for 

wave  components  of  frequency  fe,  e.g.,  i  *  1  cos  o  +  i  sin  a,  and  a 

5  ***x  -y 

obeys  a  probability  density  w^(a)  t  0,  aQ  -  tt/2  <  a  <  aQ  +  tt/2;  =  0 
elsewhere;  (there  are  no  g,c-waves  against  the  mean  wind  in  this  model). 
In  fact,  from  D,  Sec.  7.3,  we  have  explicitly 


<s<*  -  ■  I  5<fs  -  <2-5> 


which  when  applied  to  (2.4b)  gives  the  desired  wave  number  intensity 
spectrum  of  the  G-surface,  in  terms  of  its  two  (artificially)  desig¬ 
nated  components  (g,c).  The  quantity  WG(f$)  is  the  point-,  or  direction¬ 
less  intensity  spectrum  of  the  G  (=g+c)  surface  elevation  =  c  + 

One  choice  for  WG  is  the  Pearson-Moskowitz  spectrum  [8].  (For  a  wave 
number  spectrum  of  the  soli  ton  ripples,  see  (3.6b)  and  Section  3. IB  ff.) 

2.2  Scattering  Cross-Sections:  Definitions 


There  is  a  number  of  variations  on  the  concept  and  definitions  of 
scattering  cross  section..  It  is  therefore  important  to  define  the 
term  explicitly,  so  that  the  different  definitions  can  be  calibrated 
with  one  another,  as  we  shall  need  to  do  in  order  to  effect  comparisons 
with  both  the  analytic  results  of  others  and  the  corresponding  measure¬ 
ments  (cf.  Secs.  2.A,B  and  Sec.  3  ff.). 

We  begin  with: 
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A. 


Incoherent  Scatterinq  Cross  Sections, 
- •* -  i  nc 

Here  we  define 


;!°> 

me 


riirSh  (scattering  at  receiver  R)  •  rL  ?  -1 

inCOh  - —[(path  loss)  x  (gTgR)2]  , 


*inci 


dent(at  surface)  AREF 


(2.6) 


where  Iincoh»  incident  are*  resPect'* vely »  the  intensities  of  the  scat¬ 
tered  and  incident  fields  at  the  points  indicated.  The  basic  concept 
of  the  scattering  cross  section  (for  surfaces)  is  to  eliminate  the 
effects  of  source  level  and  propagation,  i.e.,  the  effects  of  the 
medium— which  are  handled  separately— when  computing  energy  loss,  and 
to  focus  on  the  effects  of  the  random  scattering  surface  itself.  For 
this  reason  path  loss  (absorption),  beam  pattern  gains  (gT,  gR),  signal 
levels,  and  source  and  receiver  distances  are  removed,  where  possible, 
as  (2.6)  indicates.  To  keep  ojj^  dimensionless,  a  reference  "illumina¬ 
tion"  area,  AREp,  is  employed,  whose  specific  form  is  suggested  by  the 
composite  beam  pattern  projection  on  the  reference  or  equilibrium  sur- 
face<c>=  0:  SQ.  Figure  2.2  shows  the  relevant  geometry. 


A 

Tz' 

1  !  z ' 

////  ! 

The  various  factors  in  Eq.  (2.6)  are  given  by 


(2.7a)  1 1 ncoh  (scat*  at  receiver)  =  Mx?<x>(0)  =  'intensity  of  the  (1n- 

coherent)  scattered  field  at  the 
receiver  (R),  cf.  Sec.  7.2  et  seq.; 
Sec.  8  ff . ; 


(2‘7b)  Iincident(at  Surface)  = 


^^in^^oT'  =  intensity  ^e  in" 
cident  field  at  0^,  on  the  equilibrium 


scattering  surface  SQ; 


(2.7c) 


gT,gR  =  the  aperture  "gain"  of  the  trans¬ 
mitting  and  receiving  systems,  cf. 
(6.3),  (6.7); 


(2.7d)  AREp  =  a  reference  area  on  the  equi¬ 

librium  surface  SQ,  cf.  Fig.  2.2, 
projected  by  the  composite  T  and  R 
beam  patterns.  (SeeC_ff.); 

•2a«?cT 

(2.7e)  ("path  loss"  =  e  0  0  cQT0  =  RqT  +  RqR;  cf.  Sec.  5.3,  Eq. 

*  (absorption)  (5.18);  u>0  (*2nfQ  )  is  the  (angular) 

frequency  of  the  emitted  signal;  a  is 
an  absorption  coefficient;  cQ  =  (mean) 
wave  front  speed  of  sound  in  the 
water  medium; 


(2.7f) 


RqR  =  distance  of  the  receiver  (origin)  from  0^, 
Fig  2.2; 


(2.7g) 


Rqj  =  distance  of  the  transmitter  (origin)  from  0S, 
Fig.  2.2; 


(2.7h) 


?(r,t)  =  (vector)  wave  surface  elevation,  cf. 

(2.1). 


The  definition  (2.6)  for  the  incoherent  scattering  cross  section, 

and  those  similar  to  it  ([3]-[5],  [8],  [9],  [23]-[25],  [29],  for  example) 

are  formally  independent  of  range  (RQy,  RqR).  This  is  not  an  inherent 

property  of  the  definition,  however;  it  is  a  direct  consequence  of  the 

far-field  assumption,  whereby  the  effective  coherent  scattering  area  is 

sufficiently  small  vis-a-vis  source  and  receiver  distances  (R0-[->  RoR) 

and  dimensions,  cf.  Sec.  5.5  and  Eq.  (5.31).  It  also  depends  on  the 

correlation  distances  («•  A,)  of  the  (components  of  the)  random  wave 
"COl  x  y 

surface.  Thus,  is  implicitly  a  function  of  geometry,  where  care 
must  be  taken  in  its  use  to  ensure  that  the  conditions  governing  the 
derivation  of  1  ^ ncoh  ^=MX-<X>^^ *  (2.7a),  are  obeyed* 

A  complete  definition  of  the  incoherent  scattering  cross  section 
includes  the  effects  of  diffraction 

o.  =  o(0)  +  £  °.^)  ,  (2  g 

incoh  incoh  incoh  v 

where  the  ,  k  3*  1,  are  formally  given  by  (2.6),  with 

l neon  (k>l^  incon 

replaced  by  the  scattering  intensities  I>n^  =  Mx_<x>(°)»  cf.  Secs. 

7.2C,  8.5  ff. 

Preliminary  estimates  (Sec.  3.3  ff.)  of  the  magnitude  of  the 

leading  diffraction  component,  ncoh  ( >>  °incoh^’  indicate  that  it 

is  ordinarily  considerably  smaller  than  the  components  of  the  "classical" 
"COl 

term,  ol”'oh,  so  that  in  this  study  we  shall  be  able  to  neglect  the 
diffraction  contributions. 


*The  factors  4ir  in  (2.7a),  (2.7b),  and  in  (2.9),  (2.10b)  ff.,  etc., 

arise  because  of  our  definition  of  the  green's  function,  (5.8),  and 

source  function  (5.3a)  in  the  equations  of  propagation  (here  a  Helmholtz 

medium).  Thus,  our  acoustic  field  is  =  a/4Tr,  where  a  is  derived  from 

a  green's  function  source  of  the  form  -4tt6 (R-R ’ ) 6 ( t- 1 ' ) .  However, 

~  (01 

because  of  the  particular  form  of  the  definition  of  a  used  here  and 
generally,  the  scaling  of  the  field  is  immaterial,  as  is  required  in 
any  useful  definition. 
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J5.  Coherent  Scattering  Cross  Sections,  0£q^: 

'■(n) 

The  coherent  scattering  cross  section,  o'  /,  is  formally  the  same 

-(Q)  COn 

as  01-p^o^»  (2.6),  except  that  now  the  incident  intensity  (2.7b)  becomes 


‘incident131  the  receiver)^  ‘  Ko<°)in/U,(RoT  +  "<*>>' 


(2-9) 


which  is  the  "mirror  reflection"  term.  Thus,  I 


mcid-coh 

(1/4>‘incid-incoh-  cf-  (^b>-  In  3<l<htl0"-  liScofi  is  >-eplaced  by 
the  coherent  component  M^^(0),  cf.  Sec.  7.2B,  Eq.  (8.7).  Since  only 

the  "classical"  or  (k=0)  component  of  the  scattered  field, contains  a 
potentially  coherent  contribution  (excluding  any  direct  field  which  may 
be  received  under  certain  mutual  geometries),  the  complete  coherent 
scatter  cross  section  is  now  specifically 


"(0)  =  IcSh(scattering  at  R)  RoR 


coh  I^„^dj(at  the  receiver 
incid' 


R]~AR£f  fpath  1oss  X  (9j9r)  }"  »  (2.10a) 


-<0)  .  RoR14'<RoR  *  rot>)2  "<>> 

COh  (9t9r)2  Aref  Kq ( 0 ) • n  •  (path  loss) 


(2.10b) 


The  Reference  Area,  ARER  : 

The  reference  area,  ARER,  appearing  in  the  above  definitions  of  the 
scattering  cross  sections,  (2.6),  (2.10),  while  arbitrary,  is  dependent 
on  the  beam  pattern  projections  on  SQ.  From  Section  6.5  IV,  (6.56),  our 


choice  of  reference  area  is  specifically 
aref  =  y2,  \  =  2*//AB(eoT)  , 

where  is  the  projected  area  (on  SQ)  of  the  combined  gauss ian-omni- 
directional  beam  pattern  used  specifically  in  this  study. 


(2.11) 


High-Frequency"  Scattering  Cross  Sections 


By  high-frequency"  we  mean  here  large  Rayleigh  numbers,  e.g., 

2 

Rg  =  (k0b0°g)  >:>  1»  specifically  for  the  large-scale  gravity  component 
of  the  wave  surface.  Thus,  we  have  for  the  elements  of  Rg: 

kQ  =  2tt/ =  2irf0/c0  =  wave  number  ( rad- meter' *)  of  the  (central) 

frequency  of  the  applied  (narrow-band) 
signal ; 


bQ  =  cos  eQj  +  cos  9qR‘,  cf.  Fig.  (5.1)  and  Sec.  5.4  ff.; 

2  2 

°G  =  ’  t^ie  mean-scluare  (gravity-capillary)  wave  surface 

height,  about  <£g>  =  0. 


Accordingly,  from  the  results  of  Sections  8.1,  8.3,  8.4  for  the 
scatter  intensities  Mx_<x>(°) ’  M<x>(°)  aPPlied  directly  to  the  scatter 
cross  sections  as  defined  by  (2.6),  (2.10b),  with  the  reference  area 
given  here  by  (2.11),  we  get  the  following  "high-frequency"  forms: 


a(0)  -  £  S2 

CTincoh  Rg»l  o  o 


(A2  +  a2  +  a2  \ 

If  ox  oy  oz  \ 
l  *oz'2  ) 


i  r<2°ox>2 ,  ( y2' 
,.2  .2  . 


8*bo”Gx°Gy 


K 

+  — Sj  Ni8'  (2“„  ),  ,  W,(2u  k  1 0) 

16*2  Gc  0  x,y  ,nc  c  0  0  x,y  G 


+  16!2  (2S°|x,y>inc  °s"s(2Vo|,J0) 


where  we  have  used  (8.6),  (8.7),  extended  to  (8.29),  with  ol  |  - 

(a  a  ),  a  two-dimensional  vector  in  the  reference  (x,y-)  or  S  -plane, 
ox y  oy  0 

The  terms  in  {  }q  represents  the  results  of  the  "physical-optics"  cum 
"perturbational "  approach,  cf.  Sec.  7.3  for  the  single  gravity-capillary 
wave  surface,  while  the  remaining  term  embodies  the  contribution  of  the 
soliton-ripples.  [See  (2.17)-(2.20)  below  for  explicit  descriptions 
of  the  various  elements  in  (2.13).] 

Here,  and  are  respectively  the  "tilt  factors"  for  the 
capillary  and  soli  ton-ripple  components  of  the  wave  surface. 

Specifically,  from  (7.33)  we  have  here  (<*>oT  =  tt/2) 


,(0)(o  )  3/-2!oxSl..+,  2a_0 

*gS  ^'incoh  \  1  +  ^2  + 


-<y  '  2aoz^ 


}  ♦  c2 
x  y 


S  16{aoy<4y>  +  6°oy°oz<cGy>  +Otoz  >  * 

2  2  2  -1  -2 

since  nz  =  (1  +  ^Gx  +  CGy)  =  1  -  0(10  )  =  1  for  the  accuracy  needed 

here.  Similarly,  we  get 

Ngc ’ > i ncoh* 2oox Sc  *  2<Vy  '  2,oz)2>G(2“oz)2 

'■  >64<vi  ♦  4>-  i  -<i>  • 


cf.  (8.28a).  For  the  gaussian  gravity  wave  elevations  assumed  here, 

4  4 

we  have  <t  >  =  3o  .  We  note  that  the  tilt-factor  for  the  capillary 
component,  (2.14b),  vanishes  as  eoT  -*■  w/2  (zero  grazing  angle),  because 
is  really  part  of  a  single  surface:  the  capillary  waves  do  not  "ride" 
on  the  gravity  components,  cf.  (2.1)  above.  On  the  other  hand,  the 
"tilt-factor"  N^,  (2.14a),  for  the  soliton-ripples  is  nonvanishing  (as 
0Qj  -*•  tt/2 ) ,  as  expected,  since  these  ripples,  or  hydraulic  jumps,  do  ride 
upon  the  gravity-capillary  surface,  tg.  These  "capillary"  and  ripple 
components  show  up  as  Bragg  (or  resonant)  scatter  terms  (-  wave  number 


spectra  at  k  =  2k  a  |  )  (Cl],  [2],  [2l]-[23],  etc.).  These  are  first- 

order  scatter  components,  since  o£,  <<  1,  with  small-scale  spatial 

correlations  [24].  The  specular-point  term  in  (2.13)  is  independent 
of  frequency,  generally,  in  this  "physical-optics"  limit. 

Also,  we  have  for  all  >  0  in  the  coherent  cases  (8.7),  (8.21), 


•(0) 

16„2 


-b2k2  2 

,2„  oVG 

2aoz  e 


■3ko[«ox/A  + 


O^'hi-freq", 


(2.15) 


2  2  2 

where  a  »  a  ,  o~,  since  R  »  1  here:  very  rough  surfaces  destroy 
9  co  g 

coherence,  as  expected.  [In  addition,  since  A,B  <<  1,  cf.  (6.32), 

the  second  exponential  also  ensures  the  vanishing  of  ocQh,  unless  the 

Snell  angle  (a  =  a  =  0)  is  chosen.]  Note  that  this  coherent-scatter 
ux  oy 

cross  section  depends  on  the  area  illuminated  (via  4,),  as  distinct 
from  the  incoherent  cases,  (2.13),  which  are  always  area-independent, 
subject,  of  course,  to  the  conditions  (2.20)  ff. 

In  compact  fashion,  we  can  rewrite  (2.13)  as 


o 

incoh  R( 


-  r;(°)  +  ;(°)  =;(0)  +;(0) 
^>>1  *■  g-inc  0c-inc^  S-inc  G-inc  S-inc 


(2.16) 


where  the  first  two  terms  of  (2.13)  correspond  t0  °g_jnc>  and  the  last, 
to  that  portion  of  the  scattering  cross  section  attributable  to  the 
soli  ton-ripple  component.  Our  results  (2 . 13)-(2. 16)  hold  for  arbitrary 
directions  of  illumination  and  observation. 

The  various  elements  of  (2.10),  (2.13)  are  specifically: 


i-  so 


=  mean-square,  mean,  reflection  coefficients  (=1,  for 
water/air  interfaces); 


(2.17)/  S  ,  S  =  mean-square,  mean,  shadowing  function,  see  Sec.  7.4C. 


tGx’  °Gy  =  <5Gj^*  <4/^  ^Gx  =  f?-’  etc-):  mean-square  slopes  of 
the  gravity-capillary  wave  component,  cf.  Sec.  7.4B; 

*2  =  mean-square  height  (about  =  0)  of  the 

"capillary"  wave  surface; 

i  =<4:  mean-square  height  of  the  soliton  humps; 

4  =  surface  wave  number  intensity  spectrum  of  c  (cf.  (2.3))  with 

v  L 

oo 

o?p  (ar,0)  =  f  W  ( k 1 0 )  cos(k*Ar)  ,  cf.  (8.25); 

G  c  -  V  c  ~  - 27t^ 


W<.  =  wave  number  spectrum  of  the  soli  ton-ripple 

oo 

»  ofw$  =  /  Ks(Ar|  0)e^*^r  d(Ar) 

_  00 

-  2 

and  w$  =  WSA*S  =  normalized  spectrum; 

Al  =  Eq.  (2.11); 

vbo  =  cos  eoT  +  cos  0oR,  (2.12). 

The  directional  vector  2a^  is  (cf.  Sec.  5.5,  Fig.  5.1), 


A  A  A  .  fj  1 

2«o  =  ioT  -  — oR  ’  -lx'1  *  IT -)  cos  *oT  S1"  eoT 


+  iy  f(1  +  R^>  s1n  *oT  sin  e0T  -Lo/RoR>  +lZbo 


(2.18) 


for  arbitrary  angles  of  illumination  and  observation. 


Important  special  cases  of  (2.18)  are: 


0).  Backscatter  (Rfr):  V  eoR  '  eoT;  JoR  '  ‘inT1  RoR  '  RoT; 


hR  =  hT’  ^oR  =  Vr  +  "/2;  cf-  Fig'  2'2)- 


•••  22„  -  2(ix  cos  *oT  sin  eoT  *  j.y  sin  *oT  sin  s>oT  *  lz  cos  eoT); 

(2.18a) 

(a2  +  a2  +  a2  )/a  /2  =  2/COS  e„T.  (2.18b) 

v  ox  oy  ozv  oz'  '  oT’  i0  ; 


1, 


(2.1! 


(ii).  Bistatic  at  the  Snell  Angle  (R  t  T) :  LQ^  0 

<ioR>x  =  <ioT>x  whe"  *oT  =  "/2  *  <*oR  -  ,/2); 
<IoR>y  =  <i„T>y  when  Lo=  <RoT  +  RoR>  si"  eoT; 
*-oR*z  "  -UoT>2  "hen  soR  =  9oT' 


Snell 

plane 


•■-22o  "  24  cos  eoT*  or  20  ■  4  c°s  9oT; 

( +  a*"  +  )/a  /2  =  2  COS  9  T  . 

'  ox  oy  ozv  oz  oT 


(2.1! 

(2.1' 


Otherwise,  (2.18)  is  the  general  relation. 

The  principal  assumptions  and  approximations  pertaining  to  our 
general  high-frequency  results (2. 13) ,  (2.15),  are: 


Far-field  (Fraunhofer)  geometries,  cf.  Sec.  5.5; 

Narrow-band  signals  (so  that  we  may  treat  time  parametrically 
in  the  moving  wave  surface  vis-a-vis  the  acoustic  signal); 
cf.  remarks  after  Eq.  (5.21b); 

Narrow  beams  (cf.  Sec.  6.6);  at  least  one  narrow-beam; 

Neglects  diffraction  terms:  k=l:  "Diffuse"  scatter;  2: 
multiple  scatter  (cf.  Sec.  3.3;  also  Sec.  8.5); 

Small  Rayleigh  numbers  for  the  small-scale  surfaces, 

^c,S  E  ^obo°c,S^  <<: 

The  capillary  and  soliton  surfaces  Uc»?s)  are  "small-scale," 

e.g.,  *  .  <<  Z  :  the  correlation  distance  of  the  gravity- 
c  fO  9 

wave  component  is  much  larger  than  that  of  the  capillary 
and  soliton  waves; 

The  small-scale  surfaces  are  statistically  independent  of 
the  gravity  wave  surface  component  (vide  remarks  in  Sec.  7.3A, 
Eq.  (7.19)  et  seq.); 

Both  components  of  the  wave  surface  are  essentially  homogeneous 
and  stationary,  at  least  over  the  "illuminated"  area  and  for 
times  long  compared  to  the  duration  of  the  incident  signal. 


[For  a  more  detailed  discussion,  see  Secs.  5.5,  6.6,  and  Sec.  7.J 


A.  Backscatter  (RgT) 

First,  by  (2.15)  the  coherent  backscatter  cross  section  3^ 
vanishes  for  these  high  frequencies.  On  the  other  hand,  the  incoherent 
backscatter  cross  section,  is  clearly  nonvanishing.  It  is 

obtained  at  once  from  (2.18)  in  (2.13).  On  setting  4>oT  =  n/2,  cf. 

Fig.  2.2a,  we  have 


(2.21) 


(2.21a) 


high-frequency  approximation,  which  is  recognized  as  the  "specular- 
point,"  facet-,  or  geometrical  acoustics  solution  form  of  the  full-wave 
approach  of  Bahar  [23],  [24],  [30],  [3l]  and  Barrick  et  al.  [5],  [22],  [24], 
[25].  The  second  term,  3c-lnc’  embodies  the  (first-order)  Bragg  scatter 
associated  with  the  small-scale  small  Rayleigh  number  capillary  wave  sur¬ 
face,  while  the  third  term,  ,  gives  the  Bragg  scatter  due  to  the 

5-inc  /qx 

soliton  ripples.  Whereas  the  specular-ooint  contribution  9'  '  is 
r  g-inc 

the  dominant  part  of  the  cross  section  at  moderate  angles  (0qT  <  60°), 

it  rapidly  vanishes  for  small  grazing  angles  (e_T  ■+■  ™/2),  cf.  Fig.  3a 
JO)  ""  0T 


of  [5]  ,  leaving  d^jnc  and  ®s?]nc,  of  which  only  the  latter  remains 
nonvanishing*as  eQj  -*•  u/2,  cf.  (2.14).  In  fact,  as  we  shall  presently 
see  (Sec.  3),  for  (0qT  ~  70°),  where  dg? ^ n c  4  °*  ds-i~-  dom1nates 


inc 


by  o(10  or  more  db).  It  is  this  latter  component,  ,  with  which 

♦Apart  from  the  shadowing  function  which  -*>  0,  of  course,  as  0oT-/2 


(0) 

c-inc 


we  are  chiefly  concerned  here.  Some  preliminary  numerical  evaluations 
and  comparisons  are  given  in  Secs.  3.1,  3.2,  following.  In  particular, 
from  (7.33)  we  easily  find  for  the  isotropic  cases  that  =  tf/2) 


**Gc-back  5  16  «s2  9oT  <°g(y.x>  s1"2  “oT  *  V 


oT  v  g(y=x) 


(2.22a) 


NGS-back  £  16  13  ag(y=x)  S1"4  eoT  +  6  ag(x=y)  cos  eoT  Sln  eoT  +  cos4  eoT} 


(2.22b) 


(For  the  more  general,  nonisotropic  cases  we  must  use  (2.14a,b).) 


B.  Bistatic  Scatter  in  the  Snell  Direction  (R^T) 

Here  we  use  (2.19)  in  (2. 13) [with  the  help  of  (8.14)  alternatively 
in  (2.6)],  since  artv  =  anxi  s  0,  aQz  =  cos  eoT,  to  get  directly 


e(°) 

incoh  R^T: Snell 
Rg»l 


4  /H- - ! - ? - +  T“S4  eoT  “c(0’0|0,!i 

L,8™GxoGy  Cos4  e0T 


k 

-§  cos4  eQT  o^ws(o,0|0)J  . 


At  small  grazing  angles  (eoT  -  */2)  the  facet  term  vanishes  because  of 
shadowing  and  the  Bragg  scatter  terms  likewise  go  to  zero,  not  only 
because  of  shadowing  but  also  because  their  "tilt-factors"  vanish. 

There  is  no  coherent  component,  of  course,  since  Rg>:>l,  cf.  (2.15). 

However,  off  the  Snell  angle,  although  the  facet  term  vanishes 
radpily  (as  90y  -*■  */2),  cf.  (2.15),  the  Bragg  terms  remain  (eQj  <  85°), 
where  shadowing  becomes  strong,  and  +  0  as  0Qy  -*•  ir/2,  since 
(a  ,onJ  t  0.  The  quantity  w(0,0  0)  is  proportional  to  t?,  cf.  (3.6b) 

Ua  UY  »  v» 


o 


o 
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Low-Frequency  Acoustic  Scatter  Cross  Sections 


Here  we  apply  the  results  of  Sections  8.2,  8.4,  especially  Eq. 
(8.9),  where  =  1  and  S2  =  -  1  (for  0qT  <  85°),  so  that 

l  0»  (8.9a),  for  these  cases  of  small  Rayleigh  numbers  R  c  << 

K5  9 )C| J 

We  obtain  the  following  incoherent  and  coherent  cross  sections,  from 

(2.6),  (2.10): 


R_  „  C«1 
9  *c  »S 


T  icbo°si"g(2^ko|x,yl0) 


+  n!v^(2  a  ,  .(2ak  |  JO)].  +  o|wJ2«  k 

Gc  '  -»ol  x.y'inc  c  ,*o  o|x,y'  /JG  GS-mc  S  S  —o 


(2.24) 


2  2  f  01 

where  ag  =  <<;g>,  and  (2.14)  gives  the  respective  tilt-factors,  Nv  ', 

cf.  (7.32)-(7.33)  for  details.  When  gravity  waves  are  present,  the 

contributions  of  the  "capillary"  component  and  the  soli  ton-ripples  can 

usually  be  neglected,  so  that  (2.24)  reduces  to 

*(0)  i  Vo  2 -  lnv  , 

dinc|Rg«l  "  lgir2  °gwg  (2ko^o  |x,y  ( 


(2.24a) 


Other  combinations  of  the  various  terms  of  (2.24)  are  possible,  accord¬ 
ing  to  the  presence  or  absence  of  the  soli  ton  ripples  or  the  gravity- 
capillary  wave  surface. 

Here  all  components  of  the  composite  surface  appear  as  Bragg  (or 
resonance)  scatter  terms,  as  expected.  However,  at  extreme  grazing 
angles  s  ->-0,  so  that  (2.24),  (2.24a)  vanish  because  of  shadowing. 

At  somewhat  larger  grazing  angles  ( 4>  =  10°,  say)  the  gravity-wave  com¬ 
ponent  can  be  ignorable  vis-a-vis  the  soli  ton-ripple  contribution, 
because  of  the  factor  b£  (=(cos  e  y  +  cos  0or)4)»  cf.  Eq.  (3.2c)  ff., 
which  in  turn  is  considerably  larger  than  the  "capillary"  contribution. 
The  latter  can  be  important,  however,  if  the  local  wind  conditions  are 
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Here  again  we  are  sufficiently  in  the  far-field  to  set  JiTR  =  0,  cf.  (2 

From  (6.11a),  (6.12a),  (6.17),  we  note  that  A- a  2  ~  (cos  e  T+cos  e  _)2 
,  *  ol  o R 

generally.  Since  A,B«1,  cf.  Sec.  6.3,  Eq.  (6.17),  the  exponential  t 

in  (2.25)  ensures  that  3^  is  very  small,  unless  (aQX  *  aQy  *  0). 
Accordingly,  in  special  cases  we  have: 

A.  Backscatter 


/ 

incoh  R@T 


"low- freq" 


=  s  Ro  77T  i[16  °g  cos  eoT  V0’  2ko  sin  eoT  0) 


*  s1"9ot|°)1g 

+  NGS)(eoT),„c°S*'s<0-2kosi"eoTl0>l’ 


C»i.iT»Tn,rhrKn«n 


where  now  the  tilt-factors  are  given  by  (2.14a,b)  and  ( 7 . 32 ) - (7 .33) . 

For  the  coherent  component  of  backscatter,  we  see  that,  similarly, 
.251  reduces  to 


2  .  2 

/ ,  0  o  n  _3k  sin  6  T/B(0  T) 

0^1  -  -  R  ^  s  A  ( ft  1  A  a  0 

enh  CIST  0  b  AT'9r.T^  C0S  ®nTe  i 


coh  R@T  v2  no  °  lv  or 

"low-freq" 


Rg  «  1, 


which  reduces  still  further  with  the  help  of  (6.11)  in  (6.18),  to 


AO)  i  i  „  2  .  P  P  2  2  -k 

c0h  0  (aT  cos^  0oT  +  bT  sin  e0T^  2  cos  9( 


-3k2  tan2  eoy/2A<a2  cos2  QqJ  +  b2  sin2  0qT). 


(2.27a 


This  in  turn  reveals  the  explicit  dependence  on  0  ,.  In  particular, 

2  IQ)  0 

A  =  2AT/RQy  «l,  so  that  here  is  always  small,  unless  0Qy  =  0 

(vertical  incidence),  cf.  Fig.  2.2a.  Then  (2.27),  (2.27a)  become* 


coh  vert 
9oT=0 


tW  ■  °>  ‘rth  ■  R,  “  '• 

IT  T  3 


(2.28 


since  S  =  i  and  ^  =  1:  there  is  no  shadowing  at  vertical  incidence 
and  |i?  |  =  1  for  these  water/air  interfaces.  Of  course,  for  large 

0  /  Q\ 

Rayleigh  numbers  (Rg  >>1),  -*-0  according  to  (2.15)  as  expected: 

when  the  surface  is  rough,  coherence  is  destroyed. 


*For  aT,  by,  see  Eqs.  (6.7),  (6.11a),  (6.12a);  0  <  ay  <  1;  0  <  by  <  1 
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Bistatic  Scatter  in  the  Snell  Direction 
From  (2.19)  we  have  here  aQR  =  eQT;  eoT  =  */2  (=  ^0r-7I/2)»  with 

aox  =  aoy  =  “oz  =  cos  6 oT*  Lo>  Thus*  w^th  the  he^P  (6*19)  we 
see  that  (2.24)  for  the  incoherent  scatter  cross  section  reduces  directly 

to 


with  the  help  of  (2.14a,b).  This  is  the  "low-frequency"  (R  <<1)  ana¬ 
logue  of  (2.24)  above,  showing  again  the  (first-order)  Bragg  scatter 
contributions  of  both  wave  surface  components,  e.g.,  the  gravity¬ 
capillary  surface  (6)  and  the  sol i ton- ri ppl es (S), which  ride  upon  it. 

Here  from  (2.14),  =  16  a4y  =  16  cos4  e^.  The  normalized 

spectra,  w,  are  seen  to  be  proportional  to  the  (mean-square)  correlation 
scales  of  the  respective  surfaces.  Because  of  the  small  Rayleigh  numbers 
for  the  large-scale  component  (g),  the  tilt  effects  are  negligible  here. 
As  grazing  incidence  is  approached  (e^-*- tr/2) ,  djp^oh  ^inishes  both 
because  of  the  cos4  factor  and  from  the  fact  that  S*  -*0,  ultimately. 

For  the  coherent  scatter  cross  section  in  the  Snell  direction, 

(2.25)  now  reduces  to 


(i?0  =  1),  Rg<<  1,  where  we  have  used  (6.19)  again.  Equation  (2.30) 
also  vanishes  as  eoT  -*ir/2,  as  expected.  In  the  special  case  where 
RoR  =  RoT*  simplifies  further  to 


=  S  OT  a  1.2 2  a  ,  u  2  .j.2  a  \  “H 

0  2.  cos  6oT  ^aT  C0S  9oT  bT  S1n  9oT^ 


coh  Snell  — 5 —  cos  0  -  (aT 

"low-freq"  2*  AT  01  1 


Rg  «  1- 


the  special  case  where  6oT  =  0  we  get  the  backscatter  result  (2.28) 


2.5  Analytic  Comparisons  with  Recent  Results 


Analytic  comparisons  with  earlier  work  now  provide  us  with  various 
checks  on  the  accuracy  and  assisnptions  of  the  present  analysis.  In 
addition,  they  indicate  a  variety  of  differences,  which  stem  mainly 
from  our  different  choice  of  physical  model,  as  well  as  a  number  of 
technical  differences  in  the  analysis  and  some  definitions.  These 
will  become  apparent  as  we  continue: 

A-  Tolstoy  and  Clay  [9] 

We  consider  first  incoherent  backscatter.  From  Eq.  (6.75)  of  [9], 
wherein  «'2>  =  vback  -  k„  cos  eoT,  «2/2v2  -  (a2x  +  )/2  cos2  e 

in  our  notation,  with  $0y  =  tt/2,  here,  so  that  aQX  =  0,  aQ  =  sin  e^, 
we  see  at  once  that  9 

_  tan  6oT 

...  s2  e  ^(r>) 

(back):  dT+C,inc|R  »1  S  Eq’  (6,75)|t+C  =  ~  4 

1  9  1  87TOGyaGx  cos  eoT 


tan2  e. 


(2.32 


Bo(sec4  0oT  e 


c2.,2 

’s  '2oG(x=y)’ 


which  is  just  the  first  term  of  our  more  general  result  for  truly  two- 
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component  surfaces,  (2.21)  above,  for  the  now  isotropic  (gravity)  wave 
surface  at  "high-frequencies." 

Moreover,  since  our 

ct2  +  cT  +  a2 

— - ^oy/-2--aoz  =  2f(e)T+c,  Eqs.  (6.23),  (6.25),  [9],  (2.33) 


we  readily  see  that  the  Tolstoy  and  Clay  result  (6.74)  [9]  for  general 
bi static  (incoherent)  scattering  in  the  "high-frequency"  regime 
becomes 


c2  -  on2 

5  "  2  g(x*y)  ’ 

which  is,  again,  precisely  the  first  term  of  our  result  (2.21).  ^ n  the 
earlier  work  the  effects  of  shadowing  are  neglected,  e.g..  S’2  =  s  -  1. 

Similarly,  for  the  "low-frequency"  cases  ( Rg  <<  1 )  we  have  in  our 
terminology,  rQ|j+Q  4  R//2,  cf.  (6.56),  [9],  so  that  (with/?^  rein¬ 
serted  in  (2.24))  we  get 


(bistatic): 


d 


(°)  . 

T+C,inc |Rg<<l 


Eq.  (6.58),  [9] 


rW 
0  0  0 


Wc(Jc|0), 


(2.35) 


This  is  once  more  just  the  first  term  of  our  more  general  result 
(2.24).  (Here  we  have  replaced  2f(0)y+£  by  bQ  =  cos  0^  =  cos  0qR, 
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with  y  =  bQ/2 ,  since  in  these  low-frequency  cases  we  do  not  integrate 
by  parts  (e.g.,  (7.27d),  etc.  vs.  p.  198,  [9]),  employing  essentially 
the  Eckart  approximations  n  •  2aQ/nz  =  -2aoz,  [32].) 

The  coherent  scatter  cross  sections  (with  fy+£  =  bQ/2  again,  Ay+£  = 
4j)  at  all  frequencies,  e.g.,  Eq.  (6.61)  of  Tolstoy  and  Clay  [9]  vs. 
(2.15)  above  likewise  agree,  except  in  the  directional  exponential:  our 
result  is 


3kV  ?  ? 

— 5 —  (a2  +  a2  ) 

2  '  OX  oy1 


VS. 


-2k^R2 (ot2  +  a2  ) 
o  '  ox  oy' 


T+C 


(2.36) 


The  discrepancy  arises  because  no  approximations  are  made  in  the 
evaluation  of  our  beam  pattern  integrals,  cf.  1^(20^),  (6.38),  (6.39), 
whereas  approximations  are  made  in  the  former  (T+C)  evaluation. 

In  general,  Tolstoy  and  Clay  [9]  also  employ  the  familiar  Tangent 
Plane  (Kirchoff)  approximation  in  arriving  at  their  results,  as  do  we 
here.  Since  they  assume  a  single-scale  wave  surface,  they,  of  course, 
do  not  obtain  the  capillary  and  nonvanishing  (as  0Qj  -*•  n/2)  (ripple-) 
components  with  decreasing  grazing  angle  ( 4>  =  w/2  -  eoT)»  cf.  the  second 
and  third  terms  of  (2.13)  above. 

B.  McDaniel  and  Gorman  [3],  [4] 

Using  a  two-scale  surface  model,  McDaniel  and  Gorman  [3],  [4]  omit 

the  specular-point  contribution  (viz.,  the  first  term  of  (2.21))  and 

concentrate  on  the  small-scale  capillary  wave  component,  which  is 

dominant  at  high-frequencies  and  small  grazing  angles  in  their  model 

(but  not  in  ours).  Their  evaluation  of  the  backscatter  cross  section 

for  this  capillary  term  embodies  a  Kirchoff  approximation  for  the 

large-scale  gravity  wave  surface  component,  and  a  perturbation  tech- 

nique  [requiring  small  rms  heights  (oZ.  «  1)  and  slopes  (a.v  .  <<  1)1 

v*  cx  ,y 

for  the  small-scale  component.  Moreover,  they  replace  the  "tilt-factor" 
in  (2.21)  by  an  average  over  a  local  grazing  angle  (as  suggested 
by  Bachmann  [29]),  which  contains  an  estimate  of  the  shadowing  effects. 
Their  treatment  postulates  the  capillary  waves  to  be  a  separate  wave 
structure  from  the  gravity  wave  component,  on  which  the  former  con- 
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sequently  ride  [so  that  if  this  were  the  case  the  appropriate  "tilt- 
factor"  would  have  the  form  of  Ng(s=c)»  (2.14a),  now,  rather  than  , 
(2.14b)].  Our  physical  model  [cf.  Sec.  2 . l]  treats  the  so-called 
"capillary"  waves  as  part  of  the  continuous  gravity-capillary  wave 
surface  regime,  where  the  former  evolves  naturally  into  the  latter  via 
the  (nonlinear)  mechanism  of  the  soliton-ripples  generated  locally  by 
the  near-surface  wind. 

Noting  [cf.  Sec.  7.3  and  (2.2a,b)]  that  our  tilt-factor  is 
derived  from 


/nl>G=g+c  *  16  ^JGg  ‘-o-back^G  = 


in  the  McDaniel -Gorman  notation  (cf.  (25),  [4]  and  from  (16),  [4]),  and 
that  their  wave  number  surface  spectral  density  w(2JsA)McD+g  = 
(2ir)-1Wc(kj0)Ml.d,  (7.51),  we  find  at  once  that,  formally, 

k4 

[3^  term  of  (2.21)  i  S<?>  ]  s  <(«  •  k/k0)%  2,«s(21Sl)McD+G 

7T 

(2 

=  f  ^£g  ’  i>%  ws^2-x^mcD+g  . 

Here  our  actual  ripple-surface  (S)  replaces  the  "capillary"  component, 
which  in  their  model  rides  on  the  gravity-wave  structure.  [McDaniel  and 
Gorman,  of  course,  use  wc(2kA)^CQ+g  in  place  of  ws(2Kj.)mcd+G’  where 
a2  <<  a2,  cf.  Secs.  3.1,  3.2,  and  omit  the  soliton-ripples.] 

Thus,  our  small-scale  results  (in  (2.21))  reduce  essentially  (but 
not  exactly)  to  that  of  McDaniel  and  Gorman  [3],  [4]  if  we  were  to 
replace  our  ripple-surface  by  a  similarly  generated,  now  independent, 
capillary  wave  surface  (S  ■*  c),  and  drop  the  "capillary"  component  in 
(2.21),  so  that  only  the  specular-point  scatter  term  and  the  (S  c) 
small-scale  terms  remain.  [Since  high  frequencies  and  small  grazing 
angles  are  primarily  considered  in  their  work,  the  specular-point 
scatter  contributions  may  be  dropped  as  negligible  vis-a-vis  the  other 
components,  cf.  (3.1)  ff . ]  The  various  analytic  differences  between 
the  two  treatments  arise  from  the  following: 


(ii). 


(iii). 


(iv). 


Our  approach  uses  the  Tangent  Plane  (Kirchoff)  method  on 
the  total  surface,  including  the  soliton  ripples.  This 
leads  to  a  "tilt-term"  of  the  form  (2.37)  above,  which, 
however,  reduces  here  to  the  simple  result 


in  all  such  cases,  since  n~*  =0(1),  cf.  (2.2b),  (2.14a), 
etc. 

Reference  [4]  also  outlines  a  similar  approach,  employing 
the  Kirchoff  method  on  the  total  surface,  vide  Appendix 
A,  [4].  Moreover,  a  tilt-factor  of  the  form  N^,  (2.14b) 
is  correctly  obtained  (however,  with  a  factor  k2  omitted 
in  Eq.  (A14),  [4],  because  of  a  corresponding  omission  in 
going  from  Eq.  (A10)  to  (A12)  in  the  analysis  [4]).  This 
tilt  factor,  of  course,  is  incapable  of  ensuring  a  non¬ 
vanishing  result  (as  eQT  -*•  tt/2),  as  we  have  already 
remarked  above,  (2.14b)  et  seq.,  and  also  the  discussion 
in  Sec.  2.1.  There  is  now  no  separate  ripple-term,  how¬ 
ever,  in  the  analysis  [4],  Appendix  A.  What  was  the 
separate  "capillary"  wave  surface  has  now  been  absorbed 
properly  into  the  gravity-capillary  wave  continuum.  This 
approach  does  not  require  small  slopes,  or  small  heights, 
unlike  the  perturbation  techniques  [4],  [33], 

A  jA 

We  perform  the  "tilt"-averaging  directly  on  (nG*  20^)  , 
etc.  (as  does  Kur'yanov  [21],  [33]),  without  the  ad  hoc 
introduction  of  a  local  grazing  angle,  which,  more¬ 
over,  does  not  account  for  the  general  anisotropy  of  the 
wave  surface  slopes. 

In  the  present  treatment  the  shadowing  and  reflection 
coefficients  are  handled  by  the  (approximate)  methods  of 
Sec.  7.4Q.,  based  on  Bass  and  Fuks  [1]  as  indicated.  (Un¬ 
less  we  are  dealing  with  angles  larger  than  85°  (say  88°), 

— 7  —2 

the  shadowing  effects  are  all  small  and  S  -  S  -  1. 


The  really  critical  difference  between  our  approach  and  that  of 
McDaniel  and  Gorman  [3],  [4]  lies  in  the  choice  of  physical  model,  since 
there  are  no  major  differences  resulting  from  the  /arious  analytic 
methods  employed.  We  choose  a  soli  ton-ripple  mechanism,  in  addition 
to  the  gravity-capillary  wave  surface,  as  explained  in  Sec.  2.1.  They 
employ  only  the  latter,  treating  the  "capillary"  component  as  a  sepa¬ 
rate  component  riding  on  the  gravity-wave  structure,  rather  than  as  a 
continuous,  high-frequency  part  of  this  wave  surface  as  a  whole.  They 
explain  observed  discrepancies  between  experimental  and  theoretical 
cross  sections  in  terms  of  a  near-surface  bubble  mechanism.  This  is 
certainly  a  possible  mechanism,  but  not  necessarily  the  only  one,  par¬ 
ticularly  when  (sufficient)  bubbles  are  not  present  [6].  (See  Secs. 

1  and  4.1  ff. ) 


C_.  Kur'yanov  [21],  [33] 

Kur'yanov  employs  a  two-scale  wave  surface  model,  implicitly  in 
the  form  (2.13)  (where  we  omit  the  "capillary"  term).  He  uses  the  per¬ 
turbation  technique,  which  requires  small  (rms)  wave  heights  and  slopes. 
For  the  contribution  of  the  small-scale  surface,  whatever  the  physical 
mechanism  may  be,  he  obtains  a  tilt-factor  of  the  form  (in  our  notation) 


16  <(Hg  ’So-back^^G  "  16  <^G  '  “o-back^/nz>G-Mid 


‘  16  <(V20-back)4/n^G 


(2. 


cf.  Eq.  (33.22)  of  [l],  again  since  nQ  =  ixcGx  +  iySGy  -  j_2,  cf. 

(2.2b)  here.  This,  in  turn  with  the  fact  that  Kur'yanov's  and  our 
(normalized)  wave  number  spectra  are  related  by  4n2w(kJj0)^  =  w(k JO)^, 
reduces  our  result  (the  third  term  of  (2.21)  here)  directly  to  the 
form  (2.38)  above,  with  w(2XA)McD+G  replaced  by  2TTOjWs(k1|0)K,  or  to 
Eqs.  (33.16),  (33.22)  of  [l].  Thus,  Kur'yanov's  result  and  ours  (2.21) 
are  essentially  equivalent  analytically,  although  different  methods 
(e.g.,  Kirchoff  vs.  perturbation)  are  used  to  achieve  them.  Of  course, 
physically,  Kur'yanov  says  nothing  about  the  specific  soli  ton-ripple 


mechanism  itself,  nor  does  he  consider  the  gravity-capillary  wave  sur¬ 
face  directly,  and  split  it  into  two  separate  physical  components,  in 
the  calculation  of  the  small-scale  Bragg  scatter  component. 


JL  Barrick  and  Bahar  [5],  [22],  [23],  [25],  [30],  [31] 

A  little  later  than  Tolstoy  and  Clay  [9],  and  by  a  different  ap¬ 
proach,  Barrick  [5]  obtained  the  high-frequency  or  "optimal-limit" 
contribution,  namely  the  so-called  "specular-point"  scattering  com¬ 
ponent  of  the  total  scattering  cross  section  This  is  the 

first  term  of  (2.13).  Subsequently,  by  a  more  general  approach, 

Bahar  [23],  [30],  [31]  has  shown  that  for  the  conditions  assumed  here 
(cf.  (2.20)),  the  two-component  (gravity-capillary)  portion  of  our 
("three-component")  result  (2.13)  can  be  alternatively  derived,  and 
similarly  for  the  third,  independent  component  postulated  here.  In 
these  papers  Bahar  has  considered  the  more  general  problem  of 
electromagnetic  (EH)  scattering  from  random  rough  surfaces. 

If  we  associate  formally  Bahar's  (and  Barrick's)  results  for  hori¬ 
zontal  polarization  (HH)  with  our  results  for  acoustic  scattering,  we 
see  that  the  various  components  of  (2.13)  are  essentially  identical,  in 
form,  to  the  corresponding  EM  results,  provided  that  we  note  that  the 
definitions  of  the  scattering  cross  sections  in  the  EM  case  contain  a 
factor  4ti  (sterradians) ,  cf.  (32)  of  [30]  with  (2.6),  (2.10)  here. 
Moreover,  the  relationship  between  Bahar's  and  our  surface  wave  number 
spectra  is  n 2W ( k. j_| 0 ) B  =  W(k jj0)Mi(1,  cf.  (43)  of  [30],  with  (7.39), 
(7.41)  ff. 

As  a  specific  example,  let  us  consider  again_the  case  of  back- 
scatter,  cf.  (2.21).  Thus,  for  weak  shadowing  {Sr  =  l  now)  we  have 
from  (2.13),  omitting  the  additional  soli  ton-ripple  term. 
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for  isotropic  gravity-capillary  wave  surfaces  in  which  s  = 

etc.,  cf.  (2.32).  With  very  small  slopes  the  tilt  factor  in  (2.40) 

4 

reduces  to  cos  e0j,  so  that  the  second  term  becomes  precisely  the 
horizontal  polarization  component  of  Eqs.  (47),  (69)  of  [30],  while 
the  first  term  of  (2.40)  is  just  Barrick's  result  (136),  [5],  or  Eq. 
(72)  of  [30],  where  4>qT  =  V2,  =  (0,  2kQ  sin  0qT). 

With  the  alternative  two-scale  model  in  which  the  small-scale 
component  actually  rides  upon  the  large-scale  surface,  the  tilt  factor 
has  the  form  N^,  (2.14a),  which  is  non-vanishing  as  -*•  */2.  (The 
counterpart  of  this  in  Bahar's  results  appear  in  his  shadowing  terms 
<s(PQ-HH)>>  cf.  Section  4  of  [23],  and  Eqs.  (44c)  et  seq.  therein.) 
Thus,  in  general  our  approach  yields  results  equivalent  to  the  more 
general  (and  less  explicitly  tractable)  analysis  of  Bahar  [23],  [30], 
based  on  a  "full-wave"  theory  development,  where  now  Tangent  Plane 
(i.e.,  Kirchoff)  methods  can  be  employed  and  the  small-scale  wave  sur¬ 
face  components  have  small  rms  heights  (but  not  necessarily  small 
slopes),  as  is  the  case  of  our  model  here.  We  remark,  however,  that 
in  the  case  of  ocean  wave-surfaces,  Bahar  et  al .  [25]  use  a  single 
(i.e.,  gravity-capillary)  wave  surface,  which  is  then  artificially 
split  into  two  components  for  analytic  convenience  (cf.  [23],  [24], 
also).  The  small-scale  component  cannot  then  be  considered  as  riding 
upon  the  large-scale  component,  cf.  our  discussion  in  Sec.  2.1  above. 


Finally,  we  remark  that  except  for  our  and  Tolstoy  and  Clay's  [9] 
analyses,  and  that  of  Clay  and  Medwin  [l 1  ] ,  which  explicitly  introduce 
(gaussian)  beam  patterns,  e.g.,  specific  apertures  (cf.  Sec.  6  ff.), 
the  other  treatments  cited  here  assume  uniform  beam  patterns  over  some 
solid  angle,  vanishing  outside  this  angle.  [This  is  implicit  in  the 
choice  of  reference  area,  A„cc  [(2.7d)  and  Sec.  2.2C],  in  the  defini- 
tion  of  the  scattering  cross  section,  dv  ,  cf.  (2.6),  (2.10).]  Recon¬ 
ciliation  between  these  different  choices  of  beam  pattern  is  achieved 
through  (2.11),  i.e.,  by  appropriate  choice  of  AREp  in  the  definition 


wmmm 


Since  we  are  concerned  in  this  study  primarily  with  the  possibility 
of  obtaining  reasonable  agreement  between  theory  and  experiment  without 
the  intervention  of  scattering  by  bubbles,  as  explained  above  in  Section  1, 
we  shall  consider  only  the  high-frequency,  small  grazing  angle  regime, 
where  the  discrepancies  between  earlier  theory  and  experiment  are  most 
pronounced  [l]-[4].  Moreover,  it  is  sufficient  to  consider  one  typical 
situation  numerically  in  order  to  illustrate  such  discrepancies  and 
potential  agreements. 

Accordingly,  we  shall  treat  the  following  specific  example:  (inco¬ 
herent)  backscatter,  with  0Qj  =  80°  ( <f>  =  10°  grazing  angle)  with  <t>oT  = 

V2,  fQ  =  20  kHz,  central  signal  frequency;  average  near-surface  wind 
velocity  v  =  10  m/sec  (4  20  knots),  cf.  Fig.  2.1a.  Equation  (2.21)  is 
the  appropriate  relation  here  for  the  (back-)  scattering  cross  section 
3incoh"  Furthermore,  in  this  example  we  may  ignore  the  negligible  con¬ 
tribution  of  the  specular-point  (or  facet)  scattering  term,  to  write 
(2.21)  as 
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(3.1) 


where  the  tilt-factors  N^,  are  given  by  Eqs.  (2.14a,b);  (the 
unnormalized  two-dimensional  wave  number  spectra  are  wc  $  =  °c  s*c  S’ 
respectively).  As  noted  above  (Sec.  2.3)  the  (non-zero)  initial  term 
of  (3.1)  provides  the  first-order  Bragg-scatter  contribution  of  the 
high-frequency,  or  capillary  portion  of  the  underlying  single  wave  sur¬ 
face,  while  the  last  term  gives  the  (Bragg-)  scatter  associated  with  the 
essentially  independent,  omni-directional  soliton  ripples,  or  "hydraulic 
jumps"  [cf.  Sec.  2.1],  which  ride  upon  the  gravity-capillary  wave  surface. 

Our  task  here  is  to  evaluate  (3.1)  and  compare  the  numerical  results 


with  surface  models  which  do  not  include  the  soliton  ripple  mechanism 
(cf.  Sec.  2.5). 


3.1  Backscatter  with  Soliton  Ripples:  Numerical  Examples 

In  order  to  compare  with  currently  listed  experimental  data  [4]  we 
shall  renormalize  our  wave  number  spectra  via  the  relation  W(k>|0)Mid  = 
cf*  (7-51b)*  In  addition,  for  our  specific  numerical 
example  we  need  the  following: 

j4.  Calculations:  [fQ  =  20  kHz;  0qT  =  80°;  v  =  10  m/sec] 


(1).  k  =  p-  =  2*f  Jc  =  (2ir)  20  ♦  103/ (1.5)  103  =  83.78  rad  m"1; 

0  A  0  O 


k4  =  (4. 93) 107  (rad  m'1)4 


(3.2a) 


(2).  <t>Qj  =  V2:  6  =  80  :A«nv  =  0;  anv  =  sin  0nT  =  0.985,  Ccf.  Eq.  (2.21)3 


ox 


oy 


oT 


(3.2b) 


sin  0qT  =  0.985;  sin2  0Q-j-  =  0.970;  sin4  0Ql-  =  0.941 
cos  0qT  =  0.174;  cos2  0qT  =  3.02  •  10-2;  cos4  eoT  =  g.09  •  10-4;) 


(3.2c) 


(3).  From  Cox  and  Munk  [35],  esp  Eq.  (6.76)  of  [9],  we  have  for  the 
mean-square  slopes: 


°g(x=y)  =  (3*°  +  5-12v)10‘3:  = 


5.42  •  10 


-2 


(3. 2d) 


(4).  From  (2.14a)  and  (7.56a)  we  get  here  specifically 


NGS-inc  =  16  {3(5.42)210"4(0.941) 


+  6(5.42  *  10_2)(0.970)(3.02  *  10"2)  +  9.09  *  10-4  > 
=  16  {  8.29  ‘  10"3  +  9.53  *  10'3  +  0.91  *  10‘3) 


16(1.87  *  10‘2) 


(3.2e) 


In  the  same  way  we  have  from  (2.14b) 


NGc-inc  =  16  {E(5*42*  10'2){0.970)  +  3.02-  10-2]  3.02-  10'2) 
=  16(2.50  *  10"3)  . 


(3.2f) 


£.  Wave  Number  Spectra 

For  the  set  of  soli  ton  ripples,  or  hydraulic  jumps,  which  ride  upon 
the  rough  gravity-capillary  wave  surface,  we  develop  the  following  ele¬ 
mentary  second-moment  model . 

A  reasonable  approximation  of  the  wave  shape,  or  "hump,"  is  given 
by  the  gaussian  wave  form,  for  a  typical  wavefront  traveling  in  some 
direction  <t>,  viz: 


i(r.t) 


Vse 


■4r2/L2 


Al(r-c$t)  , 


(3.3) 


fj 


where  i  =  1  cos  +  j.  sin  4>,  r  =  «^x  +  y 

r  y 


,  and  where  Ai  is  a  unit 
distance  "window"  which  moves  along  with  the  soliton,  at  speed  c$.  It 
is  nonzero  for  |r  -  c<.t|  <  2L,  where  2L  is  the  spread  of  the  hump  to, 
say,  within  e"4  (=  2%)  of  its  maximum  value,  cf.  Fig.  16,  [16],  p.  465,  etc. 

We  determine  the  second-moment  functions  for  an  (ensemble  of)  single 
typical  humps  from 


M$(Ar,t)  =  Cdlj.tj)  •  i(r1+Ar,t1+At)>  (3.4a) 

5?  •  -4r?/L2  -4(r.+Ar)2/L2 

■5U</  e  “V*  t, 

-oo  1 


£  -4Ar2/L2  »  -8(r2+rnAr)/L2 

=  2T  pL(t)  e  /  e  1  1  dr:  (3.4b) 

•  00 


where 


pl(t)  = 


l/T,  -  1  - 


2L 


cs | t j  <  2L;  =  0,  otherwise.  (3.4c) 


Since  we  are  interested  here  in  the  intensity,  we  set  t  =  0  in  (3.4c), 
e.g.,  pL(0)  =  1.  Moreover,  since  we  have  in  reality  a  sum  of  such  inde¬ 
pendent  solitons  traveling  in  all  directions  on  the  broken,  large-scale 
surface,  we  evaluate  (3.4b)  and  write  finally  in  this  isotropic  case 


M^(^r,0)  =  Oj  e 


-  2  2  _  N  _2 


2,.  .2 


;  °s =  5r  cs  **  ;  =  l/2t  Ar  =  Ax  +Ay 


(3.5) 


where  N  average  number  of  solitons  overlapping  at  any  instant  and  posi¬ 
tions.  Not  surprisingly,  from  a  gaussian  waveform  we  obtain  a  gaussian 
second-moment  function.  (Note  the  "d.c."  component  embodied  in  (3.5), 
since  <^>L  >  0.) 

From  (7.36),  (7.41)  we  obtain  the  corresponding  two-dimensional  wave 
number  spectrum 

00  1V£r  ,  ,  -[(2k  *c)2+(2k  *e)Z]/2 

Ws(k.| 0)  =  J  M$(&r,0)e  d(£r)  =  2tt*so2  e  ,  ox  S  oy  S 

(3.6a) 

which  reduces  here,  cf.  (3.2),  specifically  to 


W$(0,2ko  sineQT|0)  =  2 *o|*|  e 


2,2  -2<kolS  s1"  V 


(3.6b) 


Our  next  step  is  to  note  from  (7.51b)  that 


w$^~l°^Mid  =  2itWS^I°^McD+G 


•2,  («Me'2(VsSin6'>T)2),  (3.7) 


so  that,  writing  NgS-jnc  ;  1(*NqS-inc  ’  we  convert  the  second  (non¬ 
vanishing)  term  of  (3.1)  to  the  backscatter  cross  section  attributable 
to  the  soliton  ripples  in  the  McD-G  normalization: 


S-inc 


McD+G 


=  2  c(% T)  s2  *2  (41  e'2(Vs  S1"  e«T)' 

u  o  g5-inc  oT  o  s  5 


.1 


r 


It  is  this  result  we  next  examine  numerically  for  our  specific  example 
above. 

£.  A  Numerical  Example 

Preliminary  (though  sparse)  observations  [14],  [15],  [17],  [18]  sug¬ 
gest  the  following  ranges  of  values  for  2L,  o^: 


L  =  1. 6-2.0  cm.  :  L  =  2£g  =  /2 


L$)  =  0. 8-1.0  cm.  ^ 

J  .'.L-  =  1.13-1.41  cm. 


o $  =  0. 2-1.0  cm.  ;  -  4«  10-6  m2-10-4  m2; 

2  2 

here  L<.  is  the  "correlation  distance,"  where  Ar  =  is  such  that 
M<.  =  e~*  of  its  maximum  value. 

For  our  example  we  select  as  representative  values: 


L  =  1.8  cm.;.*.  =  0.90  cm.  ;  os  =  0. 2-1.0  cm. 


2  •  2  • 
Moreover,  at  0qT  =80°,  5-1,  and  for  the  water/air  interface,  =  1. 

Combining  (3.2a,d,e)  in  (3.8)  thus  gives  us,  with  (in  meters) 


,2  -2(ko£S  sin  V  . 


=  (0.90  *  10  )e 


_2  2  -2( (83.78) (0.90  •  10”^ ) ( 0 . 985 ) > 


=  (8.1  •  10"3)e 


-2  (a  743)' 


=  8.1  •  10  (e 


-5.-1.104, 


■  (8.1  *  10”5)(0.332)  =  2.69  *  10'5  m2,  (3.1 


the  following  expressions  for  the  soliton  scatter  cross  section,  since 
2/u  -  0.637: 


a<°> 

S-mc 


S (0.637) (4.93  •  107 ) ( 1 . 87  •  10'2)(2.69  •  10'5).oc 

2/tt  k4 


wS|McD+G 


2 

aS 

HI 

■ 

— 

■ 

.(0)  (Ah\.  with  Kirchoff 
°S  '  ' ’  corrections=+2.5  db 

4  • 

nr6 

6.32  •  10'5 

-42.0 

-39.5 

HI 

9  • 

10'6 

1.42  •  10"4 

-38.5 

-36.0 

0.4 

1.6- 

10'5 

2.53  •  10'4 

-36.0 

-33.5 

0.5 

2.5- 

10-3 

3.95  •  10'4 

-34.0 

-31.5 

0.6 

3.6- 

10~5 

5.69  •  10"4 

-32.4 

-29.9 

0.8 

6.4- 

10'5 

1.01  •  10-3 

-30.1 

-27.5 

1.0 

1.0 

10'4 

1.58  10‘3 

-28.0 

-25.5 

[The  Kirchoff  correction  used  here  is  taken  from  Fig.  13  of  [4];  see  the 
discussion  in  Section  V  of  [4].]  This  correction  factor  arises  because  of 
the  failure  of  the  "flatness"  condition  (i.e.,  large  radius  of  curvature, 
b0k0p  >  1)  as  grazing  angles  (0o-j-  >  60°)  are  approached.] 

Experimental  results  for  the  frequency  (f  *  20  kHz),  grazing  angle 
(eoT  =  80°),  and  wind  speed  (v  *  10  m/sec  =  20  knots)  chosen  here  (cf. 

A  above)  are  cited  in  Table  3.2  below: 


Table  3.2:  Experimental  Results*  [at  f  =  20  kHz;  6qT  =  80°;  v  =  10  m/sec] 


Source 

s(0) 

back-incoh 

a 

Galubin  [36];  also  Fig.  1.26  of  [2] 

-30  db  (19 

knots) 

1. 

Lilly  and  McConnell  [37] 

-29  db  (21 

knots) 

3. 

Hoover  and  Kaprocki  [38] 

-27.5  db  (20.5  knots) 

Roderick  [6],  at  0Qy  =  81° 

-30(±l)db(20 

knots) 

♦See  Fig.  4  of  [4]. 


On  the  assumption  that  the  data  of  Table  3.2  were  obtained  in  an 
essentially  bubble-free  environment  (there  was  no  "ground-truth"  to  es¬ 
tablish  this  fact  one  way  or  the  other),  we  see  on  comparison  with 


Table  3.1  that  for  values  of  >  0.4  an  acceptable  agreement  (within  a 
few  db)  between  theory  and  experiment  is  obtained  when  the  scattering 
mechanism  is  the  soli  ton  ripples.  [As  shown  below  in  Section  3.2, 


scattering  from  the  capillary  component  of  the  single  wave  surface  is 
o(10+  db)  smaller.] 

It  has  also  been  observed  that  the  scatter  cross  section  (for  com¬ 
paratively  small  grazing  angles  (0qT  >  60°))  increases  with  frequency 
(Fig.  1.26  of  [2],  and  [3]).  Let  us  consider  the  cases  f  =  10  kHz, 

20  kHz  and  compare  with  Galybin's  data  (Fig.  1.26  of  [2]).  From  (3.8) 
we  write 


dS-inc(@  20  kHz> 


dS-inc<9  10  kHz> 


z4.e- 1.104+0.276 


24  e-0.828 


7.0  =  8.4  db.  (3.13) 


This  compares  acceptably  with  a  '<7.5  db  estimated  difference  from  Fig. 
1.26,  [2],  and  with  Roderick's  recent  results  [6], 


3.2  Backscatter  from  Capillary  Waves 

It  remains  to  examine  the  effects  of  the  capillary  component, 
dc-inc’  in  the  total  backscatter  cross  section  (3.1)  at  these  high  fre¬ 
quencies  and  small  grazing  angles.  Here  the  physical  mechanism  is  dif¬ 
ferent:  the  capillary  extension  of  the  gravity  wave  surface  is  directional 
and  dispersive,  with  dispersion  governed  by  surface  tension  forces,  unlike 
the  soli  ton-ripple  surface  above,  which  rides  upon  the  former.  In  what 
follows  we  consider  very  briefly  several  spectral  models  of  the  single 


wave  surface  (g  +  c).  From  these  we  estimate  o  ,  the  rms  capillary  wave 


(Or 

height,  and  the  back-scatter  cross  section  in  (3.1) 


/V.  Phillips'  Spectrum  [39] 

From  Sec.  4.5  of  [39]  we  have  Phillips  capillary  wave  number  spectrum. 


*(k)c  =  -^4,  a0-w/2  <  a  <  w/2+aQ;  B'  *  1.5  •  10  2; 

trk  ~ 


(3.14) 


cf.  Eq.  4.5.9,  [39].  Since  (2ir)2y(Jk)  =  Wa(ii0)Mid  *  2*a(k|0)|fcD+G1 


cf.  (7.51),  we  see  that  <3^  in  (3.1)  becomes 


where  we  have  used  (3.2f),  and  Ng?jnc  =  16Ng^»  with  (3.14),  (3.2),  etc. 
This  back-scatter  cross  section  is  much  smaller  than  our  soli  ton-ripple 
results  (Table  3.1)  and  the  empirical  data  (Table  3.2).  Furthermore, 
it  is  independent  of  frequency,  nor  does  it  take  into  account  the  disper¬ 
sive  character  of  the  capillary  waves,  as  used  here.  However,  it  agrees 
quite  well  with  Bahar  et  al.'s  results  (-55  db)  in  Fig.  lib,  [25],  for 
the  case  of  horizontal  polarization  (which  is  the  EM  analogue  of  our 
acoustical  case  here).  [We  must  subtract  11  db  (=  4ir)  from  in 

Fig.  lib,  because  of  Bahar's  definition  of  cf.  (2.40)  above.] 

J3.  Brown's  Spectrum  [40];  (Bahar  et  al.,  [25]) 

Since  what  primarily  is  important  in  the  wave  number  spectra  appearing 

in  the  scattering  cross  sections,  cf.  (3.1),  is  not  specific  spectral 

o 

shape  but  rather  the  intensity  (a  c)  and  spectral  spread  (k  i  <■),  we 
shall  use  Brown's  spectrum  (Eqs.  (23),  (24)  of  [25];  [40])  to  determine 
a2  from  Eq.  (26)  of  [25]  and  apply  the  result  to  (3.8),  with  = 

C  ri 

1.27  cm  again,  cf.  (3.10),  as  a  reasonable  correlation  distance  (e  )  for 
these  capillary  waves.  Taking  kd  =  1.0  rad  cm"1  as  the  wave  number  at 
which  the  (continuous)  gravity-capillary  wave  surface  splits  into  "high"- 
and  "low"-frequency  components  ([25],  discussion  and  Table  1  therein), 
we  obtain 


a2  •  _B_  _  4.6  ■  10  J 

c  2k2  2(10+2  rad  m"1)2 


2.3 


10"7  m2: 


4.8  •  10 


cm 


(3.16) 


Applying  (3.16)  with  l  -  1.27  cm,  f_  =  20  kHz,  (v  =  10  m/sec)  in  (3.8), 

c(0l  0 

with  (3.2f)  replacing  N^'  since  the  wave  surface  is  a  single  surface 
[e.g.,  the  capillary  structure  is  part  of  the  wave  surface,  not  a 
separate  entity  riding  on  the  gravity-wave  component  (cf.  Sec.  2.1)], 
we  get  from  (3.12a) 


4-)nc  =  2-10  4  =  4-86 


10 


-7  _ 


■63.1  db, 


(3.17) 


which  is  clearly  too  small. 


Alternatively,  using  Brown's  spectrum,  Wy,  Eq.  (23)  of  [25], 
directly  in  (3.1),  with  (3.2f)  and  the  conversion  relation  WMid  = 


B+B=Brown 


,  cf.  (7.51c),  we  get 


a<°L  =  —  =  (0.637)(4.6  •  10"3)  (2.50  •  10"3)  =  7.3  •  10' 

c-inc  it  gc 


-51.4  db, 
(3.18) 


which  is  -40.4  db  in  Bahar  et  al.'s  [25]  definition  of  scatter  cross 
section,  which  in  turn  is  4"  *  6|^oh  *iere»  (2. 6).  Comparing  this 
with  o^HH\  Fig.  lib,  [25],  we  see  that  (3.18)  is  about  3.5  db  larger 
than  the  cross  section  for  horizontal  polarization, ,  the  radar  ana¬ 
logue  of  the  acoustic  scatter  cross  section.  (The  difference  is  attribut¬ 
able  to  the  particular  "tilt-factor"  used  here.) 

In  any  case,  these  models  give  noticeably  too  small  results  (in  this 
second  instance  independent  of  frequency),  and  all  lack  the  spectral 
character  appropriate  to  the  particular  dispersive  nature  of  capillary 
waves,  cf.  Eqs.  (7.52)  et  seq. 

C.  McDaniel's  and  Gorman's  Results  [3],  [4] 

In  two  recent  papers,  [3],  [4],  McDaniel  and  Gorman  have  obtained 
back  scatter  cross  sections,  which  are  -42  db  (including  the  Kirchoff 
corrections)  for  our  example  above  (.A),  cf.  Fig.  4  of  [4  ],  and  are 
noticeably  below  the  empirical  results,  cf.  Table  3.2.  Their  very  ex¬ 
tensive  analysis  takes  into  account  the  dispersive  character  of  the  small- 
scale,  or  capillary  component,  but  does  not  include  the  suggested  inde¬ 
pendent  soli  ton-ripple  mechanism.  Moreover,  they  also  show  that 
diffraction  effects  (including  multiple  scatter)  are  quite  negligible 
(Sec.  VI,  [  4]),  as  we  note  below  in  Sec.  2.3  also.  We  refer  the  reader 
to  these  papers  [3],  [4]  for  details;  (see  also  the  discussions  in  Secs. 

1,  4  here). 


3.3  Backscatter  from  Diffraction  Terms  (k  >  1) 

Applying  (8.17)  to  our  definition  (2.6)  we  see  at  once  that  the  back- 
(and,  in  fact,  omni-directional)  scatter  cross  section  for  the  diffuse 
diffraction  component  (k  =  1)  in  our  theory  [41 ]- [43]  becomes  (cf.  Sec. 

8.6  ff.) 


41 


y 


AIL  ■  T7T  <«glJgl>  *  «clJcl>  *  <5IlJSl>}  • 


(3.19) 


' Ja  1  “  l^xx^yy  "  ^xyU’ 


which  is  put  to  the  same  scale  as  the  principal  contribution  ^coh’ 
(3.8),  for  comparison  with  the  data  of  Table  3.2  and  the  example  of  3.1£, 
by  the  relation  2ird{J>  h  =  3^J+G,  cf.  (3.7). 

Some  simple  numerical  estimates  of  (3.19)  may  be  obtained  by  re- 

p  p  | 

placing  <C  |J  |>  by  <c  >|j  ’|.  We  suggest  that 


r.y, 

C2=  O(10_1  m2); 

|JJ  •  O(10‘6  ro"2) 

<£  Ij  L>  - 

V 

9 

Q 

9  9 

m 

m 

O(10'5  m2l; 

|jcI  =  odo0-io_1) 

•••  <i  ij  ic>  - 

(3.20a) 


;io‘5-io'6) 


0(2  -10“5  m2);  |j$  |  =  OdO'1)  |j  l$>  -  0(lO_b)  =  O(-60  db) ; 


=  O(-50  to  -60  db); 

(3.20b) 

_  -  n/_cn  ,ik\. 


(3.20c) 


where  in  (3.20b,c)  the  curvature  is  comparatively  large.  Since 
o 

( 16 tt  )  =  -21.2  db,  we  see  that  even  if  these  scatter  terms  are  10  db,  or 
even  20  db,  too  large,  they  are  ignorable  compared  to  the  main  capillary 
(and  gravity)  contributions  3^)f  3^)  unless  there  is  a  special  align¬ 
ment  of  the  former  which  renders  it  essentially  zero,  e.g.,  at  <f>0j  =  V2, 
and  a0  =  o,  so  that  the  mean  direction  of  the  gravity-capillary  wave 
surface  is  directed  at  right  angles  to  the  beams  in  backscatter,  cf. 

Fig.  2.1a.  The  main  contributions  of  the  soliton  ripples,  3^,  however, 
dominates  the  above  in  any  case,  and  is  essentially  directional  as 
well. 

We  may  expect  the  multiple-scatter  terms,  3^2\  in  the  scattering 
cross  section  (2.8).  to  be  even  smaller  than  the  diffuse  scatter  terms 
(k  *  1)  above.  This  corresponds  to  their  negligible  contributions  as 
computed  by  perturbation  theory  (when  the  slopes  of  the  capillary  and 


wa 
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As  we  have  already  noted  in  Section  1  above,  our  principal  results 
are  fundamentally  theoretical,  combined  with  selected  numerical  compari¬ 
sons  of  our  new  model  with  earlier  theory  and  experiment.  We  summarize 
below  the  general  new  features  of  this  study  (Sec.  4.1  ff.)  and  suggest 
needed  further  theoretical  and  experimental  efforts  (Sec.  4.2  ff.). 

4.1  Principal  Results 

The  new  features  of  our  work  here  are: 

(1).  The  identification  of  a  potential  mechanism  for  the  second  component 
of  a  composite  wind-wave  surface  which  may  resolve  the  often  ob¬ 
served  discrepancy  between  theory  and  experiment  at  high  frequencies 
and  small  grazing  angles  without  our  postulating  a  bubble  mechanism 
[3],  [4]. 

This  alternative  mechanism  is  the  ensemble  of  soli  ton  ripples, 
or  hydraulic  jumps,  riding  on  the  main  gravity-capillary  wave  sur¬ 
face,  by  which  the  latter  is  itself  generated  through  the  nonlinear 
action  of  local  near-surface  winds,  or  "cats-paws"  [cf.  Sec.  2.1  and 
[14]-[18].]  These  ripples,  broken  up  into  solitons  on  the  wind- 
driven,  rough  gravity  wave  surface,  travel  with  constant  speed  and 
are  nondispersive,  unlike  the  gravity  and  capillary  waves,  cf.  Secs. 
2.1,  3.1  above,  and  are  directional,  according  to  the  local  near¬ 
surface  wind.  They  disappear  when  the  lecal  wind  (momentarily)  stops, 
but  are  otherwise  generated  much  of  the  time  in  usual  seas. 

Preliminary  numerical  estimates  (Sec.  3.1)  indicate  quantitatively 
that  this  ripple  mechanism  can  account  for  the  previously  observed 
discrepancies  between  theory  and  experiment  at  high  frequencies  and 
small  grazing  angles,  when  there  are  no  bubbles. 


r 


(2).  The  development  of  a  truly  two-component  wave  surface  model,  where 


the  above  ripple  surface  is  essentially  independent  of  the  single 
gravity-capillary  wave  surface  on  which  these  soliton-ripples  ride. 

Analytically,  the  practical,  approximate  theory  requires  three 
scattering  components:  (i),  large-scale  gravity  waves  ("low"-frequen- 
cy);  (ii),  the  capillary  continuation  of  these  gravity  waves  ("high"- 
frequency);  and  (iii),  the  independent  soli  ton  or  hydraulic  “jump" 
ripples,  generated  by  the  local  wind  action. 

The  divisions  of  the  wave  surface  into  components  (i)  and  (ii) 
is  an  insightful  technical  device  [as  Bahar  et  al .  [25]  have  noted, 
for  examplel  while  (iii)  is  a  separate  phenomenon.  The  former  yield 
specular-point  scatter  (i),  which  dominates  at  large  grazing  angles, 
and  (first-order)  Bragg  scatter  (here)  (ii),  while  (iii)  also 
generates  a  first-order  Bragg  scatter,  which,  however,  we  suggest  is 
significantly  0(10+  db)  larger  than  the  capillary  contribution  (ii), 
cf.  Sections  (3.1),  (3.2).  Of  course,  in  those  intervals  when  there 
are  no  local  interacting  winds,  the  capillary  component  dominates, 
unless  there  are  sufficient  near-surface  bubbles.  For  example,  these 
latter  can  be  comparable  to  the  soliton  ripples,  producing  o(3  db) 
increase  in  the  scattering  level,  which  are  probably  not  resolvable 
within  the  accuracy  of  previous  experiments.  See  also  [6]. 

(3).  Verification  of  the  general  analytical  agreement  between  the  various 


rincipal  approximate  results  for  the  scattering  crobS  sections 


obtained  in  the  literature  [cf.  Sec.  2]. 


Thus,  the  Kirchoff-Perturbation  approaches  of  [4],  [21],  [33], 

for  example,  and  the  author's  "Kirchoff-Kirchoff"  method  (cf.  Sec. 

7.3  )  yield  essentially  the  same  form  of  Bragg  scatter  component, 

albeit  with  somewhat  different  "tilt-factors,"  (cf.  (2.37),  (2.39)). 

In  most  instances  these  tilt-factors  are  approximately  equal, 

because  the  rms  slopes  of  the  large-scale  gravity-wave  component  are 

2  -2 

themselves  quite  small,  e.g.,  a  u  =  0(2,3  •  10  ).  An  excep- 

9-w  (n} 

tion,  however,  is  the  tilt-factor,  N^',  associated  with  the  arti¬ 
ficial  split  of  the  gravity-capillary  wave  surface  into  low-  and 


M 
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high-frequency  components:  the  tilt-factor  here  vanishes  with  zero 
grazing  angle,  unlike  the  situations  where  there  is  truly  a  second, 
separate  wave  mechanism  riding  on  the  first  surface,  as  presented 
here. 

Numerical  differences  [cf.  Sec.  3.2]  between  the  earlier  results 
stem  principally  from  the  choice  of  a  wave  spectrum  model,  which  in 
turn  depends  quite  critically  on  what  is  the  underlying  physical 
mechanism. 

Corrections  for  the  well-known  limitations  of  the  Kirchoff  or 
Tangent-Plane  method  are  small  0(1-3  db),  cf.  [4],  Section  VI,  at 
small  grazing  angles,  vis-a-vis  the  general  level  of  scattering  cross 
section  determined  through  these  earlier  models,  as  they  are  also 
when  our  present  wave  surface  model  is  employed. 


(4).  The  observation  that  the  diffraction  terms  (diffuse  scatter  (k  =  1), 
and  multiple-scatter  (k>2),  cf.  (3.20),  are  generally  well  below 
[0(15-20  db,  cf.  Sec.  3.3]  the  main  contributions  of  the  direct 
scatter  components,  e.g.,  the  specular-point  scatter  components 
(important  at  large  grazing  angles)  and  the  Bragg-scatter  terms, 
which  dominate  at  small  grazing  angles  (0qT  >  60°),  and  high-fre¬ 
quencies.  This  agrees  with  recent  results  [4]  using  different  theor¬ 
etical  methods  (i.e.,  perturbation  techniques). 


(5).  From  the  above  it  is  therefore  clear  that  the  large  discrepancies 
between  theory  and  experiment  at  these  critical  regions  of  small 
grazing  angles  and  high  frequencies  are  to  be  explained  by  the  proper 
choice  of  physical  model:  the  problem  is  not  fundamentally  analytical. 

We  emphasize  again  that  a  bubble  mechanism  can  explain  the  observed 
discrepancies  between  theory  and  experiment,  provided,  of  course,  bubbles 
are  actually  present  in  sufficient  numbers.  Moreover,  if  bubbles  are  so 
present,  their  effect  should  add  only  0(3  db  =  factor  2)  to  that  already 
produced  by  the  independent  "ripple"  surface  (on  the  assumption  both 
mechanisms  are  comparable  scatterers).  The  presence  of  bubbles,  however. 


was  unfortunately  not  quantifically  verified  in  most  of  the  experimental 
data  cited  earlier  for  comparison  with  theory,  cf.  Table  3.2  and  refs. 

We  also  emphasize  that  our  proposed  soliton-ripple  mechanism,  al¬ 
though  an  observed  phenomenon,  has  only  been  measured  to  a  limited  extent, 
particularly  on  gravity  wave  surfaces:  our  numerical  estimates  have  been 
inferred  from  very  limited  data  so  far.  Thus,  while  the  reliability  of 
the  analytical  procedures  is  not  in  question  [cf.  (3)  above],  it  has  not 
yet  been  established  that  the  soliton-ripple  mechanism  fully  explains  the 
observed  scattering  levels  in  the  absence  of  bubbles.  Further  experimental 
work  is  required  [cf.  Sec.  4.2  ff.],  to  obtain  the  needed  dimensions 
(height  and  length)  of  this  mechanism.  [We  remark  that  whenever  inde¬ 
pendent  measurements  verify  this  mechanism  quantitatively,  we  can  then 
apply  acoustic  (back-)  scattering  techniques  like  those  described  in  de¬ 
tail  here  to  obtain  the  desired  scale  and  size  of  these  ripple  effects 
in  subsequent  applications.] 

Other  new  features  of  our  analysis  include:  (1)  explicit  development 
of  the  role  of  gaussian  and  omni-directional  beam  patterns;  (2)  general 
second-moment  and  covariance  functions,  from  which  we  shall  subsequently 
determine  wave  surface  spectra  and  related  (i.e.,  doppler)  effects; 

(3)  general  narrow-band  signals  and  absorptive  phenomena;  (4)  numerical 
estimates  of  diffraction  contributions;  and  (5)  an  extensive  discussion 
of  the  various  approximations,  and  their  conditions,  which  underly  the 
analytical  results. 

Finally,  we  note  that  our  general  model  here  can  be  applied  at  once 
to  such  special  cases  which  may  occur  physically  as: 

(i) .  no  soli  ton  ripples,  just  the  gravity-capillary  wave  surface 

(no  local  winds):  one  drops  the  "S-terms"  in  our  results. 

(ii) .  soli  ton-ripples  and  bubbles:  one  may  use  the  additional 

bubble  contributions,  for  example,  as  treated  in  [3].  Their 
effects,  of  course,  will  depend  strongly  on  the  bubble  den¬ 
sity,  whether  or  not  they  are  ignorable  to,  comparable  to, 
or  dominate  the  soli  ton  ripples. 


The  following  is  a  necessarily  incomplete  list  of  suggested  next 
steps  in  the  ongoing  treatment  of  surface  scattering.  It  is  convenient 
to  separate  the  list  into  Experimental  and  Theoretical  topics,  although 
the  two  are  closely  related  here.  We  have: 

A_.  Some  Experimental  Topics 

(1) .  Measure  directly  the  rms  height  and  the  mean  duration  of  the 

soli  ton-ripples,  under  various  near-surface  conditions,  [14]- 
[18],  including  no  to  full  gravity-wave  development; 

(2) .  Mean  (back-)  scatter  intensities  under  (1),  simultaneously, 

with  particular  attention  to  the  possible  presence  of  bubbles; 

(3) .  "Forward"  and  bistatic  scatter  generally,  also  under  (1),  (2) 

above; 

(4) .  Develop  acoustic  techniques,  based  on  the  associated  theory 

here  and  below  (8),  for  measuring  the  relevant  parameters  of 
the  small-scale  scatter  components. 


JL-  Theoretical  Extensions 

(1) .  Apply  results  of  (1)A  above  to  the  analysis  developed  here; 

(2) .  Calculate  representative  (backscatter)  cross  sections  for  all 

scattering  angles  (0  <  0Qy  <  tt/2) ,  and  selected  surface  ripple 
parameters; 

(3) .  Apply  Bahar's  general  "full-wave"  approach  [44]  in  detail  to 

relate  the  former  to  the  approximate  solution  obtained  here. 

(4) .  Determine  the  doppler  spectrum  of  the  composite  wave  surface 

in  our  general  model  and  compare  with  experiment  [6]; 

(5) .  Examine  the  contributions  of  higher-order  Bragg  scatter; 

(6) .  Develop  further,  and  apply  to  quantitative  estimates,  the  con¬ 

tributions  of  the  diffraction  terms  vis-a-vis  higher-order 
approximations  in  the  direct  scatter  theory; 

(7) .  Generalization  of  the  results  to  anisotropic  surfaces; 

(8) .  Examination  of  various  corrections  to  the  Kirchoff  or  Tangent 

Plane  approaches. 

Finally,  we  stress  once  more  the  need  for  carefully  controlled  experi¬ 
mentation,  where  "ground  truth"  regarding  surface  phenomena  is  fully  ob¬ 
tained,  and  where  the  subtleties  of  the  associated  signal  processing  are 
both  understood  and  applied  [6]. 


Part. 1. 1..  Analytical  Models 


Here  we  present  the  analytical  features  of  random  wave-surface  scat¬ 
tering  when  the  incident  radiation  is  acoustic  propagation  from  a  direc¬ 
tional  source.  A  directional  receiver  is  employed,  in  the  manner  of 
Figure  2.2  above.  While  much  of  the  present  theory  is  well-known  [1], 

[2],  [8],  [9],  [21],  [22],  [33],  we  shall  also  include  here  a  number  of  new 
features  and  results.  Our  treatment,  however,  is  intentionally  concise:  we 
refer  tne  reader  to  other  references  (see  Section  1)  for  many  of  the  details, 
in  the  interest  of  brevity  and  attention  to  the  main  problem  at  hand,  viz., 
the  determination  of  back-  and  forward-scatter  intensities,  along  with  their 
attendant  assumptions  and  approximations. 

5.1  The  Received  Scattered  Field,  I:  Formulation 

For  the  total  scattered  field  a(R,t)  at  a  point  P(R,t)  in  a  linear, 
inhomogeneous,  and  possibly  random  medium,  or  on  its  boundaries,  we  can 
write  generally  [12],  [43],  [41],  in  operator  form 


S  ^  =  MO, 


(5.1) 


where  fcM^-G,)]  is  the  homogeneous  or  unscattered  field  component  and 
the  brackets  {  }  denote  the  ensemble  of  such  (random)  media.  Here  -Gy  is 
the  source  distribution  (which  includes  the  sensor  array),  and  is  the 
integral  operator  whose  kernel  is  the  green's  function  g(R,t  | R_'  ,t'  )eo  for  an 
infinite,  unbounded  medium  (except  for  source  regions),  which  is  also  homo¬ 
geneous  (but  not  necessarily  lossless).  The  operator,  is  the  field- 
renormalization  operator,  where  Q,  cf.  (5.1),  embodies  the  scattering  mech¬ 
anisms  appropriate  to  the  physical  situation.  Here,  for  example,  we  have 

Q  -  Qs  +  Qv  ,  (5.2) 

where  represents  the  scattering  effects  of  the  air-water  boundary,  or 
wave  surface,  while  Qy  describes  any  volume  inhomogeneities.  The  specific 


form  of  Mm  is,  of  course,  determined  by  the  (usually  differential)  propa¬ 
gation  equation  for  the  medium,  without  inhomogeneities. 

Next,  we  shall  make  the  following  assumptions: 

I.  The  medium  supports  an  (extended)  Helmholtz  equation  of  propa¬ 
gation,  e.g.. 


(V2-i^)a=-G 

CZQ  3t2 


(5.3a 


with  cQ  =  speed  of  wavefront  propagation  in  the  medium. 

II.  The  medium  (except  for  boundaries,  here  the  surface)  is  therefore 
homogeneous 


**  Q v  =  0- 


(5.3t 


III.  Only  local  surface  scattering  interactions  are  significant,  so 

that  Eq.  (5.1)  becomes  on  developing  the  Perturbation  Theoretical 
Series  (PTS)  [42] 


..y  „  ;  «„  +  N„Qs“h  • 


(5.3c 


a  =  oiH  +  cij  ;  dj  = 


(5.3c 


The  approximation  (5.3c),  which  neglects  coupling  of  the  scattered  field 
to  the  scattering  elements  (embodied  in  Q^),  e.g.,  terms  o([MooQ]v  ')  or 
higher  are  dropped,  is  a  form  of  Born  approximation.  This  appears  to  be 
eminently  acceptable  in  practice,  since  scattering  from  wave  to  wave  over 
more  than  neighboring  wave  structures  is  negligible  provided,  for  c  = 
cQ  +  Cj(z),  that  vc  -  0:  only  "nearest  neighbor"  wave  surfaces  will 
support  multiple  scatter  (diffraction)  effects.  Consequently,  we  call 


50 


Here  f  is  the  central  frequency  of  the  original  (narrow-band)  signal, 
and  2r,  are  (inward)  unit  vectors,  with  j^R  along  the  main  axis  of  the 
(complex)  receiving  beam,  Q^,  cf.  Fig.  5.1.  (We  note  also  that  the  beam 
patterns  are  functions  of  frequency,  generally.  Sec.  6.1).  For  the  omni¬ 
directional  receiver  used  in  the  experiments  [6],  =  Q^(f0),  a  complex 

constant,  [cf.  Eq.  (6.14)  ff.]. 

5.2  The  Unscattered  Field 

For  the  Helmholtz  media  (obeying  (5.3a))  it  is  well-known  that 


M„(R.t|R',f)  =  -  J  df  J(  )R1  t.  dR'  gjR.tlR'.f); 

with  gM  =  <S(t-t'-p/c  )/4irp  ,  p=|R'-R| 

so  that  the  homogeneous,  i.e.,  unscattered  field  component  <*H,  (5.5), 
becomes  at  once  the  familiar  result  [41],  [42] 


«H(R,t)  =  /  GT(t-p/CQ,RJ  )  47(p  ,  Gy  C  vt 
VT 


(5.8) 


(5.9) 


where  VT  is  the  source  domain  and  R  lies  outside  VT.  In  the  far-field 
T  —  T 

we  get  in  straightforward  fashion  [34]  for  narrow-band  signals: 


s(t-iT-R/c„)ds 


;  4&  Ur,  Oi(dr-lvTls0/*1 '  S0IS1n(s))e 


(iT  -  R/lRt ;  ■  R) 

where  now  CL  Is  the  generalized  beam  pattern  [12]  (with  R'  .  £  ff.) 

s  (i  • i  )*C/c 

QT  ■  I  Vi-so>sin<sl<>e  0  '  °  &  •  <s0  =  2'1fo>- 

VT 

in  which  is  a  "steering"  (unit)  vector,  cf.  (5.7a).  Here 
=  Vin^LP  =  Ao-t{-(t*i)} 


(5.10) 


(5.10a) 


(5.11) 


(5.12) 


swr** 


I 


is  the  amplitude  spectrum  of  the  (here  narrow-band)  signal  S.  (t)  =  A  s ( t ) ; 

oo  —  p  +  in  o 

=  /  e  (  )dt,  Re(s)2>0).  The  quantity  Aq  is  the  (peak)  amplitude 
of  Sin(t).  e.g.  <S?n>  =  <s^>,  <s^>  =  1,  where  <>  =  time  average  (as 


T  -*  “),  or  over  the  finite  interval  (tQ,  tQ+T  <  =°)  if  S.jn(t)  is  of  finite 
duration. 

Like  AR,  (5.6),  Ay  is  the  aperture  weighting  associated  with  the 
transmitter,  at  each  element  d?  where  the  signal  is  applied.  In  the  gen¬ 
eral  case  different  signals  may  be  applied  at  different  locations  (s)  on 
the  aperture:  hence,  S.  =  AS.  (sjs),  cf.  (5.12).  However,  we  shall 
assume  with  little  loss  of  generality  that  the  same  signal  drives  each  ele¬ 
ment,  so  that  Sin(s|£)  =  Sin(s).  Then  (5.11)  reduces  to 


■  Sin<s>  f  VW  6  0  "T 
'T 


so(iT-joT>-i/c 


dC 


Sin^s%^T"^oT;fo^  =  -Tfo/co’  etC-^ 


(5.13 

(5.13 


where  now  (Xj.  is  the  (complex)  beam  pattern  of  the  transmitting  aperture, 
in  these  far-field  situations.*  Applying  (5.13)  to  (5.10)  then  allows 
us  to  write  for  the  (as  yet)  unscattered  field  the  not  unexpected  result 


(5.14 


subject  to  the  uniform  drive  condition  S - n ( s | ^ )  =  Sin(s),  with  dy  given 
by  (5.13b)  as  defined  in  (5.13a). 


5.3  Absorption  and  Doppler 

Although  so  far  we  have  assumed  that  the  medium  of  propagation  is 
lossless,  in  reality  there  is  some  attenuation,  mainly  because  of  molecular 
absorption  at  the  frequencies  employed  here  (<40  kHz).  In  general. 


*In  the  "near-field"  or  Fresnel  regions  the  "beam  pattern"  becomes 

range-dependent,  as  a  consequence  of  the  quadratic  terms  in  the  expansion 

of  p  =  I  R-C |  ,  (R1  -*■  £).  To  avoid  this,  a  useful  far-field  condition  is 
—  —  —  **» 

1tl"max^xo<<R’  w*1ere  *-max  1S  the  "*ar9est  dimension  of  the  aperture  and  XQ 
(=c0/fQ)  is  the  wavelength  of  the  (central)  frequency  of  the  postulated 
narrow-band  signals. 


absorption  affects  the  beam  patterns  and  waveforms.  But  with  narrow-band 
excitations  the  general  relations  simplify  greatly,  reducing  to  a  simple 
attenuation  factor. 

To  see  this  we  use  the  transform  YQ  _  =  ^{g^}  of  the  green's  func¬ 
tion  solution  to  the  basic  propagation  equation  for  relaxation  absorption 
[42],  viz. 


c0  3t 


5tt’  JJR' 


so  that  Yq  obeys 


[(1  *  t  s)v2  -  sZ/c l  ]Y'  .  •  S  ,  ;  Y  (jg.slR'.f).  =  Y'(R,s|R').e‘ 


0  J  0, 


For  this  it  is  readily  found  ([41],  I)  that 


4irp  /I  +  tqxS 


*  Tox  =  “x/co  ‘ 


For  typical  oceans  the  relaxation  time  tqx  of  the  medium  is  tqx  =  t^q 

-  o(10-6  secs),  so  that  for  weak  absorption,  e.g.,  (t  Js|)2<<  1,  f  v  4 

ox  max 

(=  |s|/2tt)  is  0(40  kHz).  For  frequencies  less  than  0(40  kHz)  the  medium 
is  essentially  purely  absorptive  and  nondispersive.  In  fact,  we  can  re¬ 
write  (5.16)  as 


jr-  -  st* 
co 


{1'(1'Toxs)"2}ps/co 


✓1  +  t„ys 
ox 


^Yo,“^Helmhol 


u-l' 


(tox/2co)pS‘ 


-0)nPTnx/2cn 
°  ”  °  *5, 


since  for  narrow  band  signals  s  -*•  2itif0,  where “lylg^  =  Eq.  (5.8)}  gives 
(Yq  00)Helm’  the  ^rst  ^actor  of  (5.17a).  For  fQ  =  20  kHz,  ✓1+t()x| sj  = 
1.06  =  1,  and  for  f  =  40  kHz,  /1+tqx|s|  =  1.12  =  1,  also:  the  s-dependent 
denominator  is  effectively  unity.  Accordingly,  for  narrow-band  signals 
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and  positions  in  the  far-field  of  the  source  we  simply  scale  the  signal 


amplitude  A  by  the  factor  exp(-a^R),  viz. 


a  =  t  /2c 
ox'  o 


0(3.7  -10"10  sec2/meter)  .  (5.18) 


Equation  (5.18)  shows  the  well-known  result  of  attenuation  exponentially 
dependent  on  distance  and  on  the  square  of  the  frequency  for  media  where 
relaxation  absorption  governs  propagation,  in  the  manner  of  Eq.  (5.15) 
above. 

The  effects  of  doppler  are,  however,  more  complex  than  the  ("low- 
frequency")  attenuation.  Doppler  produces  a  frequency  modulation  of  the 
original  signal.  The  principal  mechanisms  for  this  in  underwater  studies 
are  platform  movement  and  the  motion  of  the  ocean  wave-surface.  As  long 


as  the  doppler  velocities  are  comparatively  small  [e.g.,  v/cQ  <  0.10,  say]. 


a  practical  theory  is  possible,  cf.  [34].  In  our  present  study,  cf.  Fig. 
2.2,  both  the  transmitting  and  receiving  platforms  are  effectively  sta¬ 
tionary  in  space,  so  that  only  the  wave  motion  modulates  the  incident 
radiation  in  the  course  of  scattering.  The  effect  of  this  is  to  introduce 
a  time-delay  AtdS,  in  the  time-variable  quantities  which  appear  in  the 
expressions  for  the  received,  scattered  field,  as  noted  below  in  Eqs. 

5.22.  This  time  delay  is  a  function  of  a  time-variable  position,  depen¬ 
dent  on  the  moving  wave  surface. 


5.4  The  Received  Scattered  Field,  II:  Canonical  Forms  and  Geometry 

Before  we  go  on  to  consider  explicit  wave-surface  models,  which  will 
dictate  the  specific  nature  of  our  attack  on  the  problem  [cf.  Sec.  6  ff.], 
let  us  obtain  a  canonical  series  development  for  the  received,  surface- 
scattered  field  represented  here  by  (5.5). 

This  is  done  by  observing  that  the  field  renormal ization  operator, 
n,  (=  MJ)^),  cf.  (5.1),  can  be  expanded  in  some  suitable  hierarchy  of 
(linear)  operators  of  increasing  complexity,  as  represented  by  the  sum 


■n  ■  j. 

«  £=Q 


(5.19) 


since  the  scattering  process  is  linear.  The  quantity  a^  is  an  expansion 
coefficient,  whose  specific  form  is  determined  jointly  by  the  physical 
nature  of  the  scattering  surface  and  by  the  particular  method  with  which 
a  specific  expansion  is  to  be  carried  out. 

For  example,  in  the  usual  perturbation  technique  (Perturbation 
Method  h  P.M.)  here,  the  boundary,  e.g.,  scattering  surface  z,  is  per¬ 
turbed  or  expanded  with  respect  to  some  appropriate  reference  surface  S 
[4].  In  recent  treatments  [33],  [4]  where  a  two-scale  wave  surface  is 
postulated,  the  reference  surface  is  chosen  to  be  the  large-scale,  or 
gravity-wave  component.  The  perturbation  is  then  embodied  in  the  expan¬ 
sion  of  inT)^,  where  is  defined  on  £  and  expanded  on  S,  cf.  Sec. 

I  of  [4],  while  Q<-  is  likewise  defined  now*  on  S.  Thus,  (5.19)  becomes 

no  =  I  a  A.  Qs.  (5.19a) 

£=0  b 

£  £, 

The  expansion  coefficients,  a0,  become  now  a„  =  a  =  (k  a)  ,  with  the  con- 
2  Z  *  x  o  £ 

dition  kQa^  <<  i,  a ^  =  <e*>,  which  is  required  for  practical  convergence 

of  series  like  (5.19).  In  addition,  the  radii  of  curvature  of  the  small- 

scale  surface  component  must  be  large  compared  to  Ucl ,  the  magnitude  of 

the  elevations  of  the  small-scale  surface  component  on  S,  so  that  the 

expansion  (of  Mm)  is  unique.  This  is  equivalent  to  requiring  that  the 

slopes  of  cc  be  small,  also,  a  condition,  however,  not  generally  met 

here  [cf.  Remarks,  Sec.  l]. 

Accordingly,  we  employ  an  alternate  method  of  expansion,  based  on 
the  structure  of  the  (local)  surface  scatter  operator  Q<.  (on  z),  rather 
than  on  the  expanded  projection  of  the  (non-local)  integral  operator 
on  S.  The  expansion  coefficients,  a^,  are  now  based  on  certain  statis¬ 
tics  of  the  numbers  of  different  orders  of  (multiple)  scattering  inter¬ 
actions.  This  type  of  expansion  turns  out  to  be  essential ly  independent 
of  both  local  elevation  and  slope  conditions,  cf.  above.  Instead  of 
(5.19a)  we  have 


♦In  a  one-scale  wave  surface,  S  =  SQ  =  <s>  =  0:  the  reference  sur 
face  is  the  equilibrium  {=  mean)  wave  surface. 


CO 


00 


(5.1 


=  I  bA^S*3  ’  or  QS  =  £  b£QS* 

H=Q  *  *  3  fc=Q  %  5 


The  kernels,  or  weightings,  of  these  surface  scatter  operators  Q^-l 

"Till  ■* 

are  given  by  Q£  J CS.' , t * ) ,  where  .R'  €  e,  the  (composite)  wave  surface 

and  t'  is  a  doppler-delayed  epoch.  As  we  shall  see  presently,  cf.  Sec. 

7,  the  lower  orders  (*.=0,1)  generally  represent  direct  scattering  com¬ 
ponents  here,  while  the  higher  orders  (*  >  2)  embody  the  multiple-scatter, 
or  diffraction  contributions. 

The  integral  operator  M^,  cf.  (5.8),  applied  to  Q<-,  which  is  itself 
applied  to  aH  (5.14)  on  the  actual  wave  surface,  E,  in  the  same  way  aH 
was  obtained,  gives  us  in  the  far-field  of  the  surface,  the  scattered 
field  a » .  This  becomes  with  the  help  of  (5.7)  in  (5.5)  the  following 
general  (complex)  result: 


:.T(t)  =  l  X^(t) 
1  *=0 


where,  canonically,  we  find  that,  along  with  (5.18)  for  absorption. 


Ao  f  n  Sin^2^i ^FS 
1 


If  lest 
•|ro,e  27T 


(5.21a) 


FS  }  =  L  ^ -R -oR’V^-T‘^oT;VQS 
o 


-sT0-s(iT-iR)-r'/c0  -a^(iRTHl+|RR+£l )  dSc 


(4it  )^RTRRn2 


(5.21b) 


where  (r',t')  are  co'drdinates  associated  with  the  wave  surface  E.  (Here 
/  £  \ 

F'  Ms  a  degenerate  form  of  what  we  call  the  Total  Surface  Spreading 
Function  (TSSF)  F£  '(s'  ,s-s'|e.),  where  we  regard  the  moving 
wave  surface  as  a  time-variable  filter.  The  slow  (temporal)  varia¬ 
tions  of  the  surface  (elevation)  c  vis-a-vis  the  acoustic  signal  permit 


this  parametric"  form  of  the  TSSF  where  now  £  is  explicitly  time-variabl 
as  are  all  factors  which  contain  £  and  its  various  spatial  derivatives. 
The  geometry  of  the  general  (far-field)  source-surface-receiver  configura 
tion  is  shown  in  Fig.  5.1. 
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The  various  elements  of  (5.21)  are  described  below.  We  have 


=  effective  "illuminated"  surface  area,  by  joint  beam  pro¬ 
jection  in  SQ,  cf.  Sec.  6.3,  Fig.  6.1  ff.  (5.22) 

dS  =  dr  =  dx'dy'  =  dxdy  on  S  :  t  =  0;  with  dE  =  dS  /n  ;  (5.22a) 


|c‘  ~  X  +  JL  1  »  where  jr‘  is  always  on  the  wave  surface,  E ;  (5.22b) 

r  =  ixx  +  j^y;  [(x1  ,y* )-system  is  simply  a  translation  of  the 

(x.y)-system]  (5.22c) 

To  =  (RqT  +  r0r)/C0:  time-delay  from  0T  -*■  0$  -  0R:  (5.22d) 

ioT'  -*0R  5  -oT.oR^^oT*  :  u"U  vectors;  (s-22e) 

1T,  i^  =  r/|rt  R| :  unit  vectors,  to  projection  of  scattering 

point  on  E  upon  SQ;  (5.22f) 

rt-  rr.  rot”*  =  IJSr^  *  I-rI  ’  etc-»  distances;  (5.22g) 

t'  =  t  -  At^s  =  doppler-shifted  epoch;  (jr ’ , t ' )  are  on  the 

wave  surface  e;  (5.22h) 

^s=iT*r'(r,t)  =]T.{r  +£(r,t-RR/c0)}  :  RR  =  RoR  “  ioR  *  r**(5-22i ) 

L(n.t)  =  (vector)  elevation  of  the  wave  surface,  vs.  z  =  0, 

the  plane  SQ;  ( 5 . 22 j ) 

a  =  t0/2cq  *  absorption  coefficient  (sec  /meter),  cf.  (5.18); 

nz  -  (1  ♦  cj  +  V  ?y  =  lx’  fy  :  surface  slopes  (5.23a) 

n_ =  (outward)  normal  to  surface  E,  =  (it  +  it  -  i  )n  ; 

*  x  ’’j  y  *** z  z 

(inward  normal  =  -£  )  (5.23b) 

Vqt,  VqR  =  beam-steering  wave  numbers,  cf.  (5.13b),  (5.7a);  (5.23c) 

RTR  =  *t)+  ~z^T"hR^ ’  c^'  ^9-  vector  distance  between 

0j  and  0R.  (5.23d) 

From  Fig.  5.1,  it  is  readily  shown  that  the  unit  vectors  j-j.,  jR  become 
for  these  bi static  configurations 


(5.24, 


j.T U  =  l-Sj I  =  J.x  C0S*T  sineT  +  1y  sin<j>T  sin0T  + -z  C0Se 
I  °o  3 


lR|s  =  VIJ?RI  =  -Uxrt  cos^t  sineT  +  j.y(RT  sin*T  sineT  -LQ)+  izhR}/R, 

(5.241 


where 


Rr  =  (R2  sin2eT  +  L2  -  2R.j.L0si n<}>T  sineT  +  h2)1*,  and 


*(o)T  =  *(o)R 


-  ir/2. 


(5.24( 


[For  joT,  JqR  we  simply  replace  Ry  by  RQy,  4>j  +  <t>0-p ,  etc.  in  (5.24).]  We 
also  find  it  convenient  to  write 


2a  5  JT '  ip 


=  £  (1  +  Rt/Rr)  C0S4,T  sineT 

+  jy  {(1  +  rt/rr)  sin0T  -VRR> 

+  jz(coseT  +  coseR) , 


(5.25 


cf.  the  exponent  in  (5.21b).  [Again,  for  the  reference  vectors  joy, joR, 
_R„t,  etc.,  we  set  RT  -*■  R„T,  <j>T  $0y,  etc.  in  (5.25),  to  get  2aQ  =  1  ^  -  jQ 


~oT’  ■'!  ’  ‘oT’  yT 

cf.  (5.27)  et  seq.  below.] 

A  critically  important  simplification  of  our  canonical  results  (5.20), 
(5.21)  results  from  the  ability  to  employ  "narrow  beams"  iN.B.),  e.g.. 


beams  narrow  enough  so  that  the  mutually  "illuminated"  surface  z  {*  SQ) 


(i.e.,  the  shaded  region  in  Fig.  5.1,  for  example)  is  sufficiently  small 
that  the  spatial  geometry  (~Ry,  RR,  etc.)  for  source  and  receiver  change 


little  over  the  "illuminated"  surface  region  (~  SQ).  Thus,  we  can  replace 


Rt  by  R  T,  eT  by  e  T ,  etc.,  and  most  important,  note  that  the  angle- 

01  I  01  /  n  \ 

dependent  quantities  in  F£  ',  (5.21b),  can  also  be  replaced  by  the  (con¬ 
stant)  reference  quantities  eQy,  Q-j- ,  etc.  This  results  in  a  "factoring" 
of  the  beam  patterns,  in  that  they  now  will  depend  only  on  their  projec¬ 


tions  on  the  surface  about  0^,  and  not  on  the  coordinates  of  the  wave 


surface  away  from  0^.  [For  details,  see  Section  6  ff.] 


59 


ST5 

& 


Accordingly,  we  have  jy  ^-loT»  etc.,  so  that  the  TSSF  (5.21b)  can  now 
be  expressed  in  the  ultimately  much  more  usable  form  (where  the  results 
of  Sec.  6.2  are  included): 


FS  *  n.b.  =  77  x2pXdn  V"iRfo/co  ’  ^oR;V  ^(VVV^oT;V 

f.f.  So  (4,t)  RoTRoRnz 

N'B*  fill 

•  or )(r.'.t,)e  0  0  0  0  ---ooo  (5>26) 


where  £  =  £(r,t')  is  a  slowly-varying  function  of  time,  vis-a-vis  the 
signal  in  X^(t),  cf.  (5.21a)  and  remarks  ff.  Eq.  (5.21b). 

Important  special  geometries  of  operation  are: 

I.  Monostatic  Operation:  R  @  T:  (LQ=  0;  6qR  -*■  e^,  4>qR  =  ioR 

RR  =  RT;(RoR  =  RoT^  hR  =  hT^ 

2«nx  =  2  cos<j,oT  sineoT;  2aQy  =  2  sin^rtT  sinenT;  2an7  *  2  cosen 


oT  oT*  caoz 


with  associated  quantities  of  subsequent  interest  (cf.  (7.28)): 
2(4  +  “oy  +  4)/aoz  =  2/coseoT; 

(n/nz)  •  2^  =  2(cx  cos<f»oT  sineoT  +  cy  sin<OoTsineoT  -  coseoT); 


Specular:  0qT  =  0:  .'-2«ox  «  2aQy  =  0;  2aQZ  =  2. 

Bistatic  Operation  (R  f  T):  (LgT4  0;  Eqs.  (5.24),  (5.25),  Ry  -*•  R 

etc.,  generally) 


III.  Bistatic  Operation  at  Snell  (or  Specular)  Angles:  0 


Lo*  <RoT  *  RoR>  slneoT  insures  <ioR>y  *  <ioT>y 
♦oT  '  *oR  '  '/2  *  ,/Z  1nsures  <i0R>x  *  UoT>x 


coplanar 


eoR  *  eoT  ,nsures  UoR>z  *  -UnT)Z:  Sne"  1n  plane 


fLM 

IVM 


Cl 


Then  we  have 


2  cose. 


2(a2  +  a2  +  a2  )/a  =  2  C0S6 

v  ox  oy  oz"  oz  < 


(5.30) 


(j2/nz)  •  2a0  =  ^y[(l  +  ROT/ROR)sin0oT  “Lo/RoR]  •  2  cos0oT  =  *2  cos0oT’ 

this  last  from  (5.29).  (These  results,  cf.  (5.28),  (5.30),  agree,  as 
expected  with  Eqs.  (6.25)  et  seq.  of  Tolstoy  and  Clay  [9],  who  use  a 
geometry  with  0^,  Fig.  5.1,  as  the  primary  reference  system.] 

5.5  Discussion  and  Critique 

At  this  stage  we  have  carried  the  analysis  for  the  received  scattered 

field,  X(t),  to  the  point  where  the  central  physical  problem  is  to  de- 

(  £  } 

termine  the  analytic  structure  of  the  surface-scatter  kernels,  Q£  ', 
on  the  surface  z.  The  approximations  and  assumptions  governing  (5.26) 
in  (5.21a)  so  far,  are: 


) 


(i).  narrow-band 
signals 

(ii).  Snell 

doppler 


(iii).  far-field 
conditions 


these  permit  us  to  separate  the  signal  compon¬ 
ents  from  the  wave  surface  velocity  the  total 
surface  spreading  function  F^(s'  ,s-s' )  is 
~  S(s-s’),  to  yield  F^(s,t| ..)  in  (5.21a), 
(5.26) 


"'"max^max^o  <<:  RoT’  RoR 


,  where  a  *  max 


(rms)  correlation  distance  of  the  (illuminated) 

wave  surface;  L  =  maximum  dimension  of  the 
max 

transmitter  (receiver)  aperture.  This  permits 
us  to  define  beam  patterns,  which  are  inde¬ 
pendent  of  range. 

(iv).  small  surface 

displacements  vis-a-vis  Rqj,  RqR:  (part  of  (iii),  really) 

(v).  "narrow  beams"  ?  ? 

(N.8. )  AS0  «  RqT,  Rqr'>  beam  patterns  are  independent 

of  angle  variations  and  can  be  referred  to  a 

single  point  (0^)  on  the  "illuminated"  wave 

surface.  [If  the  "narrow-beam"  condition  is 
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.v\ Vi .  - .  v  «-• 


(vii). 


(viii), 


not  satisfied,  (5.21b)  is  to  be  used  in  place 
of  (5.26).] 

vc  =  0:  c(z)  =  cQ,  a  constant:  no  velocity  gradients  in  the 
vol ume ; 

Qy  =  0:  no  volume  inhomogeneities,  random  or  deter¬ 

ministic,  cf.  Eq.  (5.3b); 

Only  local  surface  scatter  interactions:  the  scattered 

field  does  not  couple  to  the  scattering  ele¬ 
ments,  cf.  (5.3c,d)  et  seq. 

M^,  -►  M^:  the  scattering  surface  is  sufficiently  re¬ 

moved  not  to  affect  the  pressure  distribution 
over  the  aperture  of  the  driving  signal 
source,  cf.  (5.1).  (This  condition  is  always 
obeyed  in  practice  unless  the  boundary  is 
very  close  to  the  aperture  and  the  source 
level  is  very  high  [41,  I]). 

A  non-dispersive,  lossy  medium  is  assumed,  so  that  the 

Helmholtz  equation  (5.3a)  for  the  propagating 
field  is  obeyed,  with  only  a  frequency  and 
range-dependent  attenuation  factor  applied 
to  the  emitted  signal  amplitude,  AQ,  cf. 
(5.18),  in  conjunction  with  (i)  above.  For 
practical  oceans  this  requires  f  o(40  kHz). 

Slowly  moving  surfaces  vis-a-vis  c _ and  the  time-scale  of 


the  injected  signal.  This  permits  us  to 
treat  the  time-variability  of  t(r,t')  para¬ 
metrically  in  (5.21)  and  subsequently. 


So  far,  no  boundary  conditions,  and  therefore,  boundary  approximations, 
have  been  explicitly  invoked.  These,  for  the  moment,  are  implicit  in  the 
surface-scatter  kernels,  Q^,  including  such  important  factors  as  re¬ 
flection  coefficients  and  shadowing  functions.  We  shall  consider  the 
specifically  in  Section  7.  Of  the  approximations  and  assumptions 

i  (  n\ 

used  to  obtain  the  basic  (total)  surface  spreading  function  ' ,  (5.26), 
the  most  critical  are  (i),  (ill),  (v),  (x).  The  others  are  normally  well- 
satisfied  physically.  An  important  exception  could  be  (vii),  Qy  *  0:  this 
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implies  that  bubbles,  and  particularly  bubble  layers  near  the  surface 
are  negligible  in  their  effects  on  the  surface  scatter  process.  Certainly, 
in  some  cases,  notably  when  there  is  a  great  deal  of  wave-breaking  ac¬ 
tivity,  bubbles  become  important.  But  in  many  others  they  appear  to  be 
"invisible."  (This  point  is  discussed  concisely  in  Section  1  earlier.) 

In  our  present  treatment  we  shall  exclude  bubbles,  accordingly. 

Finally,  our  present  formulation  so  far  contains  several  generali¬ 
zations  over  earlier  treatments.  These  are: 

(1) .  time-variable  wave-surfaces,  s,  with  doppler  delays,  cf.  ( 5 . 22i ) ; 

(2) .  general  (narrow-band)  signals; 

(3) .  steered  (complex)  beams; 

(4) .  a  canonical  scattering  structure  (kernel)  which  can  be  de¬ 

veloped  to  include  a  variety  of  approximation  procedures, 
including  diffraction  effects,  for  different  classes  of  wave 
surface. 

This  latter  is  considered  explicitly  in  Section  7  ff.  The  specific  time- 
variability  of  the  surface  elevation,  £,  is  critically  important.  It 
quite  naturally  provides  the  expected  "tilting"  effect  in  our  subsequent 
two-scale  theory,  whereby  the  small-scale  surface  is  "modulated"  by 
the  large-scale  or  gravity-wave  component,  without  recourse  to  the  ad 
hoc  mechanisms  used  in  previous  analyses  [Sec.  II  of  [4]  and  Sec.  Ill 
of  [29],  for  example.] 


6 .  Beam_  Pa tterns  _ and  Apertu  res 

The  aperture  structure  Ccf.  Ay,  AR,  (5.11),  (5.6)  above]  and  the 
associated  beam  patterns  play  a  critical  role  in  the  practical  applica¬ 
tion  of  the  theory.  We  must,  therefore,  examine  their  effects  in 
specific  detail.  Accordingly,  for  our  beam  patterns  in  (5.26)  we  shall 
begin  by  observing  that,  in  our  postulated  far-field  "narrow-beam  " 
(N.B.)  operation  C(iii),  (v),  (5.31)],  we  can  write* 


;  <I’i0T>VcoRoT  ;  ‘.SR  s  <i-Io>-ioR'VcoRiR  (6'1) 


♦These  results  are  achieved  by  noting  that  jj  =  (R0y-*r)/|  R0y+r|  * 
loT  +£‘  (l-ioTioT)/RoT’  etc”  in  the  far-field- 


where  j20  is  the  vector  distance  (on  SQ)  of  the  point  of  intersection 
of  from  0^  (for  R0j)-  When  the  main  axes  of  the  T-  R-beams  coin¬ 

cide  at  0j,  then  x0  =  0  (and  R^R  =  RqR),  and  the  beams  are  said  to  "over¬ 
lap  at  maximum  gain."  (It  is  assumed  that  the  beam-maxima  lie  along  iy, 
joRl.)  Here  ,aoT ,  j»qR ,  are  dyadics  (=  second-rank  tensors),  specifically 


J0T  =  -  "  -oT>Vr; 


aoR‘  -  ‘  ioR'^ioR'  ’ 


in  which  Ms  the  unit  dyadic  (i,j  =  1,2,3).  TheJQ-|.,  aQR,  are 
also  called  the  Fresnel  dyadics  for  j^,,  respectively.  The  con¬ 
ditions  for  (6.1)  are 

llLax/RoT<<1;  lir=r-I„l/R0R'“l-  (6- 

To  obtain  the  elements  of  ,  we  use  (5.24)  with  <)>T  ->•  <(>oT,  4>R  ■+  <t>oR 
Rt  RqT,  etc.  therein.  We  call  (6.2)  the  "narrow-beam"  (N.B. )  conditions 
as  we  shall  see  in  Sec.  6.3  ff. 

6.1  Beam  Patterns  and  Projections 


We  select  as  a  convenient  and  reasonable  beam  pattern  structure 
the  so-called  (general)  ga us si an  beam  pattern  (see  (6.5)  below): 


-  iu  *(2”pA*  v+i$T 


where  g  3  beam  gain  (>  0)  and  A  is  the  associated  aperture  dyadic 


A  =  A 


axy  axz 


a  a 
xz  yz 


aw„  =  Aa,  I  det  aM  0;  A  >  0, 
yz  —  “* 

b2 


and  where  is  a  constant  phase  here,  for  the  moment.  Since  v=o(L), 

A  is  0(1  )  and  is  thus  a  measure  of  the  physical  area  of  an  equivalent 
planar  aperture;  (the  actual  aperture  may,  of  course,  be  three-dimen¬ 
sional).  The  "beam"  0 (v)  described  by  (6.3)  is  not  necessarily  N.B. 

When  the  beam  patterns  fall  on  the  reference  surface  SQ  of  the 
wave  surface  z,  cf.  Fig.  5.1,  the  wave  number  arguments  of  (6.3)  become 


\»TtT 


>  .  Ti  * 


those  given  in  (6.1)  for  these  now  far-field  situations  and  N.B.  con¬ 
ditions.  Accordingly,  we  have 

A  /  -f  )  -  a  ff  )e  ^VRoT^-*^oT ‘-^+1$T^fo^  .  ^ 

ViT(^T-^oT»fo)  '  9T(fo>e  ;  / 
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-oT  5  ioT  '  iT  '  -0T  V 
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(6.5) 
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where  jjyAy  =J)T,  aRAR  =  J\R  have  the  generic  form  (6.4). 

We  chall  confine  our  attention  here  to  elliptical  or  circular 
(conical)  beams.  Then,  all  the  off-diagonal  elements  of  a  in  (6.4) 
vanish,  and  a  /  b  for  elliptical  cones  and  a  =  b  for  cones  with  cir¬ 
cular  cross-sections.*  Thus,  we  write 


►(6.6) 


*In  fact,  a  in  (6.4)  represents  an  ellipsoid  in  the  volume,  any  plane 
section  of  which  is  an  ellipsoid.  Setting  the  off-diagonal  elements  of  a 
equal  to  zero  (and  adjusting  a  -*•  a1,  b  ■+■  b'  accordingly)  is  equivalent  to 
a  diagonalizing  or  principal-axis  transformation.  Our  specification  of ^a 
by  (6.7)  simultaneously,  in  each  sperture  coordinate  system  (Oy,  0R)  usu¬ 
ally  requires  that  (Oy,  0S  ,  0$  ,  0R)  form  a  plane  perpendicular  to  the 
(x,y)-plane,  cf.  Fig.  (5. lT,  e.g.,  <!>oT=0,it/2  (•*•  4>R=*/2,ir,  cf.  (5.24a)). 

More  generally,  we  can  always  choose  the  orientation  of  the  coordinate 
system  for  the  beam  at  0y  to  insure  the  vanishing  of  the  off-diagonal 
terms  in  jT,  but  we  cannot,  then,  simultaneously  require  this  for  the  off- 
diagonal  terms  of  aR,  when  the  condition  of  the  perpendicular  plane  is 
not  satisfied:  the  projection  of  the  (elliptical)  receiving  beam  on  SQ 
will  be  a  "tilted"  ellipse  in  the  O^coordlnate  system.  An  exception 
is  the  case  where  one  or  both  apertures  are  point  sources,  cf.  (6.14)  ff. 


Moreover,  for  analytic  convenience  here  and  without  any  fundamental 
reduction  in  overall  generality,  we  shall  henceforth  require  the  receiving 
beam  (unless  it  is  omni-directional,  cf.  (6.14))  to  be  so  oriented  that 
°Sj’  »o  =  a^on9  with  and  ®R  ^orm  a  p^ane  Perpendicular 

to  the  basic  xy-reference  plane  of  Oy,  cf.  Fig.  5.1  (with  <!>oy  =  0  or  *12). 
(This  is  also  consistent  with  the  experimental  configuration  of  the  associ¬ 
ated  experiments  [6].)  Thus,  (6.7)  applies  for  bothdy,QR,  (6.5),  (6.6), 
simultaneously,*  where  now  a  -*•  ay,  aR;  b  -*■  by,  bR,  for  the  appropriate 
(diagonal)  elements  of  jay  R.  We  call  such  beams  "perpendicular  co- 
planar"  (=  J-coplanar)  beams. 

With  omni-directional  beams,  which  are  produced  by  (single)  point- 
sources,  we  have  for  the  source  weighting.  Ay  =  gy(fQ)  6(j;-j^)  at  a  point 
so  that  the  resulting  beam  pattern  is 
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0  .  -  2ifi  vj, 

/  9T«(e-  i0)e27r1~^  dc  -  gTe 

Vy 


(6.8) 


Consequently,  with  v, given  by  firWVV  'iR'-WV'o  t‘f-  <5-26) 
with  r  =  0] ,  we  obtain  at  once  the  following  beam  patterns 


omm 


*  9T«o>e 
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ik0(lR(r)-i0R)*(J4yhT-hR]) 


(6.9) 
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Unlike  the  gaussian  beams  above,  these  omni-directional  beams  are  "broad- 
beams"  (B.B.)  since  (6.2)  does  not  generally  apply:  iy^-j^y.  etc., 
cannot  usually  be  approximated  as  in  (6. 1); (however,  see  (6.14)  ff.  in 
the  case  of  random  surfaces). 


*See  footnote  on  page  64. 
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Since  £  has  no  z-component,  we  can  set  the  z-components  of  AQy, 

A)R*  equa^  t0  zer0  without  changing  the  projected  ellipses  jr  •  Aqj  •  r, 
(l"I0) 'J'qr  *  (JC"JT0)  on  s0-  For  additional  simplification  we  choose  the 
coordinate  system  (0S)  so  that  the  y-axes  of  0T,  0$  coincide,  i.e., 

4>0j  =  it /2  (•*•  4>oR  =  tt).  Using  (5.24)  we  readily  find  for  (6.1a)  that  now 


(4>0t  = 
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with 


RoR  ■«RoTsin9oT-y2  *  hl)H 


(6.10a) 
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Consequently,  we  see  that,  for  (6.7), 
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0  a^cos  0QT+b|cos  0QTsin  6qT 


0  aT(eoT) 
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and  so  that 


*We  note  that  det  aQy  =  0,  but  det  a^  f  0  (unless  R  @  T,  cf.  mono¬ 
static  operation);  also  j0j*^0y  =  0,  etc.;  cf.  (6.1a):  minimizes  the 
scalar  ioT-aoT-JoT  (=0),  etc. 
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P  0). 
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Similarly,  we  get  (r^  =  0): 
1 
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where  now 
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From  (6.11),  (6.12)  we  see  at  once  that  the  T  and  R  beam  patterns  form  (a 
continuum  of)  elliptical  projections  on  the  mean  surface  SQ  (=<t>  =  0),  e.g.. 


I*i0T'X|s  =  x2  +  aT^9oT^yZ  =  constant  °)  »  (6.13a) 

i°o 

I  X.’ioR*j:|s  =  x2  +  aR(eoT  L,hR)y2  =  constant  {>  0);  j;o  =  0.  (6.13b) 
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6.2  Remarks 

A  variety  of  features  of  these  beam  patterns  is  to  be  noted: 

(1) .  The  gaussian  beam  patterns  [(6.3),  (6.5),  (6.6)],  which  have 

no  "side  lobes,"  can  be  well  approximated  in  practice  by  the 
beam  patterns  produced  by  parametric  transmitters  (and  re¬ 
ceivers),  which  likewise  have  ignorable  side  lobes. 

(2) .  The  use  of  these  side-lobeless  patterns  provides  significant 

analytic  and  computational  simplification,  particularly  in 
the  evaluation  of  scatter  intensities,  cf.  Sec.  8  ff.  These 


"lobeless"  patterns  also  enable  us  to  distinguish,  and  avoid 
confusion  with,  the  diffraction-grating  patterns  produced  by 
scattering  from  any  swell  components  in  the  wave  surface  (a 
point  noted  earlier  by  Tolstoy  and  Clay,  p.  197,  [9]). 

(3).  A  far-field  condition  is  also  required  here,  cf.  (6.1),  to 
enable  us  to  replace  the  various  angle-dependent  quantities 
(e.g. » ja)  with  "central"  values  (o^)  defined  at  0,.  (and  0S , ) , 
and  thus  greatly  simplify  the  subsequent  evaluations.  Other¬ 
wise,  the  exact  expressions  for  jR,  etc.  must  be  used,* 
as  in  the  case  of  the  omni-directional  beams  (6.9)  with  bi- 
static  operation*  (R  f  T). 

With  mono-static  operation  (R  @  T),  however,  J^iz(hp-hR)  = 
£0  =  0  (by  choosing  0T  at  ^  =  0)  and  then  GtT|omni  =  9T; 

^R | omni  "  gR* 

Also,  if  either  transmitter  or  receiver  employs  a  "narrow- 
beam"  (N.B.)  while  the  other  is  "broad-beam"  (B.B.)  like 
| omni  in  for  examPle»  the  narrow-beam  is  always  con¬ 

trolling  in  the  product  0^1..  Then  v  =  -^r-vqR  = 

Eq.  (6.1),  (jcq=0)  specifically  in  (6.9),  e.g. ,  J.r(x)-1qR  = 

^•i0R)fo/coRoR’  so  that  now* 

-  n  tf  Jko(^^oR*%)/RoR 
omni  '  gR(fo)e  (| 

RjR  =  Eq.  (5.23d),  RqR  =  Eq.  (6.10a), 

when,  say,  dp  is  described  by  (6.5).  Moreover,  it  is  no 
longer  necessary  that  Op,  0^,  0R  form  a  plane  perpendicular 
to  the  xy-plane,  since  the  projections  of  an  omnidirectional 
beam  on  a  plane  are  always  circles. 


♦Actually,  if  the  maximum  correlation  distance  of  the  (here  random) 
wave  surface  is  small  compared  to  the  minimum  distances  to  the  trans¬ 
mitter  and  receiver,  the  "illuminated"  portion  of  the  surface  is  split 
up  into  essentially  independent  regimes,  each  of  which  is  in  the  far- 
field,  so  that  Eq.  (6.14)  for  omnidirectional  beams  can  be  employed. 


In  effect,  the  far-field  condition  here  requires  that  the 
/areaj  of  joint  "illumination"  by  the  transmitter  (and  "viewing," 
by  the  receiving  beam),  over  which  the  resultant  beam  pattern 
magnitude  is  significant,  be  small  compared  to  the  distances 
RoT’  RoR‘  Quantitatively,  this  may  be  expressed  by 


/area,  =  Ar.  (=  /ad^)  <<  R  T,  R 


(6.15) 


where  a,  d  are  respectively  the  major  and  minor  semi-axes  of 
the  elliptical  projection  on  SQ,  at  some  level,  -Nbdb,  down  from 
the  maximum  joint  beam  gain  9j9r>  cf.  (6.17),  and  III,  Sec.  6.5  ff 
A  reasonable  generalization  of  (6.5),  (6.6)  is  to  replace  the 
constant  phase  terms  $-[■,  fR  by  the  linear  terms  (in  v): 


$T  =  2™- *T(fQ);  *r  =  2»v*iR(f0)» 


(6.16) 


where  now  is  given  by  (6.1)  in  these  N.B.  cases,  and  R 
has  the  dimension  [L]. 

(5).  Me  note  the  dependence  of  the  beam  patterns  on  frequenc 


this  is  a  direct  consequence  of  the  definition  of  beam  pattern 
as  the  spatial  fourier  transform  of  the  (physical)  aperture 
weightings  A^C^),  AR(T],  ),  cf.  (5.6),  (5.12),  (5.13).  [Tolstoy 
and  Clay,  however  ([9],  Sec.  6.3  et  seq.),  do  not  explicitly 
include  this  frequency  dependence,  which  leads  to  a  somewhat 
different  interpretation  of  scattering  intensity.]  Moreover, 
as  (6.3),  (6.5),  (6.6)  show,  the  beam  patterns  behave  as  ex¬ 
pected  vis-a-vis  aperture  size  and  frequency.  Then,  as  AT, 

Ar  cf.  (6.5),  (6.6),  are  made  larger,  the  beam  patterns  become 
sharper  or  "narrower."  Similarly,  as  the  signal  frequency  (~f  ) 
is  increased,  so  also  do  the  beam  patterns  contract:  as  wave¬ 
length  (a  )  is  decreased,  the  (fixed)  physical  aperture  becomes 
(acoustically)  larger.  In  addition,  the  beam  gains  [gT(fQ), 
gR(f0)],  and  phase  parameters  (<j> y ,  4> R) ,  are  also  frequency- 
dependent. 
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6.3  The  "Narrow-Beam"  Case 

In  order  to  simplify  the  evaluation  of  the  Total  Surface  Spreading 
Functions  (TSSF's),  (5.26),  as  noted  above  we  must  use  beams  (or  at  least 
one  beam,  here  Gtj.)  which  are  sufficiently  "narrow"  that  their  effective 
projections  on  the  sea  surface  (specifically  on  SQ),  are  invariant  of 
angle  variations  over  that  projected  area,  cf.  remark  (3)  above.  Sec.  6.2. 
These  projections,  of  course,  depend  on  the  surface  distance  r,  cf. 

Fig.  5.1. 

Here  we  shall  give  the  conditions  relating  aperture  size  (~Ay,  AR), 
frequency  (~kQ),  and  distance  (RqT,  RqR)  which  permit  this  replacement 
of  angle-quantities  by  these  single  values  established  at  0^  on  SQ.  This 
is  what  we  mean  by  "factoring  of  the  beam  patterns."  Accordingly,  we 
begin  with  the  elliptical , -L-coplanar  beams  of  Sec.  6.1,  where  j;0  =  0 
so  that  0,j  ,  coincide.  From  (6.5)  and  (6.6),  with  (6.11),  (6.12), 
we  can  wrile  for  the  combined  gaussian  beam  patterns 


(6.17) 


where  (6.16)  is  used,  and  where 


B  "  MeoT^ 


'oT 


+  aR(eoTl  LcfhR^ 
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(**  A) , 


(6.17a) 


c  =  CioT  \£tVRoT  +  (^oR  ’  Vx/RoR; 

(6.17b) 

D  =  ^oT  **Vy/RoT  +  (aoR  '  ~R^y^RoR 
2  2 

with  a^,  aR  given  by  (6.11a),  (6.12a). 

We  distinguish  two  principal  geometries: 


A. 


Backscatter  Geometry:  R@T  (identical  beams) 

A  =  2At/RqT  ;  B  =  2ay(eoT)AT/R^T;  B  =  ay(e0-|-)A- 


(6.18) 


In  addition,  with  circular  (conical)  beams  =  a^  and  thus  a^(6oj)  = 
ayCOS^eoT,  cf.  (6.11a). 

B.  Forward  Scatter  Geometry:  Specular  Direction  (identical  beams) 

p  p 

In  this  case  (5.29)  applies  and  again  ap(e0jl  •••)  =  a^e^), 

(6.11a),  with 

A  -  At(1/r2t  *  1/R|R);  B  =  a2(eoT)AT(l/R2T  +  1/R*R)  =  a^lA  (6.19) 

for  these  identical  beams.  Here,  however,  RQ^  t  Rqj  usually. 

For  the  somewhat  more  specialized  beam  choices  of  the  Roderick 
experiments  [6],  where  the  receiving  beam  is  omni-directional,  cf.  (6.9) 
and  (6.14),  we  have  (6.17)  for  the  combined  beams,  where  now^  = 

J0r  *  «B-|-r/r0r>  Lj  t  0-  Thus,  the  parameters  (6.17a,b)  of  the  combined 
beams  are  now 


A  »  At/R^t  ;  B  »  ; 

C  =  ^oT'~tVRoT  +  VfWs1n<l’L; 

D  '  (ioT-iT>y/RoT  *  I1 


(6.20) 


For  the -L-coplanar  condition,  4>^  =  0,  consistent  with  <{> qT  =  n/2,  cf. 

Fig.  5.1.  However,  since  the  receiving  beam  is  omnidirectional,  (6.14), 
this  condition  can  be  relaxed,  and  does  not  necessarily  vanish.  For 
the  back-  and  forward-scatter  geometries  of  (6.18),  (6.19),  we  find  that 
(6.20)  reduces  to: 


C.  Backscatter  Geometry  (Gauss  x  Omni):  <t>,  =  0 


A  =  AT/RoT;  B  =  aT^9oT^A;  C  =  *Tx/RoT’ 


D  -  i'Tycos^9oT/RoT  + 


Tx  oT’ 


(6.21) 
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D.  Forward  Scatter  Geometries  (Gauss  x  Omni), 
in  the  Specular  Direction 

Here  we  have  <t>^  =  0,  <j>Qy  =  ir/2,  4>qR  =  n,  and  A,  B,  C,  D  are  readily 
found  to  reduce  to  (6.21),  provided  we  invoke  (5.29),  as  required  for  the 
specular  direction.  Off  the  specular  direction,  of  course,  C  and  D 
differ  from  (6.21). 

Our  next  step  is  to  relate  an  effective  beam  width,  ab,  to:  (i) 
the  beam  parameters  (Ay,  kQ,  ...);  (ii)  the  "narrow-beam"  conditions  (6.2) 
above.  This  is  done  by  choosing  some  distance  (r.  )  on  S„  where  the 

maximum  projection  of  the  beam  falls  to  b_i  of  its  maximum  value  (gy), 
(6.3),  in  the  manner  of  Fig.  6.1. 


IU.T:  v 


b-max 


ABOoy) 


Figure  6.1.  Geometry  of  effective  beam  width,  ABy,  in  the  far  field, 
rb-max^oT<<:  lowing  illumination  by  "short"  pulses. 
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For  the  gaussian  beam  (6.5)  we  get  (in  the  far-field) 


^T^Rb"“oT’f°^  =  9t  6XP  ^oAT^xb+aT^6oT^2yb^RoT}  =  9T/b' 


(6.22) 


Taking  the  largest  projection  on  SQ,  viz.  yb,  (aT<  1)  and  solving 
(6.22)  for  Yb/RoT.  =  0,  we  get  directly 

*b  5  Vmax'RoT  "  12  109  WAl4(9oT>ko)‘S<<  '•  <6-23> 


where  we  recall  the  far-field  condition  (6.2).  From  the  far-field 
geometry  of  Fig.  6.1  it  is  easily  seen  that  the  effective  beam  width  is 
now  (in  the  far-field) 


iBt,b  1  2*bcos8oT  ' 


(6.24) 


subject  to  the  condition 

_  AB,  . 

*b  5  W/RoT  ’  /2  109  b'\  1  koaT<8oT>  1  5^  <<1>  b  *  l-  <6-24a> 

2  2  2  ^ 

Here  (6.11a)  gives  aj(90j)  =  aT  coseoT*(cos  0Qy  +  by/aT)  sin  0Qy)2 

(=  aT  coseQy  for  circular  beams,  ay  =  by). 


E.  Example 

These  data  are  obtained  from  experiment  [6].  We  have 


ASy  b  =  5°  9  -10  db  (:b  =  10)  =  8.73-  10"2  rad. 
0qt  =  83°  =  minimum  grazing  angle  used. 


•  -yb  =  0.358 . 


(6.25a) 


-  _p 

which  is  not  very  small.  We  should  have  y^  <  0(10  )  comfortably. 

Since  hy/cos0Qy  =  RQy,  we  easily  find  that  for  hy  =  28.7  m,  RqT  =  235  m, 
and  rb  =  84.3  m,  from  (6.24a).  (This  is  approximate,  since  yb,  (6.25a), 
does  not  obey  (6.24a)  very  well.)  The  effective  bandwidth  of  5°  at  -10  db 
(from  gy)  is  measured  at  20  kHz.  This  permits  us  to  estimate  the  equi¬ 
valent  aperture  cross  section  Ay,  under  the  assumption  of  fully  circular 
beams  [aT  =  by  (=  1)].  Since 


AT  =  2(aTkoA6T,b>  1o9eb 


(6.25b) 


from  (6.24a),  we  find  that  at  f  =  20  kHz,  or  xQ  =  (3/40)  m,  »*«  Ay  =  6.16 
(b  =  10),  cf.  Fig.  6.1. 

Now  yb  ~  0.36,  cf.  (6.25a),  does  not  very  well  satisfy  the  far-field 
requirements.  To  overcome  this  critical  defect  (which  otherwise  greatly 
increases  the  practical  complexity  of  our  results,  vide  Sec.  6.2,  and 
above),  we  can  achieve  the  desired  effectively  narrow  beams  by  using 
suitably  "short"  signal  pulses,  so  that  only  a  comparatively  small  por¬ 


tion  of  the  projected  beam  area  on  SQ  is  illuminated  at  any  given  in¬ 


stant,  cf.  Fig.  6.1.  If  ts  is  the  pulse  duration,  then 


r.  -*•  Ar.  =  2t  c  /cose  Ar.  /R  T<<  1, 
b-max  b-max  s  o'  oT  b-max'  oT  ’ 


(6.26a) 


as  required,  where  2tscq  is  the  pulse  distance  in  the  medium.  The  new 


effective  beam  width  is 


AB  =  Ar.  cose  -/R  -  =  2t  c  /RrtT. 
s  b-max  o T  oT  s  o  o T 


(6.26b) 


Thus,  for  ts  =  10'  sec,  Arb_max  =  24.6  m,  and  Arb_max/R  T  =  0.10, 


which  is  better.  Still  better  are  pulses  with  t$  =  4  10  secs,  say. 


as  then  the  far-field  condition  becomes  Arb-max^RoT  =  <<:^’  accePta^^y- 

However,  there  is  a  limit  to  how  short  a  signal  duration  can  be 
tolerated  before  our  postulated  "narrow-band"  condition  [(5.31),  (i)] 
breaks  down  and  we  must  then  account  for  the  fact  that  our  apertures  are 
frequency-variable,  i.e.,  they  act  like  (linear)  filters,  cf.  (5.6).  We 
discuss  this  point  in  Sec.  6.4  following. 


6.4  Beam  Convolutions 


To  appreciate  the  filtering  action  of  the  aperture  in  its  response 
to  signals  of  finite  bandwidth,  let  us  consider  (5.9)  and  (5.10)  once 
more,  now  with  driving  signals  of  arbitrary  temporal  structure.  We  con¬ 
sider  for  the  moment  the  unscattered  field  aH>  Instead  of  (5.10)  in 


the  frequency  domain  (  now  with  sQ  ■+■  s),  we  can  write  in  the  time  domain 


the  temporal  convolution 


►  •  »  .H  *>*'•  •  ’  •  »  '  * 


.  '  -.  Vi"* 


/  dT  /  AT(C,t-T+t*(Cj)S.n(T,c)dC;  t*  =  lrU cQ  -  RqT/c 
-oo  VT  ' 

1  (6.27a) 


Alternatively,  this  may  be  expressed,  for  S.  (t,£)  =  S.n(r):  uniform 

1  n  ***  1 1 1 

drive,  as 


a..  - 


^  e  °°  Sin(rfr^(s'^T  ~oT^/co,s^ 


(6.27b) 


where  S^n  (^fj)  is  the  (amplitude)  spectrum  of  the  applied  signal,  cf. 
(5.12),  (5.13).  Equations  (6.27a)  and  (6.27b)  may  be  described,  func¬ 
tionally,  as 


(6.27a) :  ~  aperture  weighting  ©  driving  signal  waveform;  \ 

((6.27c) 

(6.27b) :  <*H  ~F.T.(beam  pattern  x  driving  signal's  spectrum)df,  > 


where  ©denotes  (here  temporal)  convolution,  and  F.T.,  the  fourier  transform. 
For  narrow-band  signals  (about  sQ  =  2"f0)  we  have  Sin(s-g-*-|-Q-~--)  = 

Sin^2rM*  where  Sin  ^  I s'l  *  Af<<fo'  Then  the  9enera1  relation  (6.27b) 
becomes 


‘ott-W'o’  f  ,  ,5'  ,r,  ,  ,,  ..  .  ^‘-Wds 

e  ;Br  s1n<2TT>  <M[iTMoT]<5  +so)/co:5o)e 

1  (6.28) 


2tt 


Provided  S^n  is  spectrally  (considerably)  narrower*  than  Qj»Gl-|-  may  here 
be  considered  essentially  constant  (in  s')  vis-a-vis  Sin,  and  so  (6.28) 
reduces  to 


°H  '(MtlT-ioT]5o/colso>  L  Si„(s/2'i)eS(t  *oT 

Brl 

which  is  just  the  form  of  (5.14)  earlier,  as  expected, 
band"  postulate  of  our  analysis,  generally. 


/co> 


ds 

2  Tri  ’ 


This  is  the 


(6.29) 

"narrow- 


*We  may  set  k  -*•  k',  ( k * < ( < )kQ ) ,  in  (6.5)  and  specifically  compare 
Q.y  with  S^n  to  establish  a  workable  quantification  of  the  term  "narrow" 


Accordingly,  as  long  as  the  pulsed  signals  of  Sec.  6.3,  cf.  (6.26a), 
are  not  too  "short,"  i.e.,  as  long  as  t  <  ofCtp  we  can  "factor" 
the  beam  pattern  and  signal  spectrally  according  to  (6.29).  The  "scannin 
of  the  beam  projection  on  S.  in  Fig.  6.1  is  then  automatically  accounted 


5.26)  in  (5.21a)  through  the  spectral  structure  of  these  narrow 


band  "short"  signals  embodied  in  S 


6.5  Beam  Integrals  for  Scatter  Intensities 


When  the  narrow-band  condition  for  the  signals  is  obeyed,  the  beam 
pattern Q.  and  the  signal  spectrum,  S^n,  "factor,"  as  in  (6.29),  (5.21a), 
(5.26).  When  in  addition  we  have  gaussian  and/or  omni-directional  beam 
patterns  of  the  type  (6.5),  (6.6),  and  particularly  here,  (6.17),  in  the 
far-field,  it  is  possible  to  evaluate  a  variety  of  integrals  over  SQ, 
which  arise  in  the  evaluation  of  the  coherent  and  incoherent  scatter 
intensities  (cf.  Section  8  ff.).  We  summarize  these  results  here. 

The  integrals  in  question  involve  the  beam  pattern  products  &tQr  : 
1TR.  With  the  help  of  (6.1)  in  (6.16),  we  see  that  (as  long  as  one  of 
the  beams  is  "narrow"  in  the  sense  of  (6.2),  now  with  r  =  0) 


*T  4  ko!-  ioT  •  VRoT;  V  korJoR-.VRc 


so  that  from  (6.5),  (6.6),  we  may  write  directly 


W  =  ^TR  ‘  9T9Re 


'  K£-'!tr’£+1i‘o£‘£tr 


where  from  (6.11),  (6.12) 


AT  Ap 

=  -  A  +  — —  A 
D2  -*oT  +  D2  -oR 

RoT  KoR 


at  ar 

R2  R2 

RoT  KoR 


aXT>!!r +  *rKt’ 


W 


Also,  we  have 


-TR  =  A)T  ’  VRoT  +  SoR  ’~R/RoR  =  bTy  =  D  = 


bTx  =  C 


Eq.  (6.17b) 
or 

Eq.  (6.20) 


*Tx/RoT  +  lRx/RoR 

V0^  Win°oTcos9oT  .  .  (RoTSin9oT-y 2  .  .  <RoTs1"9oT-y  Vfc 

- ; -  1 - p -  Vror  - ^5 - 

KoT  oR  oR 


(6.33) 


from  (6.10),  where  we  set  the  z-component  of  b-^  to  zero,  since  generally 
jr  lies  only  in  the  (x,y)-plane. 

The  integrals  in  question  are 


I2(o^,D)  =  /  "Q^r^  /“  d(Ar)0RT(r1+Ar)*e 


Ar  s  r2  -  jTj 


(6.34) 


(6.35) 


I3(4r)  =  /  ORTlljlORjtr^irl'drj 


(6.36) 


(Note  that  since  .rj.Ar  *iave  n0  z-components ,  we  can  set  (a Q)  *  0  through¬ 
out  in  what  follows.) 

Since  Q^,  (6.31),  are  (symmetric)  quadratic  forms  exponentially,  at 
most, over  infinite  intervals,  we  can  use  the  well-known  result,  Eq.  (7.26) 
of  [45]  (in  vector,  dyadic  notation): 
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iC‘U-  iu*A*u  /0  xn/2  -  k‘A-1*S 

1(0  =  /  e*-  2-—dju  =  -^VTe  2-~  “ 

[U]  [det  A P 


(6.37) 


-_1  -.1  ^ 
where  A  -  dyadic  form  of  the  inverse  square  matrix  A  ,  e.g.,  c*u  -+  £u; 

A  mm  —  *'*'  M  «%  %»  M 

£*A*u_=  uAu^,  where  the  elements  of  the  matrix  A  are  the  various  (xx,  xy,  ...) 
components  of  the  dyadic  A.  Thus,  det  £  =  det  A,  etc.  Applying  (6.37) 
successively  to  (6.34)  with  (6.31)  yields  after  considerable  matrix  manipu¬ 
lation  (on  replacing  dyadics  by  matrices,  as  convenient)*: 


.  _  (VR)22.2  _4(V?R)(E-1*iHG-1H)(2«0-krR) 

I  -  . - -  e 


2  /det  GE 


(6.38) 


where 


£  E  J^TR  +  -M  = 


A  +  d. 


12 

B  +  d. 


=  E;  D„  = 

«»n 


dll  d12 


d12  d22 


=  jDn  (6.38a) 


det  Z  =  AB  +  (Bd^  +  Ad22)  +  (d1jd22  '  d12'  >  0; 


(det  E)_i  ; 


(6.38b) 


1 5  £tr  -  *  i  • 


H  R  ML"1  *  i'^TR  *  tt  • 


(6.38c) 

(6.38d) 


and  only  the  (x  and  y)  components  of  2aQ-  £tr  are  used  here. 

Equation  (6.38)  is  an  exact  result.  Note  that  the  exponent  of  (6.38) 


♦Although  JBTR  in  (6.31)  is  diagonal,  DM  in  (6.34)  is  not--hence  the 
utility  of  the  matrix  technique  here. 


r 


ft 


is  independent  of  frequency.*  This  occurs  here  because  of  the  ex 


frequency  dependence  of  the  beam  patternsin  our  formulation  [(6.3)-(6.6) , 

and  (6.31)],  which  is  required,  realistically. 

For  I,,  (6.35),  we  find  that  since  DM  =  0,  •*•  E  =  6x0,  and  G  =  i  BTD, 
1  -1  ->-1  --i  -M  -  r  k  o  £  9TK 

H.  =  21,,  so  that  §.  =  2§.yR  *  2g,  ,  and  since  det  (|-  B^R)  =  (1/2)  det^BTp 

we  get  directly  from  (6.38) 


“a  00  *  ik  ?ct  *Ar 

I1C2b> =  /  ^RT(ii)dri  /  d(^ART(ii+^:) e 


(9t9r) 

1  K  k^AB 


2  2  it2  "  l^^o"^TR^TR^22o“-TR^ 


^(So'-TR^ 


since  det  B,TR  =  AB,  cf.  (6.32b).  Equation  (6.39)  is  also  exact. 
For  I3>  (6.36),  we  get  from  (6.31) 


9T9R)2*  -  |k2  [( ax ) 2 A+  ( Ay ) 2 B]  - i  k0bTR •  at 
k2  M 


We  note  the  explicit  role  of  the  phase  terms  (~ilTR  (6.33))  [cf.  (6.5), 
(6.6);  (6.16),  (6.17)]  in  the  above  results  (6 . 38) - (6 - 40) .  Again,  only 
the  x-  and  y-components  of  2a0-^R,  b^  (6.33),  are  used  here.  Finally, 
because  JE^.r  is  diagonal,  there  are  no  "xy"-components  in  the  exponents 
of  (6.39),  (6.40). 


A.  Approximations 

Let  us  consider  the  second  (square)  matrix  in  the  exponent  of  (6.38),  I 

HC1*.  *  <£'1SrR  +  Jlnti'^irR  -  •  <6-41 


♦This  agrees  with  Tolstoy  and  Clay  [9],  cf.  Eq.  (6.74)  et  seq.  therein. 

-4  .0 

However,  we  note  here  the  kQ  factor  in  (6.38).  rather  than  the  kQ  factor 

of  Eq.  (6.74),  [9].  The  formulation  of  Tolstoy  and  Clay  [9]  does  not  in¬ 
clude  an  explicit  frequency  dependence  (~kQ).  (The  insensitivity  of  the 
exponent  to  frequency  stems  generically  from  its  dimensionless  character.) 


dll’  d22>>  A*B»  cf*  (6-38a)  •**I>>^TR;  and  i~:i 


(6.42a) 


JL  >>^jr’  =  dild22  "  dl2>  ^ 


••  I  »  B. 


-  -TR 


*  *  HGH_” 1  =  °(BTRr1)  «  E'1. 


(6.42b) 


And  we  have 

••  det  GE  =  det  EBTR  =  det  E  det  BTR  =  AB  det  DM  =  AB(d11d22  -  dj2).  (6.42c) 


Consequently,  (6.38)  now  reduces  to 


.  .  (9t9r)22it2 

A2Uo  -TR  1  oT '  .4  An< j  .  .2 

ko  /AB^dlld22"d12 

e"  2^  ^-o~-TR  ^M1  ^  o~— TR^ 

det  _D|^j  -  diid22_dl2  ^  ^  * 

(6.43) 


which  has  "xy"-components  in  the  exponent,  since  is  not  usually 
diagonal,  cf.  (6.39),  (6.40)  above.  The  inequalities  above  are  reasonable 
provided  the  grazing  angle  4>qT  (=Tt/2-e  j)  is  not  too  small,  i.e.,  eoT  is 
not  too  close  to  ir/2.  We  must,  of  course,  test  (6.42a)-(6.42c)  for  actual 
geometries,  as  used  in  the  experiments  of  Roderick  [6]. 


A.l.  Change  of  Dimensions 

2 

Finally,  we  note  that  if  is  absorbed  into  A,  B,  and  kQ  into  C,  D, 
of  (6.17),  (6.17a),  i.e..  A,  B  =  0(L"^),  C,  D  =  0(L_1)  dimensionally,  then 
Ij,  (6.39),  and  I2,  (6.43),  replace  the  kA  factor  by  k^,  while  I^,  (6.40), 


»  J  7  ^  /  y  7  \  v  •  iv  /  »  ,  w  r  —  vw  w  •  i  w  q  ■  w  w  vw  i  Q*  1  w  *3*  \  w  *  ' 

has  no  k^  in  the  denominator.  Furthermore,  in  (6.39)  we  must  then  also 
0  2 

insert  a  factor  k£  in  the  exponent.  In  this  form  the  dependence  of  the 


effective  aperture  size  on  frequency  is  disguised. 


1 


B.  Explicit  Forms:  JbyR,  2ctQ-b>TR 


The  explicit  general  forms  for  b^  are  given  here  by  (6.33),  for  N.B. 
patterns,  in  the  usual  far-field  state.  For  the  more  general  quantity 
22o-  bjR  we  combine  (6.33)  and  (5.25),  to  obtain 


2-V~TR  =-x  +  RoT/RoR^COS<,loTS1neoT  "  * *Tx/RoT  +  £Rx/FW} 

+  iy  id  *  RoT/RoR>s1n*oTsi"eoT  -L</RoR 


['• 


<RoTs1ne0rl J 


\>R 


J  *r/RoR 


’  Rx  oR' 

Arvcos^eoT+2Tzsi 

\ 


neoTcos0oT 


(RoTsin6oT-^hR£l 


\>R 


")i 


(6.44) 


Equation  (6.44)  can  be  simplified  somewhat,  depending  on  the  geometry 
of  the  (N.B.)  arrays  employed.  For  a  vertical  array  we  have 


(i).  Vertical  array: 


*  Vz  ; 


(6.45a) 


(ii).  Horizontal  array:  ^  +  Xyly  ’» 


(6.45b) 


and  for  a  combined  horizontal -vertical  configuration  (to  form  a  beam  of 
effectively  cylindrical  shape),  we  can  superpose  (i)  and  (ii)  to  get 


(iii).  Cylindrical  Beam:  X  ~  iy  +  =  ix£x  +  ±yly  +  Xzlz 


(6.45c) 


=  i  I  +  i  i  ,  (i  =  0: 

— y  X  —7  7*  '  V 


-xx  -z  z’  '  y 
horizontal  array  axis  on  x). 


(6.45d) 


Accordingly,  we  can  set  i  =  0  in  (6.44)  for  "cylindrical"  beams. 

In  addition,  we  note  that  the  maximum  effective  dimensions  of  the 
array  are  small  compared  to  the  fundamental  distances  RQy,  Rqj,  cf.  Fig. 
5.1,  as  required  by  the  far-field  or  Fraunhofer  condition  imposed  through¬ 
out  the  present  analysis.  This  means  that  &x  z/r0t,RoR <<:  * »  so  that 
for  transmitting  beams  directed  in  the  yz-plane  (e.g.,  4>oT  =  n/2),  we  see 
that  (6.44)  can  be  further  simplified  to 


we  recall  that  the  area  of  this  ellipse  (on  <£>  =  0)  is 


=  na  a  , 
d  xy 


so  that  we  have  directly  from  (6.49)  in  (6.50)  here 


■■ 


/ABTeT 


(6.51) 


A 


b 


2 


As  we  shall  see  presently  (Sec.  8),  the  factor  2tt//AB  shows  up  in  the 
calculation  of  scattering  intensities.  From  (6.51),  accordingly,  let  us 
choose  b  =  e,  or  log  b  =  1,  i.e.,  we  select  the  e~*  level  of  pattern 
relative  to  its  maximum: 


A,  =  _ 2jr _ 

1  /ABleT 


(6.52) 


Thus,  is  the  projected  joint  beam  area  on  SQ  at  the  e~*  level.  This 
particular  analytic  result,  of  course,  stems  from  our  original  choice 
of  gaussian  beam  pattern,  cf.  (6.5),  (6.6),  (6.17). 


D.  The  Reference  Surface  Area,  AREp 

The  reference  surface  areas,  AREp,  which  are  employed  in  the  defini¬ 
tion  of  the  scattering  cross  sections  (2.6),  (2.19a,b),  are  arbitrary  and 
may  be  chosen  conveniently  to  simplify  the  result.  In  conventional  prac¬ 
tice  (cf.  [1 ]- [21 ] »  except  [9])  the  area  AREE  on  the  reference  plane  SQ: 
<C>  =  0  jointly  "illuminated"  by  the  projected  beam  patterns,  Q.Rj,  is 
determined  under  the  assumption  thatQRy,  as  projected,  is  uniform  over 


the  reference  area  and  zero  outside  it. 

rZ  /  ^  rZ 


The  associated  area,  AREp,  is 


such  that  AREE»  (>>lc  $)»  where  5.g^c  s  are  respectively  the  correla¬ 
tion  distances  of  the  large-  and  small-scale  wave  surface  components,  cf. 
Sec.  3.  Thus,  since  from  (6.3),  (6.17),  and  Sec.  6.3  above,  we  have 
specifically* 


on  S. 


=  w 


-  £{Ax2+By2)+ik0(Cx+Dy) 


(6.53) 


so  that  the  beam-pattern  integral  1,(0),  which  appears  in  (8.1),  (8.8), 

(01  J 

etc.  for  the  calculation  of  <x>(0),  cf'  (2*6a),  becomes 


♦With  (6.17a,b)  for  A,  ...,  D,  where  we  have  absorbed  the  k^-factor 

-2  0 

into  the  A,  B,  which  in  turn  now  have  the  dimensions  [L  ]. 
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00  2  2 

M°)|  =  /  I^htI  dri  =  ^t9r)  AREF:  uniform  beams 

J  uniform  -® 


(6.54a) 


=  ^g.^)2^:  gauss  beams,  (6.36)  )  (6.54b) 


Here  ^  =  2tt//AB(0) ,  (6.52),  namely,  the  projected  area  on  SQ  of  the 
gaussian  beam  patterns,  (2.11),  (6.17),  etc. 

Similarly,  we  find  that  the  associated  integral  [in  (8.1)]  becomes 
(with  the  help  of  (6.37)) 


V^oVxy)  5  ^  V-E  At<  Vir)*e 


2ik0v‘r  -  5&r-SM'^-r 

★o  i 


d(Ar) 


(gTgR)  aref 


\  o  gy  gy  '  o  u  gx  gy  JV 


(uniform  beams)  (6.55a) 


=  J  VgT9R^ 


2  2* 


bo°gy°gy 


I.  j,  f!i|Z .  ^11 

(  2bo  L  °gx  °gy  J) 


(gaussian  beams).  (6.55b) 


Accordingly,  our  choice  of  reference  area  is 


Aref  =  V2 


(6.56) 


here.  As  we  see  in  Secs.  2. 3-2. 5,  this  brings  our  results  into  agreement 
with  corresponding  portions  of  earlier  results,  viz.,  various  separate 
and  joint  determinations  of  the  large-scale  and  small-scale  scattering 
cross  sections. 


6.6  Summary  Remarks 

In  Section  6  we  have  developed  the  role  of  the  beam  pattern  in  con¬ 
siderable  analytic  detail.  Particular  attention  has  been  given  to  the 
gaussian  pattern  (6.3),  because  it  approximates  well  the  patterns  pro¬ 
duced  by  parametric  transducers.  Included  also  is  the  omni-directional 


pattern,  (6.8),  (6.9),  which  is  used  in  the  receiver  in  the  accompanying 
experiments  [6].  The  product  beam  patterns, Cl^R,  are  always  "narrow-beam" 
from  a  combined  choice  of  geometry  [cf.  Fig.  5.1]  and  injected  signal 
waveform  [cf.  Sec.  6.3]. 

We  note  the  following  in  summary: 


( 1 ) .  Projections: 


(i).  JC0  =  0:  axes  of  the  beams  coincide  at  0^, 


(iii). 


(2).  Narrow-Beams  (N.B. 


cf.  Fig.  5.1. 

(ii).  Elliptical  beams,  which  are  elliptical  on 
S(:<£>  =  0).  Here  0T,  0S,  0R  form  a  plane 
perpendicular  to  the  plane  of  S,  e.g., 

<t>0T  =  0»7l/2 »  and-*<j>oR  =  tt/2,ti. 

(iii).  Since  the  receiving  beam  is  omnidirectional 
in  the  experiments,  we  can  relax  the  planar 
constraint  of  (ii)  above,  cf.  footnote  ff. 
Eq.  (6.6). 

(iv).  The  projected  (composite)  beam  area  on  the 
reference  surface  SQ  (:<£>  =  0)  at  the  e'1 
level  vis-a-vis  the  maximum  (at  x=y=0)  is 
^  =  2*//AB,  cf.  (6.52). 

(2) .  Narrow-Beams  (N.B.):  this  requires  far-field  geometries  and 

possibly  gated  (modulated)  carriers  [Sec.  6.3]. 

(3) .  Explicit  dependence  of  beam  patterns  on  frequency  [Sec.  6.2]. 

(4) .  Explicit  results  for  N.B.  cases  in  the  forward  and  backscatter 

regimes  [Sec.  6.3] . 

( 5 )  .  Concept  of  "effective  bandwidth"  [ Sec .  6.3]. 

(6) .  The  role  of  the  Fourier  transform  of  the  beam  pattern  with  the 

driving  signal  spectrum  [Sec.  6.4],  or  when  beam  and  signal 
waveform  are  factorable,  cf.  (6.29).  (The  gated  signal  should 
not  be  so  short,  i.e.,  its  spectrum  so  broad,  that  (spectral) 
factorization  is  not  possible.  "Beam-scanning"  of  the  (wave) 
surface  is  automatically  accounted  for  here.) 

(7) .  Various  beam  integrals  needed  in  the  evaluation  of  scatter 

intensities  are  carried  out  [Sec.  6.5],  which  include  the 
effects  of  the  phase  terms  in  the  complex  beam  patterns,  cf. 
(6.3). 
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In  general,  the  above  results  are  independent  of  propagation  and  geometry 
provided 

(i).  the  medium  is  not  dispersive,  e.g.,  f  <  o(40  kHz),  only 
absorptive  at  worst; 

(ii).  vc  =  0:  negligible  velocity  gradients.  [However,  the  effects 
of  vc  t  0  can  be  accounted  for  by  suitable  reconfiguration 
into  an  equivalent  geometry--to  be  considered  in  a  later  study.] 

(iii).  far-field  conditions  are  valid,  so  that  one  can  speak  properly 
of  beam  patterns,  i.e.,  Fresnel  corrections  are  negligible. 

With  the  results  of  Section  6  we  are  now  able  to  obtain  explicit 
expressions  for  the  desired  scatter  intensities  [Sec.  8],  once  the  ap¬ 
propriate  local  scattering  model  is  implemented  [Sec.  7  ff .] .  See  also 
Section  5.5  for  a  summary  of  the  basic  structural  conditions. 


7_.  Surfacg  Models:  General.  Scatter  Intensity 

From  (5.26)  in  (5.20),  (5.21)  we  can  now  write  the  received  wave¬ 
forms  X(t)  whose  intensities  we  wish  to  determine  for  the  back-  and 
forward-scatter  cross  sections  which  are  the  principal  aim  of  the  present 
study.  Specifically,  we  have  (the  complex)  waveform 


X(t)  -  l  X(k)(t)  =  l  A  /  S  (s/2*i)F(k)(sU(r,t),..|f)e 
k=0  k=0  0  Br.  in  b  0 


st  ds 


2u  i 

(7.1) 


where  the  Total  Surface  Spreading  Function  (TSSF)  i 


s  now 


:00  = 


Fr7  =  f  0^TQ(k)(i,*t,)e 

So 


■au)ocoTo-  ^to)'2a0+coTo^  dxdy 
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RoTRoRnz 

(7.2) 


subject  to  the  various  approximations  and  assumptions  noted  in  Secs.  5.5, 
6.6  above.  The  beam  pattern  product  is  given  explicitly  by  (6.17), 
(6.31)  here,  cf.  Remarks  in  Sec.  6.2  also,  and  t’  =  t-At.c,  r‘  =  r+c, 
cf.  (5.22a,b). 
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7.1  Wave  Surface  Models 

As  before,  cf.  Sections  2.1,  3.1,  we  shall  employ  a  two-component  wave 
surface  model,  consisting  of  the  single  gravity-capillary  wave  surface  as 
one  component  (6),  on  which  rides  the  sol iton-ripple  contribution  (S), 
as  sketched  in  Fig.  2.1  above.  For  later  analytical  convenience  (cf. 

Sec.  8.6)  we  shall  in  addition  split  the  former  surface  into  a  sum  of  large- 
scale  "gravity"  (G)  and  small-scale  "capillary"  (S)  contributions,  so  that 
the  surface  elevation  at  (r,t)  on  the  reference  surface,  (SQ)  is,  in  detail, 

=  iz^g(r,t)  +  ec(r,t)]  +  DG(r,t)ss(r,t)  ;  with  ^  £$  =  n^s 

(7.3) 

where  is  the  normal  to  the  G  (=  g  +  c)  surface,  viz. 


nr  =  Uj.v^Y  +  i  “  i  )n  1 

"** «  x  a  y  z  / 
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(7.4a) 


However,  for  the  present  we  use  the  two-component  composite  surface 

£(r«t)  =  izCG(r,t)  +  nG(l»t)cs(r,t)  , 

reserving  to  Sec.  8.6  ff.  the  further  dichotomy  of  cG  into  a  "gravity" 
and  a  "capillary"  (c)  component. 


In  our  alternative  approach,  as  noted  in  Sec.  5.4  above,  we  obtain  the 
following  specific  results  for  surface  operator  kernels  ([41],  II): 


(7.4b) 


^1:|QS0)  =  rl0,|/Yo°)<v(0)>RL  ^  t 


(7.5) 
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Here  Hq  is  a  plane-wave  reflection  coefficient  (for  the  water-air  interface) 
and  s  is  a  shadowing  function  (0  <  S  <  l),  whose  statistical  properties  we 
shall  comment  upon  presently,  cf.  Sec.  7.7,  C  ff.  The  subscripts  g,  c  in¬ 
dicate  the  (moving)  surface  upon  which  the  quantities  in  question  are  to  be 
evaluated.  Our  specific  result  (7.5)  is  based  on  the  Tangent  Plane  (T.P.), 
or  Kirchoff  approximation,  cf.  Chapter  7,  Sections  19,  20,  [1],  in  particular 
the  discussions  in  Sec.  19.2.  (See  also,  Bahar  [23].)  In  (7.6)  the 
are  the  fluctuations  in  the  densities  of  the  repsective  k-coupled  scattering 
elements  on  £  and  [A  full  treatment  of  the  cases  k>l  is  described  in 
Middleton  [41,  II,  III]. 

Applying  (7.5)  to  (7.2)  then  gives  us  directly  for  the  TSSF 
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where  all  the  conditions  and  approximations  of  Sec.  5.5  above  are  in  force. 
The  result  (7.7),  with  (7.1),  is  a  generalization  of  Tolstoy  and  Clay([9], 
Eq.  (6.19),  before  integration  by  parts  therein,  cf.  (7.27)  ff.  also),  with 
the  inclusion  of  one  or  more  of  the  following: 

(i).  doppler,  in  c(r,t'); 

(ii).  general  apertures; 

(iii).  general  narrow  band  signals;  ( 

(iv).  shadowing; 

(v).  absorption. 

Similarly,  we  find  from  (7.6)  in  (7.2)  that  the  TSSF  for  the  diffrac¬ 
tion  terms  (k^l)  is 
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where,  as  above,  x,  is  given  by  (7.3) 


7.2  The  Second-Moment  Functions  of  the  Received  Scattered  Field,  X ( t 


The  second-moment  function  of  the  received  scattered  field  is  the  sum 
of  the  (k^O)  second-moment  functions,  viz. 

Hx(t)  =  l  Mjk)(T)  =  l  \  Re  <X(k)(t1)X(k)(t2)>^  ,  t  =  t2-t1  ,  (7.9 


since  all  k  (>o)  components  of  X,  cf.  (7.1),  are  statistically  independent 
and  postulated  here  to  be  stationary,  as  well.  [The  averages  <  >  are  over 
the  epochs  of  the  injected  signal.] 

Since  Sin(s/2"i)  -*  Sin(s/2in  !<t> )  here  in  (7.1),  and  since 

S*+s 

<Sin{sl/27ri^)Sin(s2/27li^)%  =  lSin(sl/2l,i),Z<5(4rii)  ’  (7-1 

with  the  signal  independent  of  the  scatter,  we  readily  find  on  applying  (7.1), 
with  (7.7),  (7.8),  to  (7.9),  that  the  various  second-moment  functions  of 
the  received  scattered  field  are  now  specifically: 

A.  Total  "Classical"  Component  k  =  0: 

iMl-  mJ0)(t)  =  Re  {G(1)k2Ko(T).ne  |  f0,«0U)} ,  (7.1 


where 
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for  narrow 
band  signals; 

(7.1 


and,  for  the  present. 
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(7.11c 


where  at  -  r^-r^,  t  =  tg-tj,  and  here  the  composite  wave  surface  c,  (7.3), 
is  represented  explicitly  by 

-1  =  ^rV^/S)’  ^2  =  ^Il+^’ti+T_RR2/co);  RR/co  =  ^RoR~-ioR*~^co’ 

(7.12 

this  last  from  (5.22i),  including  doppler,  where  also  we  note  again  the  time- 
dependence  of In  (7.11c)  we  have  made  the  further  assumption  that  the 
shadowing  and  reflection  coefficients  are  essentially  independent  of  the 
surface  elevation  here.  As  the  exhaustive  analysis  of  Bass  and  Fuks  shows 
([1],  Sections  22,  23),  this  is  reasonable:  the  shadowing  function  depends 
only  on  the  statistics  of  the  slopes  of  the  surface,  in  the  case  of  the  in¬ 
tensity  calculations  here,  cf.  [l],  5,  p.  297.  This  is  also  approximately  true 
for  mean  amplitudes,  <X>,  as  long  as  the  shadowing  is  not  too  heavy,  [l],  5, 
p.  297  again. 

B.  The  Coherent  Component  (k  =  0): 


Our  result  (7.  II)  includes  the  coherent  component,  <X>,  if  any,  of  the 
received  scattered  field.  [We  exclude  throughout  any  direct  (coherent) 
propagation,  cf.  Fig.  5.1.  This  is  very  easily  added,  if  needed.]  The 
second-moment  function  for  <X>  is  readily  obtained  on  setting  x  -*■  t  +tQ  in 
(7.11),  etc.,  and  letting  tQ  -*■  °°.  The  result  is 
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Equation  (7.13a)  follows  from  (7.11c)  (as  tQ->  «)  because  we  have  made  the 
usual  reasonable  assumption  that  the  random  wave  surface  (as  an  ensemble) 
is  stationary  and  homogeneous  (i.e.,  independent  of  (r^tj)),  at  least  over 
periods  comparable  to  the  signal  duration. 

From  (7.12)  and  (7.13)  it  follows  that  the  total  second-moment  function 
consist,  of  the  sum  of  an  incoherent  and  an  (independent)  coherent  part: 


us=ei=  «i0,<T  >  *  "«>(’). 

where  now,  specifically. 
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Eq.  (7.13a) 


C.  Higher-Order  (Diffraction)  Terms  (k>l 


In  a  similar  way  we  obtain  the  second-moment  functions  for  the  dif 
fraction  components  (Iq>1),  with  the  help  of  (7.8)  in  (7.1)  and  (7.9). 
have  directly 
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and  the  averages  are  carried  out  over  all  random  variables: 

etc.,  Rq,  etc.  [The  averages  over  the  reflection  and  shadowing  coefficients 

(/?Q,s)  are  independent  here  of  the  other  averages,  as  noted  earlier,  cf. 

(7.12)  et  seq.,  above.]  Since  these  diffraction  components  can  contain  no 
coherent  components--these  are  observed  only  in  the  "classical"  (k=0)  case, 
cf.  (7.14)--we  have  written  M^(t)  =  Mx-<X>( T) *  in  (7.15).  Generally, 
these  diffraction  terms  (k>l)  are  comparatively  small,  vis-a-vis  the  classical 
components  (k=0),  except  possibly  for  certain  transmitter-receiver  geometries 
and  directional  wave  spectra  where  the  classical  components  may  vanish.  [See 
the  remarks  in  Sec.  3.] 

7.3  Explicit  Surface  Scatter  Statistics:  Two-Scale  Models 


Before  we  evaluate  '  in  the  above  expressions  for  the  second  moment 
functions,  let  us  introduce  an  approximation  for  ( n/n2 )q+s ,  cf.  (7.7),  (7.11c), 
(7.13a). 

Let  be  the  surface  expansion  operator  of  a  truly  two-scale  surface, 
where  the  second-scale  surface  here  is  n^,  cf.  Fig.  2.1,  which  rides  on 
the  single  gravity-capillary  wave  surface,  jzcG,  cf.  (7.46).  Thus,  we  can 
write 

C2 

£G  ■  (1  +  ‘sV*T  *  T  <V*r> +  •••  >  I7-1 

provided  <|nr*jLp?~>  <  1,  so  that  the  series  is  (stochastically)  convergent 
(at  all  r,t).  Physically,  this  means  that  <jng*v|  >  <<$°  and  <£<.>  «1, 

(m.s.),  which  is  not  difficult  to  achieve  here.  Accordingly,  applying 
(7.17)  to  n/nz|Q+s  gives 

(li/nz>G+S  =  1  (^/nz}G  =  (2./nz}G  +  ^2G*^(VnzG)  +  •••  =  <2/n2>G*  {7  A 
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so  that 


<[W(^6]  >:  =  0(0S°GxoGxx)  =  °(10  *10u*10u’1)  =  0(10  \lO'J)  «  1. 

A.  The  Coherent  Component  (k=0) 


Next,  let  us  consider  (7.13a)  and  note  that  the  Rayleigh  number, 
o 

tko0S^coseoT+coseoR^  is  usually  small  compared  to  unity  for  the  soli  ton- 
waves.  Thus,  using  (7.18)  and  (7.3)  or  (7.4b)  we  can  write 


MRT)(°°I***)  "l{(!V9)6  e  °  0zCG(1+2ik0-VJQG  "  T^o'-CG  )24+’ 


First,  we  note  that  for  Cq  a  gaussian  random  process,  Cq  and  its  (first) 
derivatives  are  statistically  independent,  so  that  exp(2ikQa0£Q)  and 
(n •2aQ)/nz)G  are  likewise  independent.  Second,  to  a  good  approximation  we  may 
treat  the  gravity  and  sol i ton-ripple  waves  as  statistically  independent.* 

Using  (7.20)  and  the  fact  that  (n*2o5o)/nz)6  =  2a0X^Gx+2a()ycGy-2a02,  we  see 
that  (7.19)  reduces  directly  to 
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Fl^2aoko^G  =  ^  >  *  characteristic  function  of  (~Cq); 

and,  with  the  help  of  (7.4),  we  see  at  once  that 
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This  latter  may  be  considerably  simplified  if  we  remember  that  ^x,  obey 
a  symmetrical  pdf  (since  ^  is  likewise  symmetrically  distributed.  Then 
(7.22)  reduces  to 


♦Of  course,  this  is  not  strictly  true,  as  the  configuration  of  the  large 
scale  waves  influences  the  local  wind  action  which  produces  the  ripple 
structures  which  ride  upon  them,  [16]. 
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so  that  (7.21)  becomes 
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generally,  the  latter  when  ^^S^G-coh  <<  as  1b  usually  is.  This  [Eq.  (7.23)] 
is  a  new  result,  for  general  beam  patterns  where  (7.13b)  gives  AR-p,  where 
specifically  now 
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is  the  coherent  "tilt-factor."  This  result  is  to  be  compared  with  (6.43), 

p.  205  and  (6.41a)  of  [9].  Here  we  have  the  additional  effects  of  the 

tilted  ripple  surface  (cJ,  which  somewhat  further  decrease  coherence.  But 
2  2  (01 

since  k0osNG-coh  s  sma^  compared  to  unity,  we  can  usually  set  the  ex¬ 
ponential  equal  to  one  in  (7.23b).  Also,  in  place  of  the  angular  illumina¬ 
tion  function  f(e),  (6.25)  [ 9 ]  we  have  2aQZ  (=-2yt+c) ,  as  a  result  of  the 
direct  evaluation  of  (7.7)  in  (7.19). 

The  beam  pattern  function,  AdT,  insures  a  vanishingly  small  value  of 

(01 

M^j'foo]...)  for  scatter  angles  off  the  specular  direction,  since  from 
(7.13b)  and  (6.39)  specifically, 


Abt  a  (i  5)2  (9j9r)^ 


exp  (- 


2otoy~b1 

B 


(7.25 


o 

(where  we  have  absorbed  the  k^-factor  into  the  A,  B's  (and  kQ  into  b^.R) , 
e.g.,  At,  Ad  are  o(L2),  viz.,  are  effective  aperture  areas  now,  cf.  Ia  of 


Sec.  6.5;  also  note  the  factor  k£  in  the  exponent  now).  We  have  replaced 
(i?  S)j.  by  S  (referred  to  the  reference  plane,  SQ,  of  the  gravity  wave 
component)  with  the  help  of  (7.17),  using  the  fact  that  I<cs^g*^o5>I « 
cf.  (7.18a).  Since  A"*  =  (R^/Ay  +  RqR/Ar)»  cf.  (6.32a,b),  is  large,  only 
when  2aox-bypx,  etc.  is  small  will  ART  not  effectively  vanish.  (Here  we 
use  (6.44)-(6.46)  for  2a0-ijR  in  various  configurations.)  Thus,  for  the 
near-specular  directions  only  does  ARy  differ  from  zero  significantly. 

Furthermore,  for  the  gaussian  gravity  wave  surface  we  have  specifically 
the  familiar  result 
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As  expected,  for  rough  surfaces  (2kQ  Voz>  >>1,FlG^°  and  the  coherent 
component  vanishes. 


B.  The  Incoherent  Component  (k=0) 

Our  next  concern  is  the  incoherent  component  (7.14a),  for  which  we 
need  now  M^P(t|...),  Eq.  (7.11c).  At  this  point  we  can  simplify  (7.11c) 
by  the  following  series  of  modifications: 


1).  We  begin  by  expanding  the  exponential  term  in  (7.11c),  using 
(7.3)  viz.: 


2ik  a  U-_-  c  .) 

0-0'-=2  ~1  _  „  O  OZ'  2G  1G'  r.  .  \ 

e  -  e  [1  +  2ik  (a2Gc2S’-lGclS' 


+  ko ^ 2 “o  -G 1  ^ 2^o *^G2 ^  C2S  Hs  "  ko^2£xj 3JgS^  +  •*'].  (7-27a) 


Using  the  postulated  homogeneity  and  stationarity  of  each  component  of 
the  wave  surface,  we  get 


=  e 


2ik  a 
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(  ^2G"C1G^ 


Ce 


•ko(2*o^G)2^S 


+  2lko^o *^2g2^2S'2g1  +  ^ao’&lG^2^*^ C2SMS  +  •••]  • 


2).  Next,  we  reinsert  the  factors  (n*2a0)/nz)g  back  into  the  integrand 
of  Art  in  (7.11c)  and  then  integrate  each  factor  in  ART  by  parts  (one  over 
rj  and  one  over  at  withj^  held  fixed,  since  Cq,  are  stationary  and  homo¬ 
geneous  here,  in  the  manner  of  p.  198  of  [9],  but  only  for  the  leading  term 
of  the  series  (7.27b),  to  achieve  the  desired  simplification.  When  the 
"illuminated"  portion  of  the  surface  is  large  compared  to  (the  square  of  the 
[acoustic])  wavelength,  e.g.,  >>Xq»  or  k£A^>>  1»  or  more  precisely 

>>1,  cf.  (6.52),  we  can  neglect  the  terms  containing  the  derivatives 
of  the  beam  patterns.  The  result  is,  for  the  first  term  in  (7.27b),  where 
we  use  also  (7.18): 


We  emphasize  that  Eq.  (7.27c),  et  seq.,  applies  basically  in  the  "high- 
frequency"  regimes,  where  (2kQaozaG) »  1.  In  the  alternative  "low- 
frequency"  cases  (2k0a0oG)2<<l,  however,  we  may  approximate  the  normal 
factor  (~n)  by  writing 


since  the  slopes  (cgx»5gy)  are  here  small.  This  is  Eckart's  approximation 
[32].  Then,  integration  by  parts  is  not  required. 


3).  Applying  (7.27c)  to  (7.27b)  in  (7.11c)  allows  us  to  write  for 
the  factors  averaged  over  the  various  surface  components  Cqj,  cGx*  ...  etc:* 

♦This  procedure  is  equivalent  to  a  direct  expansion  of  the  exponential 
containing  the  soli  ton-ripple  component  in  (7.7)  and  then  integrating 
(loT"loR^nz^G+S=G  by  parts  ^or  the  first  term  the  resulting  series, 
where  again  we  require  A$k2  »1,  or  k^  »1,  cf.  [9],  p.  198.  This  is 
then  followed  by  the  indicated  calculation  of  the  second-moment  in  question. 


x 


\  r°oxtoo/aoz  1  p  .  .^Vao'WiiG' 

A,  'L  «„z/2  i  2G 


ko^l*2V^nzl  V-2>2^/nz2  V2^*£iG^22o‘£2g)C2Sc1S+-  •  >g,S 

where  specifically  we  replace  the  (interior)  exponential  term  in  the  series 
(7.27b)  by  unity,  since  here  4k^o;0*  r[z)2ss^<<l,  and  where  now 

F2G  '  F2<-2i“ozko'2Vo>G  =  exp  f-boko0GP  -  “G^-'^’ioR^o^1' 

2  2  < 
°G  *<  ’g>  1  bo  "  cos"oT+cos'’oR' 

is  the  second-order  characteristic  function  of  the  gaussian  gravity-wave 
component.  Here  is  the  normalized  directional  space-time  covariance 
function  of  Cq,  viz. 

°o 

°GPG  "  lV~r,T+~r*'ioR/,C0^  =  ^  ^  C1GC2GW2^1G’C2G’^’T  ^d?lG^2G*  ^ 

where,  of  course,  <c^>  *  <Cq  =  0.  The  quantity  Ar’i^/Cg  represents 
a  doppler  delay  produced  by  the  moving  surface  [34],  (a  delay  which  we  can 
usually  ignore  henceforth). 

4).  Since  the  correlation  distance, &q,  for  the  large-scale  surface, 
dG^G’'--^  s  e_1*  is  large  vis-a-vis  £c>  p(*s,...)  =  e"1,  for  the  small- 
scale  surface,  e.g. ,  »*s,  so  that  c2Q  =  for  all  Ar  -  ^s(«  i^) , 

we  may  set  *  ^1G  in  (7.28).  The  result  is 
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where  in  more  compact  form  we  write 


J  '2G  r'o"G-inc^so,''S’  '~G  'S' 


\  r  .  x  u  / 


with  now 


(7.32) 


the  incoherent  "tilt"  factor,  cf.  (7.24).  With  the  help  of  (7.4)  we 


find  that 


m(0)/  \ 

NGS  (Vinc 


'^°‘ox^x't'^aoycy~^aoz^ 

1  (mK> 


(7.33a) 


r(a  C  +a  ?  -a  ) 

v  ox  x  oy  y  oz' 

k  (l+c2+s2) 
v  x  y; 


-/  =  16  4-inc*  [general  ItfT 

•  g  (bistatic,  etc. )  ] 

(7.33b) 


and  specifically,  cf.  (5.28),  (5.30),  with  symmetrical  pdfs  of  cGx,  cGy 
[cf.  (7.54)  ff . ]: 


N^°^(e  )  I 

GS  of  inc-back 


foT=7,/2 


' 16  |si"4v<^> 


l+?x+y  g 


2  2 

+  6  sin  e_rcos  e, 


4 

+  cos  e 


Nvu'fe  1 
G  '  oT'inc-spec 


16  cos  e 


VT*oR=1t/2 


,  (7.33c) 


(7.33d) 


In  fact,  for  the  general  case  here  (with  <f>Qj  =  V2)  we  get 


^ ( a  ). 

GS  v'-o/inc 


=  a4  / 

.,/2  “A 


»oT**/2 


""fe),  *  Kw) 


e 

(7.33e) 


(which  are  evaluated  in  (7.66)  ff.). 


5).  Applying  (7.31a)  to  (7.11),  with  (7.13)  as  it  becomes  (7.23b) 
allows  us  to  write  specifically  the  second-moment  function  (7.14)  for 
the  incoherent  component  (k  =  0): 


M<°^(r) 


■  Re<  koG(1)K„(')fne’1V[HRT)M"-> 


where  now 


mrt>(t  I- 
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-  (2ctQZ)2e  0  02  G  0  S  G'coh  Art]> 


>(|A  ^ 

0  1  ' 


t  2  ,  2  .  2  \  2 
a  +a  +a  \ 

ox  oy  oz  1 

<  »OZ«  I 


)  F2('2lkoaoz,2lkoaoz)G  +  k0NG-incKS(^r,T,")j~ld(-r) 


(7.34a) 


— 2  — 2  »  .  -2ik  a  -Ar 

ART  =  "o  S  /  ^Te 


.*T  e  °~°  ~  drjd( Ar)  =  RQ  S  I^k^),  (7.34b) 


with  Ij  given  explicitly  by  (6.39)  here  for  the  general  gaussian  beam  patterns 
of  Sec.  6.  Following  Bass  and  Fuks  [1],  Sections  22,  23,  we  make  the  indicated 
approximations  of  removing  the  averages  of  the  reflection  coefficients  and 
shadowing  functions  from  under  the  integral  operations,  with  the  observation 
that  these  quantities  are  slowly  varying  over  the  regions  in  which  is 
noticeably  different  from  zero.  (In  fact,  we  can  set  #2  =  7T^  =  i  here  for  air- 
water  interfaces,  but  note  that  s2  t  ?,  usually,  cf.  pp.  288,  291  vs. 
pp.  308,  315  of  [1].) 

6).  At  this  point  we  take  advantage  of  the  fact  that  F^q  and  de¬ 
pend  only  on  the  distance  difference  Ar,  because  of  homogeneity.  Thus  we 
can  integrate  over  n  in  (7.34a),  using  (6.36),  (6.40).  In  addition,  be- 
cause  of  the  short  correlation  distance  ts  vis-a-vis  the  domain  of  "illumin¬ 
ated"  wave  surface,  i.e.,  the  region  on  SQ,  Fig.  5.1,  where  Q^j  is  non- 
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negligible,  we  may  set  Ar  =  0  in  I3(Ar),  (6.40),  e.g.,  since  A,  B,  (6.32), 
are  small  compared  to  unity.  The  resulting  expression  for  (t | . . . ) , 
(7.34a),  becomes  (in  (these  "high-frequency")  relations: 


The  first  term  of  (7.35)  embodies  the  (non-diffractive)  scattering 
attributable  to  the  large-scale  or  gravity-wave  surface,  where  the  Kirchoff 
conditions  are  assumed  to  hold  (cf.  remarks  after  (7.6)).  Here  the  "high- 
frequency"  condition  C( k0b0aG) >:>1  ^  is  assumed  to  apply  cf.  remarks  in  2). 
above.  Otherwise,  in  the  "low  frequency"  cases  (kQb0Og)2  <<1.  (7.27d)  is 
used  directly,  and  the  factor  ^aox+aoy+0toz^aoz^is  rePlaced  by  (2“oZ)2 
in  the  first  term  of  (7.35). 

The  second  term  of  (7.35)  gives  the  scattering  due  to  the  small-scale, 
or  soliton  structure,  which  rides  independently  on  the  large-scale  surface. 
This  is  generally  nonvanishing,  since  Ng^|nc>  0  as  6qT  w/2,  as  we  shall 
see,  cf.  (7.56),  because  of  the  tilted  surface  Ug).  This  component  we 
shall  call  a  perturbational  conponent,  because  of  its  small  Rayleigh 
number  usually,  even  though  we  do  not  employ  the  standard  perturbation 
theory  to  derive  it,  cf.  [4].  Cln  fact,  if  we  were  dealing  with  capillary 
waves  (c^)  riding  on  the  gravity  wave  surface,  the  slopes  of  these  capillary 
waves  are  steep,  cf.  Figs.  4.10,  4.17  of  Phillips  [27],  so  that  conventional 
perturbation  techniques  are  not  strictly  valid,  even  though  the  elevations 
themselves  are  suitably  small,  cf.  remarks.  Sec.  1.3 


C.  Remarks  on  Spectra  and  Covariance  Functions 

The  last  integral  of  (7.35)  is  recognized  as  a  form  of  wave-number¬ 
time  intensity  spectrum,  here  of  the  small-scale  surface  Different 
forms  of  such  spectra,  depending  on  the  details  of  the  definition,  can 
be  obtained  from  the  fundamental  quantity,  the  covariance  function, 
here. 


Accordingly,  let  us  consider  the  following  forms  and  definitions  (for 

two-dimensional  wave  numbers  k  =  (k  ,k  )): 

x  y 


Wa(k| T )  =  /  K  (Ar,x)e~’~r  d(Ar), 


(7.36 


for  surface  (a).  Here  the  fundamental  definition  is  the  space-time 
covariance  function: 


Ka(£r^)  H  <^a<Il*t1Ka(^r*t1+T>  =  hac2a 


<C^>  =  0,  (7.37 


where,  as  before,  we  have  assumed  that  the  surface  elevation,  is  both 

d 

homogeneous  and  stationary.  The  transform  relation  corresponding  to  (7.36) 
is  easily  found  by  multiplying  both  sides  of  (7.36)  by  (2it)  exp(-ik*Ar), 
integrating  over  jk,  observing  that 


j"  e!i-(ir'-5.r)  _iL. 


(7.38 


The  result  is 


(7.39 


so  that,  for  example,  the  mean-square  surface  elevation  of  is,  from  (7.37), 

a 


=  K  (0,0)  =  /  W  (k|0)  ~^2=  , 

*  3  -  3~  [Zv)Z  3 

•**P  a (0 ,0 )  =  Ka(0,0)/o2  =  1, 


(7.40 


(7.40 


when  our  spectral  defintion  (7.36)  is  used,  as  we  shall  do  so  here  henceforth 
Consequently,  the  last  integral  in  (7.35)  is  seen  to  be 

/  Ks(Ar,T+...)e2lk<^0  ^  d(4r)  =  W^k^x)  =  a|ws(2kQa0|T  ),  (7.41 
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the  wave-number- time  spectrum  of  the  small-scale  surface;  w$  is  the  normalized 
spectral  density,  both  defined  according  to  (7.36)  above.  We  shall  use 
this  result  subsequently,  in  Section  8. 

Next,  let  us  consider  the  Wiener-Khintchine  forms  (Sec.  3.2-2,  [45]) 


J£(f)=  2/  Ka(T)e‘la)tdr  =2/  dr  /  Ka(Ar,T)ei^^r"iwTd(Ar) 


W  =  2nf 


(7.42i 


with  the  inverse  relation 


Ka(T)  ■  j  I  W( f)e1WTdf  =  /  K  (Ar.^e^’-dtAr). 


2  1  "a 


(7 -42t 


Also,  in  terms  of  W  (k|t),  etc.,  we  have  the  transform  pairs 

3  **" 


W  (k,w)  =  /  W  (k!T)e‘iuTdT  =  /  /  K  (Ar,T)e1'-*^  ia)Td(Ar)dT 

“  -OO  -00 

0)  -  2nf , 
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3  ~  3  2tt  (2tt)Z 


so  that 


2  2 

a  a 


Ka(0,0) 


J  J  W  (k^)  ^3  =  f  W  (k|0)  -^-y  , 
-  3  (2jt  y  -  (20Z 


by  (7.40). 

Similarly,  using  v.  =Jc/2tt,  f  =w/2tt  as  the  basic  variables  for  spatial 
frequency  and  frequency,  respectively,  we  can  write 

*a(*lT>"/  Ka(Ar,T)e27Ti^r  d(Ar)  =  W  (Jc|T);  v=Jc/2,,  (7. 

“OO 

•*•((  (Ar,t)  =  /  f/a(v|i)e  2it^-  — r  dy  * 


(7. 


so  that,  again, 

<£>  =  K  (0,0)  -  I  Va(v|0)dv  , 

—  CO 

.•.»a(s|0)  -  wa{k|0)  .  o)  »  /*  «a(k,„)  |m  . 


(7.46a) 

(7.46b) 


Also,  we  have 


f)  -  J  va(^.lT)e  1a)T  dt. 


(7.47a) 


(7.47b) 


(7.47c) 


Thus,  (7.47b,c)  are  W-K  transform  pairs. 
Other  forms  are  also  used: 


I.  Bass  and  Fuks.  [13.  Eg-  (3-21): 


(p->Ar;  t-m-  );  WB+p-+Pa 
WB+F^Wa| B+F 


Ka^»T)  =  °apa(^r*T)E  °a^  /  Wa(£»«)B+Fe1&r'-“1“T 


(7.48a) 


Wa^k,w^B+F  =  — f  /pa(^r.T)e'i^’“+ia,Td(Ar)dT. 


Therefore,  from  (7.43a,b),  we  have 


(2rr )3Wa(k;,u))B+p  =  ^•Wa(k,w);  since  Pa(Ar»T)  =  Pa(-Ar ,-r ) . 


(7.48b) 


(7.49) 


or  W,(k,o.). 


T-o-  *Uk,u). 
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II.  Phillips.  C271;  Sec.  4.1: 

(r+Ar;  Za(Ar,T)p  =  =  Ka(^f»T) 


Z->k;  n-Hj 


x(k.«.)  -  /"  /  Z  (4r,T)p.W1i-S  4S!i«l 


(2ir)‘ 


(7.50a) 

(7.50b) 


-  W  (l^)/(2ff)3o?.  etc. 


X  =  W(k,u))p  =  Wa(k,,w)/(2n)3 

t(k)  =  H(k)p  «  (2w)-2Wa(klO)  ^  •*wa(i»<*>)B+F  *  Wa(k,w)p/a; 
*(«)  -  /  W(k,u>)pdk  =  Wa(ui)/2TT 

III.  McDaniels  and  Gorman  [41,  Eg.  (16): 

^a^^McD+G  27  ^  vs* 

•oo 

wa(-l0)McD+S  '  <2-)*lwa<ilO>M1d* 

IV.  Bahar  [23],  Eg.  (14b): 


(7.51) 


(7.51a) 

(7.51b) 


Wa^~l°^B+B  *  71  Wa^°^Mid 


'•Wa(k|O)B+0  =  *  Wa(k|0)ffcD+G 


(7.51c) 

(7.51d) 


D.  Capillary  Wave-Number  Spectra 

For  the  record  it  is  instructive  to  summarize  the  specific  results 
for  capillary  wave-number  spectra,  which  are  needed  in  subsequent  inves¬ 
tigations. 

Proceeding  from  the  defining  relation  (7.36)  we  have 


ik.ar 


Wvkyl°> 5  J  °>e  --  "<£> 


cap' 


(7.52a) 


.  1 


J  I  v  (f  )  <6(v-K_ (o)/2n)>.  df  ,  V  =  K  /2b,  f.  >  0, 

C  CS  —  —s  a  5  —  —  5  5 


(7.52b) 


where  W  is  the  point  (intensity)-spectrum  of  the  (pure)  capillary  surface. 
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c  ,  assumed  here  to  homogeneous  and  stationary,  and  «  is  the  wave  direction. 
Since  in  polar  coordinates 


5  = 


F0(u»s))«(V«). 


2nf  s , 


(7.53) 


where  Kc  and  f  are  related  by  the  dispersion  law  Kc  *  2tt/x  =  a*  w  , 
with  ac  =  Pwater/-3»  where  pwat.«r  is  the  density  of  water  (=  1  gm/cnr)  and 
J  =  surface  tension  (=  74  dynes-cm),  we  find  that 


6  (v 


i,2n^'2). 


fs  > o . 


(7.54) 


Thus,  (7.52b)  becomes  directly,  with  k  -J  +  k^  ,  <p  =  tan-1(k  /k  ): 

X  jr  V  y  * 

3  /r 

"cap<k-*J°>  ’  rm  ‘'c(fs-Ks/W/2>  <5<V5)>5  •  fs*  °-  (7-55) 

4ac 


2 

A  similar  calculation  for  (pure)  gravity  waves,  where  K$  =  u^/g  is 
the  dispersion  law,  gives 


yyo)  ■  I  9  K‘3/2  »g(£  as9)  <s(*v-i)>.  .  f5>  0. 

In  the  above  we  have 


<S((f>y-a)>-  =  w^(ci-a0);  -n/2  <  a-aQ  <  tt/ 2 

=  0,  elsewhere 


(7.56) 


(7.57) 


where  Wj(o-aQ)  is  the  pdf  of  wavefront  directions.  No  backward  waves  are 
permitted  here:  (fg  >  0)  and  (a-aQ)  <  n/2. 

For  our  results  above  K$  =  2n\>  =  2k0sineoT,  <j>oT  =  n/2,  e.g.,  kx  =  0, 
ky  =  2k0sineoT  and  K$  =  (kx+ky)\  with  =  n/2.  Thus,  if  aQ  =  0,  so 
that  the  transmitting  and  receiving  beams  are  cross-wind  to  the  mean  wave 
direction  («  =  0),  Wg(k,$jO)  =  0,  since  w^ ( tt/2)  =  0  here. 

A  variety  of  point  spectra,  w,  is  available.  If  we  choose  the  Pearson- 
Moskowitz  spectrum  [8],  we  can  write  specifically 


St:* 
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yM.(fs)  =  T  e  i  a  ■  8.10  •  10-V;  b  =  0.74g4/v4;  {J  g 

“s  v  =  m/sec,  g  =  m/sec^. 

For  the  purely  capillary  cases  we  have,  accordingly,  us  =  (Ks/ac)^, 
while  for  the  (purely)  gravity  wave  situations,  =  /Ksg. 

In  the  general  case  the  dispersion  relation  of  the  combined  gravity- 
capillary  wave  is 

“s=  (9Ks  °r  fs  =  ST  (3KS  *  p-  KI>S*  (7-5 

so  that 

«(v  -  FL))  =  6(f  .  JL  [gK  +  ^  at2),  (7.61 

s  s  2  /a  s  s  s  c 

in  (7.53)  et  seq. 

7.4  Extensions 

Here  we  summarize  various  additional  results,  needed  in  our  analysis 
above.  These  are: 


A.  The  Diffraction  Terms  (k^  1) 

These  are  obtained  from  (7.16)  in  (7.15),  and  clearly  depend  critically 
on  the  statistics  of  <Av To  date  we  have  only  evaluated  the 
case  k=l :  diffuse  (single-point)  diffraction.  For  this  we  have  found  that 
[41  II,  III] 

jjSzU  -  IJjl^Af-O),  where  |J|  =  |cxxcyy  -  cxy|-  (7.61 

Thus,  | J |  is  a  measure  of  the  surface  curvature  (at  stationary  phase 
points).  Here 


lJlG+S  *  lJG 


1 


1 


°  /l+4+sGy 


Uc!  ;  UJ, 


since  <(l+CGx+CQy)""2>  =  1  +  o(2 agx)  =1;  (see  below  the  coefficient  of 


(7.6 
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cos4eoT  in  (7.66);  also,  |Jg|<<|Jsl,  because  cQx»  are  essentially  con¬ 
stant  in  the  region  in  question  about  where  c<.  has  (local)  stationary 
phase  points.  [Also,  we  can  show  formally  that  [41,  II]] 


(K>1)  <Av[k)Av2k)>R  =  R12)(^rt1^2’t2)  (>0)’  (7‘62) 

(kl 

where  is  proportional  to  the  joint  probability  of  having  k-coupled 
scatterers  jointly  at  (r^tj)  and  (j^tg). 

We  remark,  however,  that  in  most  cases  these  diffraction  terms  (k^l) 
are  considerably  smaller  than  the  (k=0)  components,  cf.  remarks  in  Sec. 

3.3.  In  any  case,  see  Sec.  8.5  for  an  evaluation  of  Mx-<X>^’  (7.15),  (7.16), 

B.  The  Evaluation  of  the  "Tilt-Factor,"  N^j_. 


The  "tilt-factor"  Ng_jnc,  (7.33),  can  be  readily  evaluated  in  the 
isotropic  cases;  (we  reserve  the  general  anisotropic  case  to  a  subsequent 
study).  Before  we  proceed  to  an  evaluation,  however,  it  is  immediately 
evident  from  (7.33c)  that  N^jnc  >0  when  0oT  -*■  tt/2 ,  or  <f>grazing  =  0. 

This  nonvanishing  result  is  qualitatively  consistent  with  other  results 
[cf.  Eq.  (27),  Eq.  A. 14  of  [4],  and  (4)  in  (2)  of  [3];  also  Sec.  Ill  of 
[29]],  as  0oT  -*■  tt/2.  However,  these  earlier  results  use  an  ad  hoc  mech¬ 
anism:  special  modulation  of  the  grazing  angle— to  account  for  the  "tilt" 
produced  here  naturally  through  the  explicit  (slow)  time-variability  of 
nG,  cf.  (7.4)  and  the  comments  after  Eq.  (5.21). 

Since  Cq  is  gaussian,  <cQ>  =  0,  the  pdf  of  the  slopes  is  likewise, 
such  that 


wl^cGx’c6y^ 


2  ,,  2  2.-2 
1“eGx/2oGx'cGy/ZoGiy 

27TaGx°Gy  1 


E  ?Gx?Gy  =  °- 


(7.63) 


We  use  (7.63)  in  the  isotropic  case  oQx  =  o ^  for  direct  evaluation  of 
the  averages  in  (7.33c)  for  backscatter  regimes.  We  consider  the  integrals 

_  o  a  _ /9 _2 


! (0) , (2) , (4)  .  j”j  ll,y2,y4)  0*(x  +y  )/2ox=y  dxdy 


i+x  +y  2* 

2tt  -r2'2i 

=  J  dr  J  d4>[l,r2sin24>,r  si n4 4>3  - - = -  rdrd$ 

0  0  2ttox 


These  integrals  may  be  expressed  in  termsof  the  exponential  integral,  e.g.. 


,  «  -ay  a  „ 

2  Jq  T+7"  dy  =  -  %-  Et(-a  ) 

2  2 

Using  the  fact  that  a  =  1/2 o£>>  1,  here,  and  the  expansion 

-exEi  (-x)  *  2o^  l  (-l)m  m!(2 ;  x  =  a2  =  l/2o^, 
m=0 


we  find  finally  that  (7.33c)  reduces  to 


NGS^0oT^inc-back  =  16{  3(aGx^S1n^6oT+6cGxSin^eoTcos<!0oT+cos^eoT} ; 


2  \2  ._4„  .,2  _,_2„ 


(♦oT^/2)  (or-v=°ou)  * 


Gx  Gy' 


which  becomes  48ogx  (>0)  when  3^  -*•  ir/2,  demonstrating  its  nonvanishing 
nature,  as  noted  above.  Similarly,  from  (7.33b)  and  the  above  we  get 
directly 


NGS  ^oT^inc-for-spec  <t> 0j=d> 0r=tt /2  +  °^2oGx=y^C0S  eoT 


Similarly,  for  the  general  case  we  get  from  (7.33e)  in  these  isotropic 
situations: 


nGS  '  o'inc 


"oT^/2 
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Here  (7.66b),  unlike  the  backscatter  cases,  vanishes  as  e  -  -►  ^/2.  Note, 

(01  0 

however,  from  (5.24),  (5.25)  that  N}n<(  >  0,  generally,  as  eoT  -►  u/2,  cf. 
(7.66c). 


(7.65a) 


(7.65b) 


(7.66a) 


(7.66b) 


(7.66c) 


_C.  Shadowing  Functions  St  £  _  _ 

The  role  of  the  shadowing  functions  S,  £  is  extensively  discussed 
in  Sections  22,  23  of  [1].  We  summarize  the  principal  results  needed  here. 
The  precise  forms  depend  on  whether  monostatic  (ROT)  or  bistatic  (R^T) 
operation  is  involved. 

With  a  =  (tane)/ox=y  =  cot90x/°x=y*  we  have  for  ^  in  the  coherent 
component  [(k=0):  (7.13),  (7.25)] 


^tlweak  shadowing  '  t1  '  ,2  Eq'  <22'39>'  CO-  P-  288, 

a  >>x  (7.67a) 
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(7.67b) 
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| strong  shadowing  ’  iEt<-  <22-52-  tO.  P-  291}/e  .  (7.68) 


Here  we  have  specifically 


x  2 

AN(a)  =  j^{/f  e‘a<f/2  -  a ( 1  -  e[a//2])j;  g(x)  =  ~  j  e_t  dt: 

*  Jit  fl 


✓it  0 

error  function. 


(7.68a) 


These  results  also  apply  in  the  Snell  direction  for  foreward  scattering, 

R^T,  obeying  (5.29). 

2 

For  S  ,  associated  with  the  intensity  of  the  incoherent  components, 
we  find  that 


back- 

scatter 


ROT  '  1+a. 


,  Eq.  (23.29a),  p.  308,  [1] 


*  AN(a)"1  :  strong  shadowing; 

a<<  1 

=<  1  -  Aw(a)|  :  weak  shadowing,  cf.  pp.  308-311,1.1]. 
n  a»l 


(7.69a) 


(7.69b) 


foreward 

scatter 


Nv  1  Nv 


,  Eq.  (23.40),  p.  315,  Cl  ]: 


(7.69c) 


For  details,  see  subsections  4,5  of  Sec.  22,  and  pp.  304-315,  of  [1]. 


The  two  limiting  cases  which  we  consider  here  are  respectively  (1)  the 

"high-frequency",  large-scale  gravity  wave  surface  (G),  where  (k0b0oG  * 

ko^cos0oT+cos9oR^<cG>^2>>1,  and  ^  the  ,,l°w"freclljency"  condition,  where 

(k0b0OG)^<<l>  viz.,  large  and  small  Rayleigh  numbers  for  this  large-scale 

surface  wave  component.  The  small-scale,  or  sol i ton-ripple  surface  (S), 

2 

is  such  that  its  Rayleigh  number  is  always  small,  e.g.  (k  b  <<  1, 

p  p  u  u  o 

=  <c|>.  Moreover,  critical  to  the  specific  results  here  is  the  fact  that 
the  correlation  distance*  of  the  small-scale  surface  is  likewise  small  vis- 
a-vis  that  of  the  large-scale  component,  e.g.,  <<  *q.  (This  fact  ensures 

the  explicit  separability  of  the  two  components  in  the  detailed  analysis, 
cf.  Sec.  7.) 

Our  task  here  is  to  provide  specific  relations  for  the  mean  scatter 
intensities  under  a  variety  of  operating  regimes:  (1)  "high-frequency"  (G), 
backward  and  "foreward"  scatter,  both  incoherent  and  coherent;  and  (2)  the 
same  for  the  "low-frequency"  (G)  cases,  including  specular  (Snell)  and  non- 
specular  directions.  These  results,  in  turn,  are  employed  in  Section  2  to 
give  the  desired  scattering  cross-sections,  which  are  the  ultimate  analytical 
goals  of  this  initial  study.  For  the  most  part,  the  diffraction  terms  (k^l) 
are  ignorable  vis-a-vis  the  "classical"  (k=0)  contributions  studied  here. 

Before  examining  the  various  special  cases  above,  let  us  note  the 
following  general  results: 

8.1  Scattering  Intensities:  General  Forms  for  "High  Frequencies" 

From  (7.34)  and  (7.35)  we  have  directly  the  following  "high-frequency" 
form  for  the  incoherent  scatter  intensities  (t=0)  of  these  two-component 
wave  surfaces: 


1  G(l>Ko(0)1rlff2  s2  [k 


g(i)k0(oj,^  s2 

*  kSNG?Lc<%o>I3(0)WS<^okol0)]  • 


(8.1) 


♦For  the  moment  we  use  an  "isotropic"  distance  i  - 

a  y 

to  describe  this  quantity;  generally,  a  wave  surface  has  two  correlation 
distances,  *x,  % .  ,  cf.  [1],  Sec.  3. 
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•  .  •% 

k.  •  . 


where  the  term  M<x>("l--*)  ^  0  in  (7.34)  for  these  high-frequency  situations, 
cf.  A,  Sec.  7.3,  and  A ^  (=2-n/>/)\B)  is  such  that  k^4j»l,  cf.  (7.27c)  et  seq. 
and  throughout,  2aQ  =  (2a0X>2a0y) »  cf-  Ar*2«0  in  *2’  I3:  ~  has  on1y  x,y 
components. 

Since  F^q  (7.29)  in  (7.34a)  becomes  here  for  the  gaussian  gravity 
wave  surface,  on  expanding  pg(a£,0): 

„  -(ko/2)ir-£H-4r  ,2  , 

F2G  =  e  00  .  e  .  (k060«G)  »1. 

(8.2) 

where  now  the  dyadic 
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=  0;  b  = 


cos0oT+coseoR; 

(8.2a 
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;  det  Dm  =  b^ ( aQx°gy  -  Pxy)  =  bo°GxaGy  (> 


(8.2b 


see  that  I2  in  (8.1)  is  specifically 

2iknct  *Ar  -  (k2/2)Ar*DM.Ar 

WW /%»»r  •  °  ■"-*  d(Ar)  (8.3a 


;  ^(9T9r)2  27t2  ‘  |(22o_^tr) ‘Sm1  * (2“o"^tr) 

k2  /AB(9)detO^ 


(8.3b 


(with  kp  absorbed  into  A,  B,  cf.  remarks  (6.43)  ff.),  from  (6.43)  (subject 
to  the  approximations  (6.42)).  Moreover,  from  (6.40)  and  (6.52),  we  have 


I3(0)  =  (gTgR)2  it/ /AB  =  (gTgR)2A1/2. 


(8.4) 


Applying  (8.2b),  (8.3b),  (8.4)  to  (8.1)  then  yields  the  following  general 
high-frequency  form  for  the  incoherent  scatter  intensity  for  arbitrary 


directions  of  incidence  and  observation,  with  the  finite  dimensions  of  the 
transmitting  and  receiving  arrays  taken  into  account: 


”^(0)  “  G  K«(0)'- 52  ^(9^’2  , 
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/a2  +a2  +c2  V  * 
(  ox  oy  oz  I 

\  V2  / 


a_  o„  b 
Gx  Gy  o 
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Generally,  bTDv  =  0,  if  one  is  sufficiently  in  the  far-field  of  the 

I  KX  yj 

scattering  surface,  so  that  here,  cf.  Fig.  (5.1),  (8.5)  reduces  to  the 
simpler  result 


The  first  term  of  (8.6)  represents  the  "geometrical  acoustic"  (i.e.,  high- 
frequency,  k  -►  °°)  solution,  for  the  specular  point,  or  facet,  (acoustic) 
scatter  from  the  large-scale  gravity  wave  component.  As  expected,  it  is 
independent  of  frequency.  The  second  term  of  (8.6)  represents  the  "per- 
turbational"  or  small-scale  component  solution  through  first-order  Bragg 
scatter  (~Wj.),  which  is  always  nonvanishing  (kQ<  “).  [Various,  less  general 
forms  of  (8.6)  have  been  obtained  by  a  number  of  authors  recently;  this 
point  is  discussed  more  fully  in  Sections  2  and  4.] 

The  coherent  component  of  the  scatter  intensity,  M^(0),  is  similarly 
obtained  from  (7 . 23 ) - ( 7 . 26 )  in  (7.13).  The  result  is 
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which  is  valid  for  all  frequencies  (and  the  various  other  conditions  of 
Sec.  5.5  ,  Sec.  6).  Here  we  have,  cf.  (6.32), 
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cf.  (6.11a),  where  we  have  absorbed  the  k^-factor  in  (6.17)  into  A.p  AR, 
so  that  these  quantities  represent  the  effective  aperture  (or  array)  areas 
of  the  transmitter  (T)  and  receiver  (R).  For  the  high  frequency  cases 
considered  here,  b2agk2  >>1,  so  that  M^(0)  -*■  0,  as  noted  above,  cf. 
(8.1)  et  seq. 

8.2  Scattering  Intensities:  General  Forms  for  "Low  Frequencies" 


For  the  "low-frequency"  cases  we  use  (7.27d),  with  (7.23b),  (7.25), 
to  write  for  the  incoherent  scatter  intensity  (7.34)  now,  on  expanding 

F2G  in  (7.34a),  with  ^ox+aoy+aoz/aoz/2^2  rePlaced  bV  (-2aQZ)2  since  the 
condition  k2/^  >>1  may  not  hold: 
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where  we  have  integrated  over  £  to  get  I.j(Ar),  (6.36),  (6.40)  in  the 
second  term.  Since  A,  B,  (8.7a)  are  small  compared  to  the  regions  at 


where  Kg(Ar,0)  is  significant,  we  can  set  at  +  0  i 
noticeably  changing  the  result.  Accordingly,  (8.8 
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The  coherent  backscatter  scatter  intensity  ■®(0)  (8.7)  vanishes,  because 

2  2  2  2 

of  the  high-frequency  condition,  4cos eoTboko  G  >>^:  tbe  number  coherent 
"specular  points,"  or  facets,  backscattering  at  0Qy,  and  even  as  0Qy  +0 
(Snell  angle),  vanishes  because  of  the  extreme  roughness  of  the  surface. 

The  contributions  from  the  incoherent  "specular  points"  (first-term  of 
(8.10))  likewise  rapidly  vanish  as  0qT  -*■  ?r/2  (grazing),  as  their  number, 
effective  at  0qT,  becomes  vanishingly  small  as  0qT  ■*  */2.  This  is  not  the 
case,  of  course,  for  the  "tilted"  contribution, of  the  small-scale  soliton 
waves,  since  f4-inc>0,  cf‘  (7-56a*b)-  Here,  since  aox=0,  aoy=sin0oy, 
Ws(2aokol0)  “  ws(O.2kQsin0oj|O).  Particular  models  for  Ws,  (7.41),  are 
discussed  in  Section  3  preceding. 

Similarly,  in  the  "low-frequency"  cases  we  use  (8.9)-(8.9c)  (where 
adc=0),  to  get  the  incoherent  backscatter  intensity 
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which  reveals  the  expected  first-order  Bragg  scatter  terms  (~Wg,  W^). 

The  corresponding  coherent  backscatter  intensity  is  found  directly 
from  (8.7),  with  bTR  =  0,  viz.:  I  n  I 


M«><0)  V./2  ;  K“sV" 

Low-Freq  *■ 

back  B  <<:  1 , 


3kos1"2|SoT 

B 


less  eoT  =  0  (vertical  direction),  whereupon  (8.12)  reduces  at  once  to 

*  2G<1)Ko(0)1n^  ^(9T9R)2^k2,  («oT  •  0). 
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8.4  (Bistatic)  Scatter  Intensities 

With  bistatic,  or  "forward11  scatter,  the  facet  terms  vanish  at  high 
frequencies,  except  for  a  small  case  about  the  specular,  or  Snell  angle, 
where  2a„w=2a/%ws0,  cf.  (5.29a),  (5.30),  since  oqx,  ogy  «1,  cf.  (8.6), 
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where  these  facet  terms  become 
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cf.  Bass  and  Fuks,  ([1],  Sec.  20,  1,  2,  esp.  Eqs.  (20.28),  (20.32);  also 
[5],  Eq.  (13a),  where  y=0  at  the  Snell  angle).  The  complete  high-frequency 
incoherent  intensity,  (8.6),  with  forward  scatter  at  the  Snell  angle,  so 
that  (5.29a),  (5.30)  apply,  becomes  specifically  for  these  two-scale 
surface  models 


Mx.<x>(°)  <t>oT=7r/2  =  G  V^in^o  5  (9tV  A 1  [  aGx°Gy 
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(8.14) 


where  by  (7.56),  NQ?]nc  ^  16cos®2.j.  and  8qT  <V2.  For  eoT  -*•  ir/2 ,  strong 
shadowing  becomes  dominant  for  the  facet  term,  so  that  its  contribution 
vanishes.  The  small-scale  term  likewise  disappears  (See  [1],  Sec  20,  (2) 
for  conditions).  From  (8.7),  it  is  seen  that  the  coherent  term  always 
vanishes  at  these  high  frequencies. 

With  bi static  ("forward")  scatter  at  low  frequencies  (8.7),  (8.9)  apply. 
At  the  Snell  angle,  we  again  have  2aQX  *  2aQy  =  0,  bQ  =  2coseoj,  so  that 
specifically  (£„<.  =  0) 
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As  we  see  from  (7.36), 


W  (0,0|0)  =  1  K  ( Ar ,0)d(  Ajf )  =  n tj;a?  , 


where  n  =  0(10°),  and  *  is  an  (isotropic)  correlation  distance.  (See 
Section  3  for  a  discussion.) 

8.5  Scatter  Intensity  of  the  Diffuse  Diffraction  Term  (k=l) 


Applying  (7.52)  and  (7.52a)  to  each  term  of  (7.16)  in  (7.15),  where  we 
(kl  (kl  fkl  (kl 

note  that  <Avq  'av£  =  <Avg  >r<Av5  >j^  *  0,  since  the  gravity  wave  and 
soli  ton  components  are  postulated  to  be  independent,  we  get  directly, 
with  the  help  of  (7.52a): 
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O 

where  J s|  =Uxx£yy"5XylG  S  =  (absolute)  curvature  of  each  surface  (G,S), 

cf.  (7.52).  Removing  the  ~S ^  from  under  the  integral  sign  as  before,  and 
using  (7.3),  (7.4)  we  get  directly  for  (8.16a) 
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from  (6.36),  (6.40),  (6.52):  the  fact  that  13(0)  *  (9y9RrA j/2,  and  the  fact 
that  °QX»aQy  <<1  along  with  the  statistical  independence  of  cG,  and  (Cg^Cgy) 
for  gaussian  processes  (at  the  same  point  in  space-time). 

Accordingly,  from  (7.15)  the  desired  diffuse  scatter  intensity  (k=l) 


becomes 


Mx-<x><°)  *  G<1)Ko(0)in(gTgR)2  ?  tils2  «t|UGl>  +<t|Usl»]  (8. 


given. as  before  by  (7.11a).  As  expected,  these  diffuse  diffraction 


terms  (8.17)  are  (1),  independent  of  frequency,  and  (2)  independent  of 
directionality:  there  is  no  focusing  of  the  scattered  energy  or  dependence 
on  the  angles  of  incidence  and  viewing,  as  in  the  “classical"  case  (k*0). 

The  former  property,  (1),  stems  from  the  "quasi -optical"  character 
of  the  scattering  regions  around  each  stationary  phase  point  on  each  sur¬ 
face,  while  (2)  is  likewise  a  consequence  of  the  fact  that  the  local  sur¬ 
face  has  all  possible  slopes,  considered  over  the  ensemble  of  possible 
stationary  phase  points,  so  that  there  is  no  dominant  subset  of  ^specular 
points"  at  any  given  viewing  angle  (i.e.,  surface  slope).  This  is  quite 
different  from  the  high-frequency,  specular-point  or  "speckle"  situation 
in  the  "classical"  case  (k*0),  cf.  (8.6),  (8.10),  (8.14),  where  only  the 
stationary  phase-points  with  comparatively  large  local  flat  (i.e.,  tangent- 
plane)  areas  (facets)  are  considered,  sufficiently  widely  separated  that 
their  reflected  signals  undergo  large  phase  changes,  and  such  that  kQp  >>1, 
where  p(*  /[r^ |  )  is  the  effective  radius  of  curvature  (rj,  r ^  being 
the  associated  gaussian  radii ‘of  curvature).  Then  it  is  possible  to  dis¬ 
tinguish  numbers  of  such  facets  at  any  specified  angle  (apart  from  very 
small  grazing  angles  where  shadowing  conceals  the  facets),  from  which  a 
resultant,  directional  incoherent  radiation  is  observed,. [5].  In  our  pres¬ 
ent  "quasi -optical :  cases  (k>l)  we  have  kQp<  1,  rather  than  kQp  >(>)  1 
for  k*0,  and  there  is  no  preferred  direction  of  scattering  (k*l).  Direc¬ 
tional  scattering  does,  however,  appear  in  the  multiple-scatter  terms  (k>2), 
in  view  of  the  directional  character  of  the  correlation  function 
R12 ^1**1 ’12**2^*  cf.  (7.53).  See  [41,  II]  for  details. 

8.6  Extension  to  Include  Explicit  Capillary  Waves 

In  our  preceding  wave' surface  models,  cf.  (7.4b),  we  have  not  ex¬ 
plicitly  distinguished  the  "capillary  wave"  component  (cc)  from  the  "gravity 
wave"  component  Us)  in  the  (single)  continuous  wave  surface,  c6,  cf.  (7.3). 
Although  the  so-called  capillary  term  (when  present)  is  small  vis-a-vis 
the  gravity  component  and  the  separate,  ^oliton  ripples  U5),  it  cannot 
be  neglected  at  small  grazing  angles  and  high  frequencies,  where  the 
"geometrical  acoustic"  or  specular-point  scatter  term  itself  becomes  van¬ 
ishingly  small. 

Using  (7.3)  for  cQ  we  consider  first: 
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A..  The  Coherent  Component  ( k=Q) 

Here  we  have  now  for  (7.19)  the  following  expression  for  the  second- 
moment  function: 


u(° ) /  | 

kT  H- 


=K(V^)g  e  °  °ZCG  (1  +  2iko^o*^G;S  +  «W^C 

k2  *  2 

"  f  f-0 -G^?S  +  +  ^  ^ccs]  •••^|  ART' 


(8.18) 


With  cg,  cc>  Cj  taken  to  be  statistically  independent,  wi  ,h  tg+^c  =  tg 
gaussian*,  so  that  Cg,  Cgx<  etc.  are  likewise  statistically  independent, 
we  see  that  (8.18)  reduces  to 
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(8.19) 


where  now 
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S  *  C.f .  Of  (-;  )  (*  e  02  0  9),  cf.  (7.26).  (8:194 


The  "tilt-factor"  is  given  by  (7.22),  (7.22a),  (7.24).  The 

(01  u-con 

"tilt-factor"  N;u;_.  is  obtained  from 
c-con 


,(0)  . 


c-coh 


<W  >-  * 


(8.20) 


The  other  factors  are  given  as  before:  cf.  (7.25)  for  ART.  Thus,  (7.23) 
is  extended  now  to _ 

?  k2h2rV„2N<0>  +k2nV°) 

M(°)/t0|  ,  ,  \2  »  ,F  .  ,2  koboWsNG-coh  ko°cNt-coh  ( 

%("!•••)  (2aoz)  ART|FlgiG'  e  < 

where  we  can  drop  the  terms  k2c2N^,  k^N^  (<<  1)  in  the  exponent, 

0  5b  0  C  C 

♦Note  that  ;  ,  ;  are  not  individually  gaussian,  strictly,  since  ; 
is  certainly  not  gaussian.  (However,  in  practice  we  may  treat  tg  as  an 
approximately  normal  process,  since  is  small  vs.  cg,  usually,  and.’. 

;g  4  ;GJ 


(8.21) 
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2  «  2 

and  F.  •  F,r,  since  s  Thus,  for  the  coherent  component  of  scat- 
19  /Q\l«  G  9 

tering  is  now  given  by  (8.7),  with  the  exponential  term  of  (8.21), 

2  2  2  9  9 

which  is  practically  replaced  by  exp(-b*k*Cg),  °g  *  °g» 

B;  The  Incoherent  Component  (k«0) 

Our  modification  of  the  preceding  results  [Sec.  7.3B  et  seq.]  starts 
with  (7.29).  We  employ  the  aforementioned  technique  of  splitting  the 
wave-number  spectrum,  discussed  in  Sec.  2.1  above.  We  begin  by  rewriting 
the  exponent  of  (7.29)  as 


boko  G(1  '  PG5  *  bSkoCG{1  ‘  Pg  "  6c*;  5g  +  fc  *  PG;  °g  +  °c  *  °G’ 


.2l2,,2 


2  t  „2  _  „2 


(8.22) 


where  Pg(0,0)  *  1  is  truly  normalized,  but  6^  c (0,0) <  1.  Further 
rewriting  gives 


bokoaG^ 


pg)  *  boko{°g(lH5g>  +  °c(1‘6c)}; 
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(8.23) 


Now  6^,  c^.  are  properly  normalized.  Accordingly,  we  have  for  (7.29) 

F2S  *  e*p  :-boko°G(1‘cG!1  *  <exp  ['t'oko°g(’-69)1)  11  +  kobo°c(1'-pc)  * 

.*  -  c  c  0  0  c 

.  -bnkn0n(l"^)  ?  ?  9 

c  e  009  9  {i  +  k  5  o?p  +  , . 

0  0  6  C 


k2b2o2  <<  1, 
0  0  c 


(8.24) 


since  the  capillary  component  is  very  small,  even  at  high  frequencies: 
kQ  *  83.78  rad*m’*  («20  kHz),  cos  80°  *  0.174,  oc  *  0.5  nm.  This  gives 
(kftb  o  )2  -  [{83. 78) *2* 0.174*  (5*  10’4)]2  -  2.13*  10‘.4«  1,  «  5.33  10’3  at 
100  kHz,  etc. 


The  covariance  functions  Pg,  dc  are  specifically,  from  (2.3), 

k 

°G‘VAr‘T')  =  ^  Wq^'T'^e  ~  '"r  ~ J =  ^  W  (k|T’)e-1-’^  (8.25a) 

G  9  -  9  (20  -kD  9  ( 2tt  ) 

^cter*1')  =  /  wc(£JT)e  ~  — ^2  =  /  Wc(k|T '  )cos  Jc*4r  .  (8.25b) 

-*  (2tt)  kn  2ti 


Equation  (7.31)  is  now  extended  to 


Eq.  (7.31) 


+  ko  N6S-incKS^’T')  ’ 


where  cr„„  and  5.  are  independent  (cf.  remarks  following  (8.18a)  above) 

/n\  c 

and  where  is  given  by  (7.32),  (7.33),  (7.65),  (7.66),  and 


F2g  “  e 


-bVo^d-p') 
o  o  g  g 


o  o 

for  the  essentially  gaussian  "gravity-wave"  component.  [Since  Og  =  og 
here,  p^  =  pg  and  pg(0,0)  =  1.] 

With  (8.26)  as  the  argument  in  {  }  of  (7.34a)  we  may  proceed  as  in 
(7.35),  to  get  now 
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Go-inc' 
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b2  .  N<° 
0  GS 


At  this  point  we  proceed  as  in  Section  (8.1),  using  (8.25b)  to  write 
finally  the  extended  version  of  (8.6): 
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The  coherent  component  is  given  by  (8.21)  in  (8.7)  (with  og  =  oQ),  for  all 
frequencies. 
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Similarly,  we  find  that  the  'low- frequency"  version  of  (8.9)  becomes 
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where  A^Ug)  is  given  by  (8.9a,b),  etc. 

We  can  now  specialize  (8.29)-(8.31)  to  the  various  scatter  geometries, 
as  we  may  wish,  e.g.,  "backscatter,"  forward  scatter  in  the  Snell  direction, 
etc.  The  results  appear  as  the  various  cross-sections  presented  in  Sections 
2,  3  above.  Finally,  for  the  diffraction  terms,  cf.  Sec.  8.5,  the  extensions 
are  directly  made  as  given  in  Section  3.3. 
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